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Abstract. We consider the conformal decomposition of Einstein's constraint 
equations introduced by Lichnerowicz and York, on a compact manifold witli 
boundary. We first develop some technical results for the momentum con- 
straint operator under weak assumptions on the problem data, including gen- 
eralized Korn inequalities on manifolds with boundary not currently in the 
literature. We then consider the Hamiltonian constraint, and using order re- 
lations on appropriate Banach spaces we derive weak solution generalizations 
of known sub- and super-solutions (barriers). We also establish some related 
a priori L°°-bounds on any '^-solution. The barriers are combined with 
variational methods to establish existence of solutions to the Hamiltonian con- 
straint in L°° n W^'^. The result is established under weak assumptions on 
the problem data, and for scalar curvature R having any sign; non-negative R 
requires additional positivity assumptions either on the matter energy density 
or on the trace-free divergence-free part of the extrinsic curvature. Although 
the formulation is different, the result can be viewed as extending the regu- 
larity of the recent result of Maxwell on "rough" CMC solutions in W^''^ for 
k > 3/2 down to L°° n W^'^. The results for the individual constraints are 
then combined to establish existence of non-CMC solutions in W^'P, p > 3 
for the three-metric and in L"?, q = 6p/(3 -|-p) for the extrinsic curvature. 
The result is obtained using fixed-point iteration and compactness arguments 
directly, rather than by building a contraction map. The non-CMC result can 
be viewed as a type of extension of the regularity of the 1996 non-CMC result 
of Isenberg and Moncrief down toW^'^ for p > 3, and extending their result to 
R having any sign. Similarly, the result can also be viewed as type of extension 
of the recent work of Maxwell on rough solutions from the CMC case to the 
non-CMC case. Although our presentation is for 3-manifolds, the results also 
hold in higher dimensions with minor adjustments. The results should also 
extend to other cases such as closed and (fully or partially) open manifolds 
without substantial difficulty. 
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1. Introduction 

In this article, we give an analysis of the coupled Hamiltonian and momentum 
constraints in the Einstein equations on 3-dimcnsional compact manifolds with 
boundary. We consider the equations with matter sources satisfying an energy 
condition implied by the dominant energy condition in the 4-dimensional spacetime; 
the unknowns are a Riemannian three-metric and a two-index symmetric tensor. 
The equations form an under-determined system; therefore, we focus entirely on a 
standard reformulation used in both mathematical and numerical general relativity, 
called the conformal method, introduced by Lichnerowicz and York [371[5H[55] . The 
conformal method assumes that the unknown metric is known up to a scalar field 
called a conformal factor, and also assumes that the trace and a term proportional 
to the trace-free divergence-free part of the two-index symmetric tensor is known, 
leaving as unknown a term proportional to the traceless symmetrized derivative of 
a vector. Therefore, the new unknowns are a scalar and a vector field, transforming 
the original under-determined system for a metric and a symmetric tensor into a 
(potentially) well-posed elliptic system for a scalar and a vector field. See [6] for 
a recent review article. We point out just some of the quite substantial number of 
previous related works, including: the original work on the Lichnerowicz equation 
[57] ; the development of the conformal method [511 [S31 [Ml [5S] ; the initial solution 
theory for the Hamiltonian constraint [13 [33] ; the thin sandwich alternative 
to the conformal method [5l [40] ; the complete classification of CMC initial data 
[29] and the few known non-CMC results [30l [311 IE]; various technical results 
on transverse-traceless tensors and the conformal Killing operator [71 [9] ; the more 
recent development of the conformal thin sandwich formulation [56]; initial data 
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for black holes [H [10] ; initial data for Kerr- like black holes [HI [17] ; and the gluing 
approach to generating initial data [IS] . 

The conformal method gives rise to a coupled nonlinear elliptic system for the 
unknown scalar and vector fields; the trace of the symmetric tensor plays an im- 
portant role: in the case that the trace is constant (referred to as the constant 
mean curvature or CMC case), the two equations decouple, giving rise to the term 
"semi-decoupling decomposition" which is sometimes used to describe the conformal 
method 6 . In this linear equation for the unknown vector can be solved 

first, and then a semi-linear equation for the scalar field can be solved, where a 
coefficient in the nonlinearity depends quadratically on derivatives of the vector 
unknown. Almost all of the previous work on developing a solution theory for the 
constraints has focussed on the conformal decomposition in the CMC case, primar- 
ily in the case of compact manifolds without boundary [6 . A notable exception 
is the non-CMC existence and uniqueness result in Holder-classes for a particular 
physical scenario, which was established in '30J. 

In this article, we extend the solution theory for the individual and coupled 
Hamiltonian and momentum constraints on compact manifolds with boundary in 
three ways: 

(i) Some technical results, including generalized Korn inequalities, are estab- 
lished for the conformal Killing operator on compact Riemannian manifolds 
with boundary, under several different boundary condition assumptions. The 
results, which are not currently in the literature, allow us to establish well- 
posedness of the momentum constraint equation in W^''^{M.) on a manifold 
A4 with boundary, using either variational or Riesz-Schauder methods. The 
assumptions we make on the data using either method are weak enough that 
standard techniques to establish additional regularity are not available. 

(ii) Existence (and in some cases, uniqueness) results are established for weak (or 
rough) CMC solutions to Hamiltonian constraint, for weaker solution spaces 
than appeared previously in [30[ I38j. In particular, using variational meth- 
ods we establish existence of weak solutions to the Hamiltonian constraint 
in L°°{M) n W^''^{M) on compact manifolds with boundary, under assump- 
tions on the data that do not allow for the use of standard techniques to 
establish additional regularity. The variational methods we employ make use 
of (generalized) barriers for the Hamiltonian constraint equation; in ^A.Sl we 
summarize the barriers we use for different values of the Ricci scalar of the 
background metric. We also establish some related a priori L°°-bounds on any 
VF^'^-solution to the Hamiltonian constraint in H4.4\ Although such results 
are standard for semi-linear scalar problems with monotone nonlinearities (see 
for example [32]), our results hold for a class of non-monotone nonlinearities 
that includes the Hamiltonian constraint nonlinearity and appear to be new. 

(in) Existence results are established for non-CMC solutions to the coupled system 
of constraints on compact manifolds with boundary, in the setting of weaker 
(rougher) solutions spaces and for more general physical scenarios than ap- 
peared previously in |30| . In particular, we establish existence of solutions to 
the coupled Hamiltonian and momentum constraints, in W^'P{A4) for p > 3 
for the conformal factor and in W^''^{M-) for q — 6p/{3+p) for the momentum 
vector, on compact manifolds with boundary, with no restrictions on the sign 
of the scalar curvature R. For the case of non-negative R, either the matter 
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energy density or the trace-free divergence-free part of the extrinsic curvature 
must be globally positive. The technical condition on the trace of the extrinsic 
curvature (called the "near-CMC" condition in 6 ) used to produce the cou- 
pled system result in [30] is still present here, although it now involves weaker 
norms (see ij5]). In addition, this condition is only used here to construct a 
global super-solution to the Hamiltonian constraint, and is not used a second 
distinct time as part of the fixed-point argument as was needed in [30! . 

The results above imply that the weakest differentiable solutions of the Einstein 
constraint equations we have found correspond to CMC hypersurfaces with physical 
spatial metric hab and extrinsic curvature kab satisfying 

habeL°°{M)nW^'HM), kab^L^M). (1.1) 

The curvature of such data can be computed in a distributional sense, following [24] . 

There are at least four distinct, but related, motivations for establishing the ex- 
tensions outlined above. First, as outlined in [6], new results for the non-CMC case, 
beyond the case analyzed in [30j , are of great interest in both mathematical and nu- 
merical relativity. Second, there is currently substantial research activity in rough 
solutions to the Einstein evolution equations, which rest on rough/ weak solution 
results for the initial data [34] . Third, the role of boundary conditions and bounded 
domains in the solution and approximation theory is of importance particularly in 
numerical relativity; most existing results are for closed (compact without bound- 
ary), open, or only partially bounded domains. Finally, the approximation theory 
for Petrov-Galerkin-type methods (including finite element, wavelet, spectral, and 
other methods) for the constraints and similar systems previously developed in |28| 
establishes convergence of numerical solutions in very general physical situations, 
but rests on assumptions about the solution theory; the results in the present paper 
help to complete this approximation theory framework. 

An outline of the paper describing the results is as follows. 

In we give a brief outline of the notation used throughout the paper. In 
Sj2l we quickly overview the conformal decomposition, describe the classical strong 
formulation of the resulting coupled elliptic system, and then define weak formula- 
tions of the constraint equations that will allow us to develop solution theories for 
the constraints in the spaces with the weakest possible regularity. Our formulation 
allows for a mix of Dirichlet and Robin boundary conditions for modeling e.g. black 
hole and other physically important scenarios. 

In S|3l we develop some basic technical results for the momentum constraint 
equation on compact manifolds with boundary that we will need later for anal- 
ysis of the Hamiltonian constraint and for analysis of the coupled system. We 
first develop the weak formulation of the momentum constraint for a given scalar 
conformal factor in ^3.11 and then develop some preliminary results related to the 
Korn inequality in ?13.2\ In particular, we establish generalized Korn inequalities for 
the conformal Killing operator on compact manifolds with boundary under several 
boundary condition scenarios; these results do not appear to be in the literature. 
In ij3.3l we use the preliminary results from tj3.2[ together with a variational ar- 
gument, to establish existence and uniqueness of weak solutions to the momentum 
constraint in W^'^{A4) when the Dirichlet part of the boundary is non-empty, with 
assumptions on the data that do not allow for additional regularity in the sense 
described earlier. In §3. 41 we give a second (non- variational) argument for existence 
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and uniqueness again in W^''^{M-) using a Riesz-Schauder (Fredholm alternative) 
argument following. This second argument is more general than the variational ar- 
gument in the sense that the Dirichlet part of the boundary might be empty. While 
the techniques we use in the analysis of the momentum constraint are standard, 
this collection of results of the momentum constraint operator on compact domains 
with boundary apparently are not in the literature. Regularity of solutions to the 
momentum constraint is discussed briefly in §3.51 

In 3H we give a corresponding analysis of the Hamiltonian constraint on compact 
manifolds with boundary. We first develop the weak formulation of the Hamilton- 
ian constraint for a given momentum vector variable in §4.1i and establish some 
preliminary results on local and global barriers (constant sub- and super-solutions) 
for weak solutions in §4.21 The term local means that the barrier depends on the 
given momentum vector variable solution of the momentum constraint equation, 
while global means that the barrier is not local. These barriers are non-trivial ex- 
tensions of those in [29l [30] to nonlinearities with coefficients in W~^'^{A4). We 
also establish some related a priori L°°-bounds on any VF^^^-solution to the Hamil- 
tonian constraint; although such results are standard for semi-linear problems with 
monotone nonlinearities, our results hold for a class of non-monotone nonlinearities 
that includes the Hamiltonian constraint nonlinearity and appear to be new. The 
weak solution barriers are critical to extending the solution theory for the Hamil- 
tonian constraint to the weakest possible setting of W^'^{A4) in §4.5[ and are also 
key to extending the solution theory for the coupled system to weaker spaces and to 
new physical scenarios in |j5l In §4.5) we use the barriers from §4.21 together with 
a variational argument, to establish existence and uniqueness of solutions to the 
Hamiltonian constraint in the weakest possible setting of L°°{A4) r\W^''^{Ai). Due 
to the lack of Gateaux-differentiability of the nonlinearity in W^''^{Ai), the con- 
nection between the energy used for the variational argument and the Hamiltonian 
constraint as its Euler condition is non-trivial, and is established through several 
Lemmas. We note that our arguments allow for scalar curvature to have any sign, 
and the assumptions on the data are such that additional regularity is not possible 
in the sense described earlier. The results for non-negative R require an assumption 
of global positivity either on the trace-free and divergence-free part of the extrin- 
sic curvature or on the matter energy density. Although our problem formulation 
is somewhat different, the result can be viewed as extending the regularity of the 
recent result of Maxwell [35| on "rough" CMC solutions in W''''^{M) for k > 3/2 
down to L°°{A4) nW^''^{Ai). In ij4.6l we give a second (non-variational) argument 
for existence and uniqueness, using the barriers (sub-/super-solution) approach as 
in most of the earlier work [29l [30l [38] . Unlike the case of the momentum con- 
straint in i i3.4l where the non-variational technique allows for the development of a 
solution theory in same weak setting of W^'^{A4) as does the variational method, 
the Hamiltonian constraint barriers approach requires additional regularity beyond 
what the variational approach in §4.51 requires. Regularity of solutions is discussed 
briefly in §47l 

Finally, in ^we use the results for the individual constraints derived earlier to 
establish a new non-CMC result for the coupled system. In §5.11 we establish exis- 
tence of non-CMC solutions to the coupled constraints through flxed-point iteration 
and compactness arguments directly, rather than by using the Contraction Mapping 
Theorem as was done in the original work of Isenberg and Moncrief in [30' . The 
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"near-CMC" assumption on the trace of the extrinsic curvature required in [30j is 
still present here, although it now involves weaker norms (see In addition, the 
condition is only used here to construct a global super-solution to the Hamiltonian 
constraint, and is not used a second distinct time as part of the fixed-point argu- 
ment as was necessary in |30| . If a global super-solution can be constructed without 
the near-CMC assumption, then this new coupled system result would still hold for 
"far-from-CMC" scenarios. The result requires more regularity than that needed 
for the result established in [J2]and [2] for the individual constraints, with solutions 
in W^'P{A4) for p > 3 for the conformal factor and in W'^'''{Ai) for q — 6p/ (3 +p) 
for the momentum vector, but still extends the existing theory for the system in 
two distinct ways. First, although our problem formulation is somewhat different 
(bounded domains with matter), the result can be viewed as a type of extension of 
the 1996 non-CMC result of Isenberg and Moncrief in [30j from the scalar curvature 
i? = — 1 case to R having any sign (with non- negative R requiring the assumption 
that either the matter energy density or the trace-free divergence-free part of the 
extrinsic curvature be globally positive), and to weaker solution spaces. Second, 
again although the problem formulation is different, the result could be viewed as a 
type of extension of the recent rough CMC solution work of Maxwell in [38| to the 
non-CMC case. Although our presentation is for 3-manifolds, the results hold in 
higher spatial dimensions with minor adjustments, and the techniques we employ 
should extend to other cases such as closed and (fully or partially) open manifolds 
through the use of tools such as weighted Sobolev spaces. 
We summarize our results in S|6l 

1.1. Notation and conventions. Let {A4,hab) be a a Riemannian manifold, 
where A4 is a 3-dimensional, smooth, compact manifold with non-empty boundary 
dA4, and hab is a metric on Ai, that is, a symmetric, positive definite, covari- 
ant, two-index tensor on M with all components in a smooth coordinate system 
having two continuous derivatives. Latin indices denote abstract indices as they are 
defined in [50], §2.4. The metric defines an inner product on T^Ad, the vector space 
tangent to Ai at the point x ^ Ai. Denote by ft,"'' the inverse of the metric tensor 
hab, that is, hach'"^ = 6a'', where 6a'' ■ T^Ai T^Ai is the identity map. We use 
the convention that repeated indices, one upper-index and one sub-index, denote 
contraction. Let Va be the Levi-Civita connection associated with the metric hab, 
that is, the unique torsion-free connection satisfying Wahbc — 0. Let Rabc^ be the 
Riemann tensor of the connection Vq , where the sign convention used in this article 
is (VaVfc — VhVa)t;c := Rabc'^Vd- Dcnotc by Rab ■— Racb'^ the Ricci tensor and by 
R :— Rabh"'^ the Ricci curvature scalar of this connection. (Only in §2.1l we modify 
the notation for the connection Vq.) 

Indices on tensors will be raised and lowered with ft,"'' and hab, respectively. 
For example, given the tensor s"''c we denote Sabc = haaih^b^ s°'^''^c, and s"''^ = 
j^cci gO-b^^ . notice that the order of the indices is important in the case that the 
tensor Sabc or s"'"^ is not symmetric. We say that a tensor is an n-index tensor 
iff it can be transformed into a tensor Sai-.-a^ by lowering appropriate indices. We 
denote by C°° {A4 , n) the set of all smooth rt-index tensor fields on M . Given 
an arbitrary tensor s^i ' ^^bj^...;,^^ , which is an n -I- m-index tensor, we define its 
magnitude at any point x ^ Ai as the real-valued function given by 
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Integration on A4 is performed with the volume element dx associated to the metric 
hab- A norm of an arbitrary smooth tensor field s'^^"""6i - b„ on A4 can be defined 
for any 1 ^ p < cxd and for p = oo respectively using (|1.2p as follows, 



\s\Pdx 

M 



1/p 

\\s\\oc := ess sup \s\. (1.3) 

xeM 



We introduce the Lebesgue spaces L'P{A4,n), for 1 ^ p ^ oo, oi n-index tensor 
valued fields as the completion of C°°{A4, n) under the norm in Eq. ()1.3|) . and this 
norm is called the L^-norm. The Lebesgue spaces LP{A4,n) are Banach spaces; 
they are separable when I ^ p < oo and reflexive when 1 < p < oo. For the case 
p = 2 the spaces (7W , n) form a Hilbert space with the inner product and norm 
given by 



(s,r):=/ s<,,...,„r'^--"dx, ||.s|| := ^(i:^) = ||s||2. (1.4) 
Covariant derivatives of tensor fields are denoted as 

s ■— vbi,--- ,b,„s .— Vbi ■ ■ ■ yb„,s , 

where the super-script m indicates the total number of derivatives, which plays the 
role of the number \a\ for a multi-index a, in the multi-index notation used in the 
PDE literature. Using again Eq. (|1.2p introduce the real-valued function 

as the starting point for defining L^-type norms involving derivatives. One such 
norm in the vector space C°° {Ai , n) is given for any non-negative integer k and a 
real number p with 1 ^ p < oo, and separately for p = oo, as follows 



k 



1/p 



s|U,p:= > l|V™sE , ||s||fc,oo := max IIV-sIU. (1.5) 



m=0 

We introduce the Sobolev spaces W'''P{A4,n) of n-index tensor valued fields as 
the completion of C°° {A4 , n) under the norm in Eq. (jl.Sp , and this norm is called 
the W^'P-norm. The Sobolev spaces VF'^'P(7W,n) are Banach spaces; being based 
on LP(A4,n), they are separable when 1 ^ p < oo and refiexive when 1 < p < oo. 
For the case p = 2 the spaces W^'^{M,n) form a Hilbert space with the inner 
product and norm given by 

fe 

{s, r)u (^"^' ^"^)' = v/(^, (1-6) 

m— 

where we have introduced the notation 

(V"s,V"r):= / (V;„---V;,„s„,...,J(V''i---V'''"r"--'^")dx. 

J M 

Therefore we have that LP{M,n) = W°'P{M,n) and \\s\\p = \\s\\o,p. These defini- 
tions follow [27] for the case of scalar fields and [45] for the case of arbitrary tensor 
fields. These definitions can also be extended, using appropriate partitions of the 
unity and Fourier transforms, from non- negative integers k to real numbers s. For 
example see [i5] . 
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In this article we are mainly concerned with spaces of scalar-valued fields and 
vector- valued fields on , so we introduce the following notation for these special 
cases, 

C"°{M,0), C°°(A4,I), 

LP(A^,0), LP -.^ LP{M,l), 

M^'='f(7W,0), W'''" -.^W'^^PiM,!). 

However, we will not suppress the manifold from the notation of these spaces when 
this information is important in a given situation. In a similar way, we use both 
notations w and w'^ to denote a vector field. We consider in this article that the 
boundary dM of M can be divided in the following two, possible different, ways 
as follows, 

dM = dMDlidMN, OMondMN = 9, (1.7) 
dM^dMiDlidMrn, dMnndMr^ = 9. (1.8) 
Introduce the trace operators 

tro : W'P ^ W'-p-PidMn,0), tr^ : W'P ^ W'-p^P{dMn, 1), 

for s > 1/p, which are the continuous extensions to Sobolcv spaces of the operators 
defined on smooth fields given by tr^c/) :— <p\dMn, and tr joW := w\gM]Q, see for 
example gT]. Both spaces W'-^'P{dMD,0) and Wp^PidMnA) are Banach 
spaces and we denote their norms as ||(/>||5_i ^ and p jp, respectively. In 

the particular case p = 2 these spaces become Hilbert spaces and we denote their 
inner product as {(j), <t>)s-^ d and {w,w)^_i j^. We will be mainly concerned with 
the case s = 1, p = 2, and in this case we denote their inner product and norms as 
follows 

(0,0)o, = (0, (^)o''^ {w,W)iD, \\w\\i, = {w,w)]l'^, 

for aU 0, ^ G W^'^idMo, 0) and all w,w€ wi^^idMn, !)■ 

In an analogous way, introduce the trace operators trjv and tr^v and the spaces 

W''~p-p{Mn,0) and W^p-PiMMA)- We use the notation 

W^'P {0 e W^-P : trc0 = 0, tr^^(/) = 0}, W^'^ := {<f) £ W^^p : tr^^ = 0}, 

Wl^P := {w e W^'P : trjow = 0, tr^io = 0}, W],:^ := {w e W^'P : trjow = 0}, 

for the function spaces. The dual spaces of some Sobolev spaces will be denoted as 
follows: 

W-'^'P := [W'^^P']*, W^'^'P := [W^Y, W^'^'P := [VF^]*, 

W-'^'P := [W'^-P']*, Wo'^'^P := [W'^'"']*, W'^^^p [W''/]* , 

where we denote by p' the conjugate of p in the sense ^ + ^ = 1- These are Banach 
spaces with the norm 

l|s*||_fc,p:= sup ' I ^' , 

Q^r^W-p' (M,n) If llfe,p' 

where s* G W~'''P{A4,n), and the asterisk is introduced to emphasize that s* is 
a linear and bounded map s* : W^'P (A4,n) R. The product of an element 
in M e L°°(A^,0) by an element in s* G W~^'P{M,n), with I < p < oo will be 
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denoted by (su)*. Such element is a weh-defined functional in W~^'P{Ai,n). The 
proof of this statement is given in the Appendix using appropriate Gelfand triple 
structures. 

We will need order structure in some Sobolev spaces. See the Appendix for a 
review on ordered Banach spaces, where we explain the particular notation from 
this field we use in this article, and where we define the main order cones needed 
in this article: L°^, L^, and W^'^. We use both notations {u — v) Q X+ and u ^ v 
to state that the element w in a Banach space X is bigger than or equal to another 
element v in that space. The former notation specifies which Banach space the 
elements belong to and also which order cone is used in that particular Banach 
space; while both pieces of information are not explicitly displayed in the latter 
notation. We use the notation u ^ v when there is no ambiguity, otherwise we use 
the notation {u — v) G X+. We also write —{u — v) G X+ to denote u ^ v. 

Given Banach spaces X, Y and an operator A : Da C X ^ Ra C F, we denote 
by Da and Ra the domain and range of A, respectively, while Na denotes the null 
space of A. 

We recall the generalized Holder inequality (see [35], page 904), which is 
used in several places in this article, and says that given n-index tensor fields 
Si G L^''(A^, n) with i = 1, • • ■ ,k for fc G N, the pointwise product is well-defined 
a.e. in Ai, the tensor field si • ■ • G U'{M.,nk), where p = X^iLi ^-^d the 
following estimate holds, 

||si---Sfe||p ||si||pi ■•■lls/cllpfc. 

One last comment on the notation: Given a Banach space X, the elements x, 
x&X denote different elements. This notation usually appears in equations written 
in weak form, where the element x denotes the trial function and the element x 
denotes the test function. 

2. The constraint equations 

We give a quick overview of the conformal decomposition, describe the clas- 
sical strong formulation of the resulting coupled elliptic system, and then define 
weak formulations of the constraint equations that will allow us to develop solu- 
tion theories for the constraints in the spaces with the weakest possible regularity. 
Our formulation allows for a mix of Dirichlet and Robin boundary conditions for 
modeling e.g. black hole and other physically important scenarios. 

2.1. The conformal decomposition method. Let (Af, g^^) be a smooth 4-di- 
mensional spacetime, that is, M is a 4-dimensional, smooth manifold, and g^i, is a 
smooth, Lorentzian metric on M with signature (— , +, +, +). Let V^x be the Levi- 
Civita connection associated with the metric gr^j/, that is, the unique torsion- free 
connection satisfying Vo-g^j/ = 0. The Einstein equation is 

where G^jy = i?^,^ — Rg^^/'i. is the Einstein tensor, T^^ is the stress-energy tensor, 
and K, = SttG/c'*, with G the gravitation constant and c the speed of light. The 
Ricci tensor is Rf^^, — R^cru'^ and R — R^^g^^" is the Ricci scalar, where g^^" is the 
inverse of 5^^, that is gfj.ag'"^ — . The Ricci tensor is defined as a contraction 
of the Riemann tensor R^^a'^Wp ~ (V^Vi, — Vi/V^)u;(j, where is any 1-form 
on M. The stress-energy tensor T^y is assumed to be symmetric, to satisfy the 
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integrability condition V^T^'' = 0, and the dominant energy condition, that 
is, the vector —T^'^Vu is timehke and future-directed, where is any timehke and 
future-directed vector field (see [SD], page 219). In this section Greek indices /i, v, cr, 
p denote abstract spacetime indices, that is, tensorial character on the 4-dimensional 
manifold M . They are raised and lowered with g'^'^ and g^^, respectively. Later on 
Latin indices a, 6, c, d will denote tensorial character on a 3-diniensional manifold. 

The map i : M — > R is a time function iff the function t is differentiable and 
the vector field —V^t is a timelike, future-directed vector field on M . Introduce 
the hypersurface M :— {x € M : t{x) = 0}, and denote by ti^ the unit 1-form 
orthogonal to vector fields tangent to M.. By definition of the 1-form has 
the form rt^ = — aV^t, where a is a positive function such that n^n,ygf^'^ — — 1, 
which is called the lapse function. Since the lapse function is positive, the vector 
field n'^ is future-directed. Let /i^^ and /c^^ be the first and second fundamental 
forms of the hypersurface M, that is, 

The Einstein constraint equations on A4 are given by 

It is a straightforward albeit long computation to express these equations involving 
tensors on M as an equation involving tensors on A4. The result is the following 
equations, 

''R + P -kahk"'' -2kp^0, (2.1) 
D^k - bbk"^ + Kf = 0, (2.2) 

where tensors hab^ kab, T ^'^^ p on a 3-dimensional manifold are the pull-back 
on M. of the tensors h^^v, fc^^, and p on the 4-dimensional manifold M. We 
have introduced the energy density p :— n^rt^T^'' and the momentum current 
density := —h^\,n„T'^"' . We have denoted by Da the Levi-Civita connection 
associated to hab, so {A4,hab) is a 3-dimensional Riemannian manifold, with hab 
having signature (4-, 4-, -I-), and we use the notation for the inverse of the metric 
hab- Indices have been raised and lowered with ft,"'' and hab, respectively. We have 
also denoted by the Ricci scalar of curvature of the metric hab- Finally, recall 
that the constraint Eqs. (|2.ip - (|2.2p are indeed equations on hab and kab due to the 
matter fields satisfying the energy condition — + jaj"" < 0, which is implied by 
the dominant energy condition on the stress-energy tensor T^"^ in spacetime. 

Let be a positive scalar field on J\4, and decompose the extrinsic curvature 
tensor kab — Sab + habT/S, where f := kabh"'^ is the trace and then Sab is the 
traceless part of the extrinsic curvature tensor. Then, introduce now the following 
conformal rescaling: 

hab <l>^ hab, f=:T, (2.3) 

f^--<p-'"r, (2.4) 

We have introduced the Riemannian metric hab on the 3-dimensional manifold A4 , 
which determines the Levi-Civita connection Da, and so we have that Dah^c — 0. 
We have also introduced the symmetric, traceless tensor Sab, and the non-physical 
matter sources j" and p. The different powers of the conformal rescaling above are 
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carefully chosen so that the constraint Eqs. (|2.1|) - (|2.2p transform into the following 
equations 

-8A0 + '^R(j)+^ r'^cjf' - Sabs"^ (j)^'^ - 2kp (f>-'^ = 0, (2.5) 

-Dfcs"'' + + Kj"^ = 0, (2.6) 

where in equation above, and from now on, indices of unhatted fields are raised 
and lowered with ft,"'' and hab respectively. We have also introduced the Laplace- 
Beltrami operator with respect to the metric hab, acting on smooth scalar fields; 
it is defined as follows 

A0 := h'^^'DM. 

Eqs. (I2.5p - (I2.6|) can be obtained by a long, but otherwise straightforward compu- 
tation. In order to perform this calculation it is useful to recall that both Da and 
Da are connections on the manifold A^, and so they differ on a tensor field Cab'^: 
which can be computed explicitly in terms of 0, and has the form 

Cab' = 45(a^i?6) ln(0) - 2habh''''DdH^). 

We remark that the power four on the rescaling of the metric hab and M being 
3-dimensional imply that = (f>~^{^R(j) — 8A(j>), and for any other power in the 
rescaling, terms proportional to h'^^{Da4>){Db4>) / (jy^ appear in the transformation. 
Similar reasons force the power negative ten on the rescaling of the tensor s""^ 
and so terms proportional to {Da(j))/4' cancel out in Eq. (|2.6|) . Finally, the 
ratio between the conformal rescaling powers of p and j" is chosen such that the 
inequality — -I- habj""]^ < imphes the inequality — + habj°'f' < 0. 

There is one more step to convert the original constraint Eq. (I2.1|) - (l2.2p into a 
determined elliptic system of equations. This step is the following: Decompose the 
symmetric, traceless tensor Sab into a divergence-free part aab, and the symmetrized 
and traceless gradient of a vector, that is, s"''' =: a""^ + [CwY^ , where DaCr"^^ = and 
we have introduced the conformal Killing operator C acting on smooth vector 
fields and defined as follows 

[Cwf^ := D^w^ + D^w" - - {Drw") (2.7) 

3 

Therefore, the constraint Eqs. (|2.ip - (|2.2p are transformed by the conformal rescaling 
into the following equations 

-8A(/. + 'i?</> + ^ T^,/)^ - {<Jab + {Cw)ab) (fT°' + {Cwf'') " 2np<j>~^ = 0, (2.8) 

-Db{CwY^ + \<I)'^D''T + Kf = Q. (2.9) 

In the next section we interpret these equations above as partial differential equa- 
tions for the scalar field (j) and the vector field w°, while the rest of the fields 
are considered as given fields. Given a solution (f) and of Eqs. (|2.8p - (|2.9p . the 
physical metric hab and extrinsic curvature fc"'' of the hypcrsurface M. are given by 

hab = ct^^ hab, = [a''" + [CwT^] +\r^T ft'^^ 

while the matter fields are given by Eq (|2.4p . 
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2.2. Classical formulation. Beginning in tliis section, we will change the notation 
slightly from the classical notation used to introduce the conformal method in §2.11 
In particular, the Levi-Civita connection of the metric hat on the 3-dimensional 
manifold Ai will denoted by Vq rather than Da, and the Ricci scalar of hat will 
be denoted by R instead of '^R. This change will simplify the presentation in the 
remainder of the paper. 

Let {A4,h) be a 3-diniensional Riemannian manifold, where is a smooth, 
compact manifold with smooth boundary dAi, and h € C°°{A4,2) is a positive 
definite metric. Let L : C°° C°° and X : ^ be the Laplace-Beltrami 
and momentum operators, respectively, with actions on a scalar field (p £ C°° and 
a vector field w G given by 

L(j) := ~A(j), (2.10) 
(Ew)" := -VbiCw^'', (2.11) 

where Acj) :— VaV^i?!), and C denotes the conformal Killing operator defined in 
Eq. (|2.7p . We will also use the index- free notation ILw and £w. Assume that dA4 
is divided according to Eqs. (|1.7p - (|1.8p . Then, consider boundary conditions for the 
scalar equation of Dirichlet type on dAi u and of Robin type on dAi^. Analogously, 
consider boundary conditions for the vector equation of Dirichlet type on dA4 ed and 
of Robin type on dAijN- The conditions 

are needed to simplify the proofs regarding the regularity at the intersection points 
of the two types of boundaries of solutions of elliptic equations with Dirichlet-Robin 
boundary conditions. In what follows we include the cases given by dA4r> — 9 or 
OMn = 0, and dMiD = or dAdi^ = 0. 

The freely specifiable functions of the problem are a scalar function r, inter- 
preted as the trace of the physical extrinsic curvature; a symmetric, traceless, and 
divergence-free, contravariant, two-index tensor cr; the non-physical energy density 
p and the non-physical momentum current density vector j subject to the require- 
ment —p^ + habj'^j'' < 0. The term non-physical refers here to a conformal rescaled 
field, while physical refers to a conformally non-rescaled field. The requirement 
on p and j mentioned above and the particular conformal rescaling used in the 
semi-decoupling decomposition imply that the same inequality is satisfied by the 
physical energy and momentum current densities. Introduce the nonlinear opera- 
tors F : C°° X ^ and F : C°° ^ given by 

F{(j}, w) := ar(t>^ + aR(j) - apcj)-^ - a^0~^ (2.12) 
JF((/)) := 6,</)S + 6,, (2.13) 

where the coefficient functions are defined as follows 

«P := Jp, (2.14) 

^V°T, (2.15) 

Notice that the scalar coefficients a^, a,w, and Op are non-negative, while there is 
no sign restriction on a^. 
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The classical Dirichlet-Robin boundary value formulation for the semi- 
decoupHng Einstein constraint equations is the foUowing: Given the freely specifi- 
able smooth fields r, cr, p, and j in A4 and the smooth Dirichlet boundary data 
4>o on dMo and ivjo on dMio, and smooth Robin boundary data 0]v, K, scalar 
fields on dAi^ and tfw, IK, a vector and a two-index tensor fields on d/AjM, find 
a scalar field cj) and a vector field w in Ai solution of the system 

[ n ■ Wcf) + K4' = (j)N on OMn, 

( W= Wr, on dAiri, 

Lw+F{^) = OmM, , ' (2.17) 

1^ n ■ Vw + JKw= Wjpf on oMm, 

where IKw denotes the vector field {IKw}^ ~ K'^'^wi,, and n Ww denotes the vector 

The classical formulation has been done on spaces of smooth fields, which are 
not complete spaces under any known norm defined on them. This is inconvenient 
for finding solutions to PDE, because these solutions are usually found as limits 
of appropriate approximations. If a normed vector space is not complete, then a 
Cauchy sequence may not converge. In the next section we introduce the weak 
formulation of the equations above, where we rewrite the classical formulation in 
appropriate normed vector spaces which are also complete. 



2.3. Weak formulation. We present the weak formulation associated with the 
classical formulation with Eqs. (I2.16|) - (I2.17|) . We introduce one of the weakest 
forms of the constraint equations, that is, we assume the weakest regularity of the 
equation coefficients such that the equation itself is well-defined. We will be able to 
obtain existence and uniqueness results for the momentum constraint using either 
variational methods in §3.3fo r Riesz-Schauder methods in §3.41 We will also be able 
to obtain existence (and when possible, uniqueness) for the Hamiltonian constraint 
in M.SI using variational methods in this weakest setting. However, the barrier- 
based existence and uniqueness results for the Hamiltonian constraint equation 
in §4.6i and the compactness argument in §5.11 giving existence for the coupled 
system of constraints, require higher regularity on the equation coefficients, but 
we still obtain some non-CMC results for the coupled system in weaker settings 
and in more general physical situations than have been previously obtained. These 
additional assumptions are clearly stated in those sections. 

Let {A4,h) be a 3-dimensional Riemannian manifold, where is a smooth, 
compact manifold with Lipschitz boundary dA4, and h G C^{A4,2) is a positive 
definite metric. Introduce the bilinear forms 

ar^-.W^'^ xW^-^ ^R, ai(0,^) (V0,V0) + (i^tr^.(/.,tr^^)iv, (2.18) 
ajL : W^''^ X W^''^ M, ai(u;, w) := [Cw,Cw) + {KtrM^w,trM^w)iN, (2.19) 

where the Robin scalar field K £ L°°(9A^jv, 0) and the two-index, symmetric tensor 
field K e L°^{dMM, 2) satisfy the bounds 

k Wu^m ^ {Ktr^4>, tr^0)jv, V0 e (2.20) 
K|[trwHlw «S (^trww,trww)w, Vw e (2.21) 
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with k and K being non-negative constants. In this Section set the number p — 12/5, 
and then fix the source functions 

TGL^ p*(.W-]:P, aeLP{M,2), j*eW^''P, (2.22) 

where a is symmetric, traceless and divergence-free in weak sense, that is, it satisfies 
{a,Cu}_) = for all lj G Wq'^. The asterisk on the matter fields is to emphasize 

that they are elements of spaces of linear functionals. That is, p* : W^'^ — > R is 
linear and bounded, and an analogous definition holds for j* . In the Appendix it 
is shown that the spaces VK^"' C = [L^]* C VF^^'^, with 1 < g < oo, form a 
Gelfand triple, so given any element p* e Wj^^'^ there exists a sequence {pn} C 
such that 

p* (ip) := \im (pn,(p), VipeW^-P. 

— n — >oo — — 

An analogous statement holds for j*. We say that the matter fields p* and j* 
satisfy the energy condition in weak sense iff there exist sequences {pn} C 
and C such that 

[pl-3n-jn]^Ll VneN. (2.23) 

(We have required p* G W^l:^ in Eq. (P?^ instead of p* G 1^^+^^^ because j* 
must belong to W^'^ and Eq. (|2.23p must hold.) Given any function t g L^, then 
it is known that (Vr)* e Wq where the asterisk, we repeat for the last time, 
is added only to reinforce the idea that (Vr)* is a linear functional on elements 
in Wq'P , and it is not meant to indicate the adjoint operator of V. We assume 
here that (Vr)* G W^'^, which is indeed an extra assumption due to W"^^'^ C 
Wq this assumption is needed because the functional (6°)* '■— (2/3) (V^r)* is 
the source in the momentum constraint equation, which requires this particular 
type of boundary conditions. The assumptions above on r and a imply that for 
every w G W^'^ the functions Ur and a„ belong to L^/^. The assumption on the 
background metric implies that Ur is a continuous function on Ai. These three 
functions define the elements a*, and aJi, in the space W^^'^ as follows 

a*{i£) := {ar,(£), a*j^{(£) := [aR.f), al,{f) [a^.f), G W'^'^ . (2.24) 

The proof that these functionals are well-defined is based on Holder inequality, for 
example, consider the functional a*, then 

\{ar,i£)\ ||a^||| ll^lle Cs ||a^||| ||^||i,2- 

These functionals above belong to a particular class of elements in W^^''^, while 
the functionals a* := (k/4)/9* and b* ;= nj* are not restricted to such a particular 
form, and they can be any element in Wjj^'^ and W^^'^, respectively compatible 
with the energy condition. Given any two functions 0i, (f>2 G with 0i ^02, 
define the interval 

[</)!, 02] {0 G : </)! s: s: 02} c 

which is a closed, bounded set in L°°. Assume 0i > 0, and then introduce the 
nonlinear operators 

If ■■ [01,02] C X W^'P ^ W^''\ : [01,02] CL^^ W^'^^, 
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Mcf,, w) := (a,,/)5)* + {anc^y - {upcj)-^)* - (a^,0-")*, (2.25) 
/^(<^) :=(6,06)* + 6*, (2.26) 

where the product of an element (f) £ L°° by an element in a* e W^^'^^ with 
1 < g < cxo, denoted by (a<^)*, is a well-defined element in W~^''^. The proof of this 
statement is given in the Appendix using appropriate Gelfand triple structures. 
The functionals fp and f ^ are the generalizations of the functionals F and F 
defined in Eq. (|2Tl2| - ((2T5| . We remark that the operators defined in Eqs. ([2?25ll - 
p.26p are continuous but not Gateaux differentiable. They have Gateaux derivatives 
only along directions in L°°, not on the whole space L^. This fact introduces 
some technical complexity with the use of variational methods for the individual 
Hamiltonian and momentum constraints (see §3.31 and i i4.5l) . Recall that the trace 
operators 

tro : W^^^ wi^^{dMn,0), tr^ : W^^" -> W^'^'idMn, 1), 

satisfy the following property: given any element (j>]j G W^'^(dMi:,,0), there exists 
an element e W^''^ such that tr^i^ij = </)£,; analogously, given a boundary data 
element Wjd G Wp^'^{dMiD,l), there exists an element Wjo & W^'^ such that 
tr loWjo — Wjo- The elements (pn and Wjo are called here extensions of 4>d and Wjd, 
respectively. They are not uniquely determined by the boundary data. 

The weak Dirichlet-Robin boundary value formulation which is associated 
with Eqs. (|2.16p - (|2.17|) is the following: Fix Dirichlet boundary data 

0<ess inf $d L'^idMD,0) r]W^-^{dMD,0), Wid eWV^P^dMnA), 

dM]j 

with extensions inigMo 0d ^ G W^'^ and Wjo G W^'^, respectively. Choose 
the extension function 0d as a harmonic extension of the Dirichlet boundary data 
using the Laplace-Beltrami operator on A4: 

= in M, troCpD = 0d > on dMo, trjv^D = inf 0o > on OMn- 

dMo 

The maximum principle for the Laplace-Beltrami operator implies that < 
infa_A4 ^ (j)o{x), a.e. in A4. Fix Robin boundary data functionals 

G W~i'^{dMM,0), wl, G W-i'P{dM,^, 1); 

Given the extension c/jjj of the Dirichlet data (f)r> chosen above, fix any two functions 
(pi, (p2 G L°° n W^'^, with the property that < (pi ^ (p2 and satisfying (p^ G 
[01,02] n W^''^. Introduce the non-principal part operators including the Robin 
boundary conditions, 

/: [01,02] CL^ X W''P^W^^^\ /(0,i«)(^) :=/^(0,u;)((^)-0^(tr^^), 

(2.27) 

/ : [01, 02] CL^^ W^^'P, /(0)(a;) := /^(0)(u;) - wl,{trj^u), 

(2.28) 

where fp and fjp are given by Eqs. (|2.25|) - (|2.26|) . Introduce the affine spaces A^-^ 
and A^'P, which include the Dirichlet boundary conditions, as follows, 

^1-2 := 03 + _ 1^ ^ ^1,2 . _ 0^ £ w^f}, (2.29) 
A^'P := wn + W]iP := {w G W^'P : w - Wn E T^lj^}. (2.30) 
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Then, find elements 4> ^ [4>i, 4>2\ H A^'^ and to G A^^'^ solutions of 

a^{c^,^) + f{<P,w)[^) = y^eWl^\ (2.31) 

a^{w,ui)+f{(l)){ui)^Q 'iuieW]iJ''. (2.32) 

It will be convenient later on to express Eqs. ()2.31|) - ()2.32p in terms of operators 
instead of bilinear forms. Introduce the operators 

Ai^-.W^'^ ^W^^-\ A^4>{if):^a^{4>,<£), (2.33) 

: W^'P ^ W-^^^P, A^w[u£) a^{w,ui). (2.34) 

Also recall that given any (p G [0i,02] and w e W^'^ then f{(j),w) G W^^'^ and 
f{(f>) G W]^'^. Hence, the Eqs. (|2.3ip - (|2.32p written in terms of operators is the 
following: Find elements G [4>i, (j)2] H A^'^ and w G A^'^ solutions of 

AL^ + fi<l>,w) = 0, (2.35) 
A,LW + fi(f>)=0. (2.36) 

Lemma 1. Every smooth solution (p, w of the classical problem with Eqs. \2.16\) - 
\2.11^ is also a solution of the weak problem with Eqs. ll2.31\) - [2.32\) . 

Proof. (LemmaUi) Given any smooth fields (j), w solutions of Eqs. (j2.16|l - (j2.17|) . 
then the proof consists in multiplying these equations by test functions (p G W^^' 

and tJ G W],f , respectively, and then integrating by parts. In the case of Eq. ()2.16p 
one gets 

(-Ac/.,^) + (F(</), «;),<£) =0. (2.37) 
The first term on the left hand side can be rewritten as follows, 

(-A0,^) = (V0,V^) - (tr^.(n• V0),tr^^)^ 

= (V0,V^) + ([Xtrjv0-<^iv],trjv(£)^ 

= a£(0,^) - {(f) N, tr 1^ If) 1^, 

which holds for all ip G ^ where in the second line we introduce the Robin 
boundary condition on OMn, and note the integral on dAdo vanishes because the 
test function if vanishes on this part of the boundary, in the third line we introduce 
the definition of the bilinear form aj^. Now, replace this expression into Eq. (|2.37p 
and one obtains Eq. (|2.3ip . Finally, since is a solution of the classical problem, 
it can be written as (/) = (j)o + for some smooth extension (jjo of the boundary 
data (j)o, therefore G A^''^. (Again, (j>o can be constructed e.g. by harmonic 
extension.) In the case of Eq. (|2.17p one gets 

(-V • iCw),uj_) + (F(0), w) = 0. (2.38) 

The first term on the left hand side can be rewritten as follows, 

= aiL{w,uj_) - {wm,tr]MUj) jj^, 

which holds for all a; G W]if , where the first term in the second line comes from the 
symmetries of C, and the second term in that line comes from the Robin boundary 
conditions; the definition of is used to obtain the third line. Now, replace this 
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expression into Eq. p.38p and one obtains Eq. ()2.32|) . Finally, since to is a solution 
of classical problem, it can be written as w = Wjo + uj for some smooth extension 
W]D of the boundary data therefore lo e A^'^. □ 
Let us recall here that the space W]^"^ is an ordered Banach space with order 
cone Wjjj^ defined as follows: 

■■= {(p e : ^ a.e. in M}. 

The order relation is then (/)^(/)iff0 — W^_d+- In the Appendix we discuss the 
main properties of ordered Banach spaces. In particular, we show that the order 
structure implied by W]^'^ can be translated to the dual space W^^''^ as follows, 

Wn+^ := {r e W-''^ : 0*(0) ^0 V0 S }. 

Given two ordered Banach spaces X, X+ and Y, an operator A : Da C X ^ Y 
satisfies the maximum principle iff for every elements u, v E Da such that 
Au — Av G y+ it holds that u ~ v G X^ . In the particular case that the operator 
{A, Da) is linear, then it satisfies the maximum principle iff for every element 
u e X+ such that Au G 1+ it holds that u G X+. If an operator A satisfies the 
maximum principle and is invertiblc, then the inverse is a monotone increasing 
operator, a result shown in the Appendix. This last property is useful to solve 
nonlinear equations of the form Au — f{u), in the case that there exist sub- and 
super-solutions to that equation (see below for the definition). In this case there 
is a well-known existence proof technique that works for many equations of this 
type, and has been one of the main techniques used previously for the Hamiltonian 
constraint [51 [551 1301 [3H| ■ While we will exploit the fact that the construction of 
sub- and super-solutions can be done in a very weak setting, the use of the existence 
proof based directly on barriers requires additional regularity beyond what is needed 
for the barrier construction. This additional regularity assumption can be avoided 
by combining barriers with variational techniques, which we do in §4.51 

The following properties are of interest to us below. Firstly, in the Appendix 
we review results from the literature showing that the operator Aj^ defined above 
satisfies a maximum principle. Secondly, we can show that there exist sub- and 
super-solutions to Eq. (|2.35p . Given any function u G W^'^ , introduce the notation 

:= ess max{u, 0}, := —ess minju, 0}. 

An element cj)^ G W^'^ is called a sub-solution of Eq. (I2.35P iff the function 
satisfies the inequalities 

- cb-)- G W^o'^ and - + f{cb-,w)] G (2-39) 

An element (f>+ G W^''^ is called a super-solution of Eq. (12.3511 iff the scalar 
function (/)_|_ satisfies the inequalities 

- 0+)+ e W^A'' and [A^cf>+ + f{<j)+, w)] G W^l:\ (2.40) 

The sub and super-solutions of Eq. (|2.35p may depend on the choice of the vector 
field w that appears in the functional . A sub-solution (/)_ of Eq. (|2.35p is called 
global iff Eq. (I2.39P holds for every vector field G W^'^ solution of the momen- 
tum constraint Eq. (|2.36p with any source function satisfying (0— </)_) G Wjj^, and 
it is called local iff it is not global. Analogous definitions are introduced for super- 
solutions. While it will be sufficient to derive only local sub- and super-solutions to 
produce the existence and uniqueness results for the Hamiltonian constraint using 
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variational methods in §4.51 and using barrier methods in §4.61 proving results for 
the coupled system rests critically on deriving global sub- and super-solutions for 
this coupled system; we come back to this in fj5] 

3. The momentum constraint 

In this section we fix a particular scalar function cf) G L°° and consider the mo- 
mentum constraint equation (|2.36p for the vector valued function w G W^'^. The 
result is a linear elliptic system of equations for this variable w. We first develop 
the weak formulation of the momentum constraint more precisely in H3.1[ In H3.2\ 
we establish generalized Korn inequalities for the conformal Killing operator on 
compact manifolds with boundary under several boundary condition scenarios; the 
results do not appear to be in the literature. We then briefly summarize here the 
main ideas for solving the Dirichlet-Robin problems for the momentum constraint 
equation, for an appropriately given 0. We use two different methods, namely vari- 
ational methods [33l EHl [57] , and Riesz-Schauder theory for compact operators [51] . 
Both methods yield essentially the same results, since the momentum constraint 
equation is linear in the variable w. 

The variational approach is taken in ! j3.3l when the Dirichlet part of the boundary 
is non-empty, giving existence and uniqueness of weak solutions to the momentum 
constraint in W^'^. The weak assumptions on the data do not allow for the use 
of standard techniques to establish additional regularity. While the variational 
approach has no real advantage over Riesz-Schauder theory for the momentum 
constraint, it will give us some insight in its use for the Hamiltonian constraint, 
for which it will be critical. In addition, some of the supporting results are of 
interest in their own right, so we include the analysis using variational methods 
here along with the Riesz-Schauder arguments. In §3.41 we establish existence and 
uniqueness of solutions to the Dirichlet-Robin problem for the momentum con- 
straint equation using Riesz-Schauder theory for compact operators. The literature 
on Riesz-Schauder theory for systems of elliptic equations is not so clearly presented 
as it is for scalar equations, so we summarize it here. The main ideas in this method 
include establishing a Garding inequality for a bilinear form associated to the prin- 
cipal part of the equation in the appropriate function spaces, and then transforming 
the problem into one involving a Fredholm operator. Finally, regularity of solutions 
to the momentum constraint is discussed briefly in H3.5I 

3.1. Weak formulation. Let {A4,h) be a 3-dimensional Riemannian manifold, 
where A4 is a smooth, compact manifold with Lipschitz boundary dA4, and h G 
C^{A4,2) is a positive definite metric. Introduce the bilinear form 

ajL : W^'^ X W^-^ M, ai[^{w,w) {Cw,Cw) + (mrwtu, tr^w)^, (3.1) 

where the Robin tensor field K G L°° {dAA]N ,2) is symmetric and satisfies the 
bound 

K||trwHlw (^'trivi«,trivw)w, Viu G (3.2) 
and where K is a non-negative constant. Fix the functionals b*, b* G W~^'^ . Fix 
a function 4> G L°° and introduce the linear functional 

VeW^B'', V:=(b.0^)* + b*, (3.3) 

We used the subscript <f> in f to emphasize that (f) is not a variable of the problem. 
The functional f is a generalization of the functional IF defined in Eq. (|2.13p . 
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The weak Dirichlet-Robin boundary value formulation for the momentum 
constraint is the following: Fix Dirichlet and Robin boundary data 

ivn e W^^^OMb, 1), K e W-^'^idM,M, 1), (3.4) 

and introduce an extension Wjo of the Dirichlet boundary data as described in H2.3I 
Introduce the non- principal part operator including the Robin boundary conditions, 

e f4^) ■■= f^ppiyi) - w]^{^r,My^), (3.5) 

where f is given by Eq. (|3.3p : Let A^'^ be the afhne space given in Eq. (|2.30p for 
the case p = 2; Then, find an element to e A^'^ solution of 

a^{w,ui)+f^{ui)^Q Va;eT^li^ (3.6) 

It is convenient to express Eq. (j3.6p in terms of operators instead of bilinear 
forms. Introduce the operator 

Hence, Eq. (|3.6p written in terms of operators is the following: find an element 
w £ A ' solution of 

A,^w + f^ = 0. (3.7) 

Lemma 2. Every smooth solution w of the classical Eq. \2.11^ for a given smooth 
function (f> is also a solution of the Eqs. i3.6]) . 

Proof. (Lemma\^) The proof is similar to the proof of Lemma [H and we do not 
reproduce it here. □ 

3.2. Generalized Korn's inequalities. The Korn inequalities are a funda- 
mental step in proving existence of solutions to the linearized displacement-traction 
equations in elasticity. The inequalities involve the Killing operator i : W^'^ — > 
L^{Ai,2) with action {£u)ab '■= ^aUb + VbUa- There are two main inequalities, 
called "without" or "with boundary conditions", which can be described in terms 
of the bilinear form ae : W^'^ x W^'^ K with action ai{u,v) := {£u,£v). The 
former inequality says that the bilinear form ai satisfies Garding's inequality, 
that is, there exists fco > such that 

The latter inequality says that the bilinear form ai is coercive in the space W]^ 
in the case that meas{dMio) ^ 0, that is, there exists a constant fco > such that 

fcoiiMii?,2 ^ iinr yuew]r,\ 

These inequalities were first established in the case that the manifold Ai C M.^ 
and the metric hat is the Euclidean metric [3S1 , with new proofs given in [21] . 
A review of elasticity theory is nicely presented in [M] with Korn's inequalities 
discussed on Volume II, pages 10-13. See also [44]. Both types of Korn's inequalities 
for the Killing operator have been generalized to Riemannian manifolds in [12j . 

We just mention here that the Garding type inequality on the particular case 
of the spaces WQ''^{M,n) can be proven for a general class of bilinear forms called 
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strongly elliptic. See [58], exercise 22.7b, page 396. A bilinear form a : WQ''^{A4,n)x 
Wq''^{M, 7i) ^ R with action 

Jm 

is strongly elliptic iff there exists a positive constant uq such that 

for all vectors C G M'^ and all tensors Uci---c„ G R'^". An example of a strongly 
elliptic form is the bilinear form ai. 

The role played in elasticity theory by the Killing operator £ is played in the 
momentum constraint Eq. (|2.32p by the conformal Killing operator C, which is 
defined in Ea. (|2.7p . Inequalities similar to those satisfied by the Killing operator 
can be obtained for the conformal Killing operator, called here generalized Korn's 
inequalities. First notice that the bilinear form ac '■ W'q'^ x W^''^ R given by 
ac{u, v) — {Cu, Cv) is strongly elliptic, as the following calculation shows: 

Hence, a Garding type inequality is satisfied by the bilinear form ac on the Hilbert 
space W^J'^. However, this space is too small in our case where we need the same 
inequality on the space 'W^''^ . In addition, later we will need the coercivity type 
inequality for the bilinear form ac on the space . 

We first review the generalized Korn inequality without boundary conditions, 
which has been proven in [18 in the case where M. C R", with n > 3, and hab is 
the Euclidean metric. It is also shown in [TH] that the inequality does not hold for 
n — 1 where the null space of the conformal Killing operator is infinite dimensional. 
It is also mentioned in that article that the same arguments given in |121 imply 
that the generalized Korn inequality without boundary conditions also holds on a 
Riemannian manifold. We summarize these ideas in the following result. 

Lemma 3. (Garding's inequality for £) Let {M,hab) be a 3- dimensional, 
compact. Riemannian manifold, with Lipschitz boundary, and with a metric h £ 
C^{A4, 2). Then, there exists a positive constant ko such that the following inequal- 
ity holds 

ko\\u\\l2^\\u\\^ + \\Cuf yueW^''^. (3.8) 

Proof. (Lemmal^i) See [18] for the proof. □ 
Using Lemma [3] it is not difficult to establish that the same type of inequality is 
satisfied by the bilinear form a^. 

Corollary 1. (Garding's inequality for Ojj,) Let {A4,hab) be a 3- dimensional, 
compact, Riemannian manifold, with Lipschitz boundary and with a metric h £ 
C^{Ai,2). Let Oil be the bilinear form defined in Eq. \3.1\) for any tensor IK £ 
L°°{dAiiN,2). Then, there exists a positive constant ki such that the following 
inequality holds 

kiMl^i^Wuf + a,L{u,u) VnGW^'^. (3.9) 
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Remark. This resuh holds for both cases dAin ^ and dMn ~ 0, and also 
notice that the Robin tensor field K is arbitrary; we do not require this tensor to 
be positive definite. 

Proof. ( Corollary QJ ) The definition of the bilinear form in Eq. (|3.ip implies 

aiLiu, u) = \\Cuf + (mrivH,tr^M)w Vit G W^''^ . 

Recalling that IK G L°°{dAi]if, 2), then the second term on the right hand side can 
be bounded as follows: 



(l|hf + .^I|Vnf), 



for every non-zero number e. Let ki := ||_Zft'||oo/2, and then compute 



ajL(tt, u) ^ \\Cu\\'^ ^ ^i( ^ 11^ 



^ - u 



ko\\u\\l2-h[^\\u\\ 



11. 2 



Choose the number e such that ko — kie^ — fco/2, then 



ajL(M, tt) ^ -(l 



3^1 N 



u 



u 



1,2; 



and so, 



2 



11.2 



«c 1 



^ 1 



2fc^' 
2fcJ^' 



Divide by (l + and set ki 



if + at(M, n) . 
and the result is the ineauahtv (|3.9p . □ 



(2feo+3fef) 

We have not found in the literature the generalized Korn inequality with bound- 
ary conditions on only part of the manifold boundary, neither in Euclidean space 
nor in an arbitrary Riemannian manifold. This type of inequality is crucial in ^3.31 
so we proceed to establish this result. 

Lemma 4. (Coercivity of C) Let {A4,hab) be a 3- dimensional, compact, Rie- 
mannian manifold, with Lipschitz boundary such that meas{dA4io) > 0, and the 
metric h G C^(A1,2). Then, there exists a positive constant fco such that the fol- 
lowing inequality holds 

fco||M||?,2^ ll^wf yueW]^\ (3.10) 

Proof. (Lemma\^) The proof has two main parts: The first one is to show that the 
null space of the operator C : W]^ L'^{M, 2) is trivial when meas^dMio) > 0; 
the second part uses the Carding type inequality satisfied by the operator C and 
presented in Lemma [3] together with a well known argument by contradiction to 
show Eq. (|3l0ll . 



22 



M. HOLST, J. KOMMEMI, AND G. NAGY 



The first part mentioned above also consists of two steps. We first step is to 
sliow that any vector field belonging to the null space of £, vectors called confor- 
nial Killing vectors, must satisfy a particular set of ordinary differential equations 
(ODE). Indeed, assume that w° is a conformal Killing vector, so Cu = 0, and 
introduce the fields 

where we introduced the notation V[QUfc] :— (VaUfc — VbUa)/'^, and similarly we 
will denote V(QUfc) := (VaUh + Vfcita)/2. A straightforward albeit long computation 
commuting derivatives shows that a conformal Killing vector and its derivatives 
introduced above must satisfy the following equations, 

VaUb = aab + ^/3hab, (3-11) 

Va/3 = 7a, (3.12) 
2 

"^aOibc = -Rbc/ud + - l[bK]a, (3.13) 

Va76 = -3u^V,La6 - 2l3Lab ~ 6Rc(aabf, (3.14) 

where the tensor Rabc^ is the Riemann tensor of the metric connection Va, the 
tensor Rab = Racb'^ is the Ricci tensor, and we have introduced the tensor Lab '■— 
Rab — Rhab/4:. The first two equations above are the definitions of the fields aab, 
(3 and 7a. The other two equations are obtained by commuting second and third 
derivatives of the conformal Killing vector u'^. They are generalizations of the well- 
known formulas for Killing vectors (where /3 = 0, 7a = 0) which can be found for 
example in [5D], page 443. These formulas in the case of Lorentzian metrics have 
been used in [53]. Contract Eqs. p.lip - p.l4p on index a with any vector field v"", 
and the result is a system of ODE for the fields aab, (3 and 7a. From this system 
of ODE we conclude the following: If these four fields vanish at a single point in 
A4, then they vanish identically on A4. 

The second step is to show the following: If is a conformal Killing vector that 
vanishes on a two-dimensional hypersurface dMjo with meas{dMiD) > 0, then the 
vector m", and the fields aab, P and 7a vanish at any point on the hypersurface 
dA4io- This statement and the conclusion of the paragraph above will imply that 
vanishes identically on the manifold 7W. Denote by n° the unit vector field 
normal to the tangent space at each point in the manifold dA4]o, and introduce 
the first and second fundamental forms of the hypersurface dMjo as follows. 

Lb ■= hab - riaUb, Kab ■= -la"^ cnb, 

where the tensor Kab is symmetric. Denote by Da the Levi-Civita connection as- 
sociated with the two- metric lab defined on the 2-dimensional hypersurface dAdjo, 
so the connection satisfies the property that Dahc = 0. Extend the vector field 
n° to a neighborhood of the hypersurface dA4]o in the manifold A4 as the tan- 
gent vector solution to the geodesic equation nf^Van^ = with initial data n° on 
dA4.]D- Hence, the resulting vector field satisfies n°'Van^ = in a neighborhood of 
the hypersurface dAijo- Then, decompose the conformal Killing field u'^ as follows 
M° = w„n° -|- w", with TT-au" = 0. Denote by jjo evaluation at the hypersurface 
dM]D, then the condition u"'\]o =0 implies 

Un\D^Q, M"liD=0. 
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The latter condition means that {DaUn)\iD = and {DaU^)\nj = 0, while the latter 
equation together with the equation n"-n''{£u)ab = 0, which also holds on dAin, 
imply that (V„u„)|iD = 0, where we use the notation V„ := n"Va. Therefore, from 
expression VqW" = V„u„ + Dau"" + Ka°'Un we then conclude that 

I3\n = 0. 

The equation n°'lJ'{Cu)ab — implies, after a short calculation, that the equation 
IcWnUb + KcbU^ + DcUn = holds in the manifold M, and so it also holds on the 
hypersurface dA4]o- This result together with our previous results establish the 
condition (Zc''V„U6)|jd = 0. The decomposition 

and our previous results then imply that 

aab\lD = 0. 

We still have to show that the vector field 7a ~ Vq/3 vanishes on the hypersur- 
face dAijD- Since /3 is a scalar field and vanishes on dAinj, we conclude that 
the field {Da(3)\]D = 0, which impUes {Da'^nUn)\iD = 0. We now only need 
to compute the field (V„/?)|iD- The identity n^lJ'\i'[a'^b]Un = when evaluated 
on the hypersurface dMjo together with our previous results imply the equation 
(£'cV„u„)|iD = (Zc^V„£'fcM„)|jD. But we just showed that the left hand side van- 
ishes, and so then does the right hand side (?c''V„£'bU„)|iD — 0. Finally, from the 
equation n°/rf''(2V[QVf,]Uc — Rabc'^Ue) = evaluated on the hypersurface dA4io we 
conclude that 

{lc''ld''VnDaUb)\n - [Dc{ld''VnUb)]\n = 0. 

The result {lJ'^nUb)\]D = implies that the second term on the left hand side 
above vanishes, so we conclude that {lc°'ld''^nDaUb)\iD = 0, and from this equation 
one gets [V„ (-Daw")] |jd = 0. Therefore, in order to show that the field (V„/9)|jd 
vanishes we only have left to prove that the field (V„V„u„)|jd vanishes. That 
this is the case follows from the equation Vn[n'^n'^{^u)ab] — 0, which holds in the 
manifold M, and so on the hypersurface dAdjo, and an explicit computation shows 
that (V„V„u„)|jD = 0. We then conclude that 

lain = 0. 

Let us now recall that the ODE equations obtained by contracting Eqs. (|3.1ip - 
p.l4p on index a with any vector field v"" are homogeneous on the fields aab, 
(3 and "fa with vanishing initial data on the hypersurface dAijo- The solution 
vanish identically in a neighborhood of this hypersurface dAijo in the manifold 
M.. Repeating this procedure we conclude that the conformal Killing vector field 
vanishes identically in M . This result establishes that the null space of the operator 
C is trivial on the space VF^^. 

We now consider the second part of the proof of Lemma |4] using a well-known 
argument by contradiction. Assume that there exists a sequence {it„} C 
such that 

||itn||i,2 = 1, and lim = 0. 

' n — >oo 

The sequence {tt„} is bounded in W]^ which is a reflexive Banach space, so there 
exists a subsequence, also denoted as {un}, such that 

^ Uq in W]i^, and u„ ^ tto in X^, 



24 



M. HOLST, J. KOMMEMI, AND G. NAGY 



the latter statement following from the imbedding W]i^ C being compact. So 
{un} is a Cauchy sequence in L^, and by assumption the sequence {Cun} C is 
also a Cauchy sequence. The Carding inequality in Lemma [3] implies that 

fco ||lt„ - Itmll^ 2 - ItrnlP + ||£u„ - >CWm|P ^ aS n,m—>- OO, 

and so the sequence {un} is also a Cauchy sequence in W]^. We then conclude 
that 

Un ^ Uq in W]^ CUn ^ = Cuo. 

But the null space of the operator C is trivial on the space W]^ , therefore we 
conclude that the element ito = 0. However, this leads us to a contradiction from 
the hypothesis that ||itn||i,2 = 1 which implies that ||i«o||i,2 — 1 so the element 
tto ^ 0. Therefore, such sequence {iXn} does not exist, which then establishes the 
Lemma. □ 

Corollary 2. (Coercivity of a^) Let {A4,hab) be a 3- dimensional, compact, 
Riemannian manifold, with Lipschitz boundary such that meas{dAiio) > 0, and the 
metric h G C^{M,2). Let ajj, be the bilinear form defined in Eq. and assume 

that the Robin tensor IK £ L°^'{dA4pf,2) is positive definite. Then, there exists a 
positive constant ki such that the following inequality holds 

^1 l|w||?,2 ai(u,it) Vite Wl^,^ (3.15) 

Proof. ( Corollary O ) Since the Robin tensor M is positive definite, then the 
result is straightforward from Eq. p.lOp . due to the following inequalities, 

\\Cuf+k\\Ui, 

s$ \\Cuf + (mriv«,trww)iv 

^ aj[^{u, u), Vn e W ]^ . 

This inequality establishes the Corollary. □ 
Remark. It can be shown that the result in the Corollary [2] remains valid if the 
Robin tensor field K is slightly negative definite. 

3.3. Results using variational methods. The momentum constraint Eq. (|3.6[) 
can be written as the Euler condition for stationarity of a real-valued functional on 
a Banach space. Direct methods in the calculus of variations can be used to find 
the points that minimize this functional in the Banach space in the case that the 
hypersurface dAin ^ 0. The main concepts needed from the calculus of variations 
are summarized in the Appendix, where we also explain the part of the notation used 
in this Section. Since the momentum constraint is linear, we will achieve similar (in 
fact, slightly more general) results in §3.41 using Riesz-Schauder Theory. However, 
the presentation here is a guide for our variational treatment of the Hamiltonian 
constraint equation in > j4.51 and the results we assemble in this section are of interest 
in their own right. 

Let ajL : x W]i^ ^ R be a bilinear form with action defined in Eq. (|3.ip . 

and fix the functionals ft*, ft* G ' . Let Wjo & W^'^ be the extension of the 

Dirichlet data Wjo, and ib^ be the Robin data functional, both data defined in 
Eq. (|3.4p . Introduce the functional 

Jj,:W]i^^R, J^(u;) := ia^(u;,u;) + G^(u;), (3.16) 
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where the functional is given by 

G^{u}) :=.90(iOjD - io^(tr^u;) + ai(wjD,'^), 
with the functional g^{w) having the form 

g^w) ■.= {br^y{w) + b*{w), 
and we will use the notation w = Wjo + i^. 

1,2 



(3.17) 



be the func- 



Theorem 1. (Existence of a minimizer) Let Jji^ : W 
tional defined in Eq. \3.16]) . Assume that the hypersurface dAdjo =/= 0, and fix 
an extension of the Dirichlet boundary data Wjo G W ' and the Robin boundary 
data w*jj^ e W^2''^(^dA4ir~f,l). Fix the functionals 6*, b* G W jj^'"^ , and the tensor 
IK G L°°{dA4M,2) satisfying the inequality in Eq. \3. ^|) with K ^ 0. Then, there 
exists a unique element uj G W minimizer of the functional Jjj, on W , that 
is, 

Je,(u3) = inf Jjl{ui). 



Proof. (Theorem[l[) We start showing that the functional J^j, is coercive, and the 
first step is the following inequality 

JiL{i^) ^ ^ a]L{uj,uj) - \G^{uj)\. 

The linear terms in G^ can be bounded as follows: recall the notation w = w^j +a;, 
then 



\g^{w)\^\{br4>''r{w)\ + \b*{w)\ 



-1,2 



-1.2 



l"'l|l,2, 



introducing the constant Cg :— 
\g^{w)\ < 2cg ||i«||i^2 ^-c 



lo\\K\ 



-1.2 + ||&j||-i.2]/2, we then obtain 



2 ' e||-u;|i2 



1.2' ^ g S 



2e\\w^\\l^ + 2e\\u:\\l^, (3.18) 



where e is any positive constant. The second term in the functional can be 
bounded as follows 

l-w'wlti'ww)! ^ ll«'wll-i,2.w l|trw'^||i,2,w 

Co ||W^||_1,2,W ll<^lll,2 



2e 



-no;. 



ii,2> (3.19) 
where Cq is a positive constant. The third term in G^ can be bounded as follows 

\ai,{Wn-,'J^)\ < C2L |1i«2d||i,2 ||t^||l,2 

|2 , - II, .112 ^3 

By adding the inequalities in Eqs. (|3.18p - (|3.20p we obtain the bound on G^, 

\G^{i^)\^CG + ie\Ml^, (3.21) 

„2 



<f ll«'«lll,2 + f ll'^l 



'-SvlLi^2,W 
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We now consider the bilinear form at. First, the tensor K satisfies the inequaUty 
in Eq. ()3.2p with K ^ 0, so we have 

aj[^{<jj,<jj) > {Cuj,£u)). 

The assumption that the hypersurface dAijo ^ and the generahzed Korn in- 
equahty in Lemma |4] imply that there exists a positive constant kg such that 
(jCw, £-u>) ^ ko \\<^\\i 2i which together with the inequality above imply 

air.{i^,uj)^koMl2. (3.22) 

Therefore, from the inequalities ((3?2T|) - ((3?22)) we obtain 

Ji(a;) ^ (fco - 3e)||u;||2 2 -cg. 

Choosing e small enough we have established that Jjj, is coercive in W^^. 

We now show that the functional Jjl is Isc^, . Let {a;„} C W],^ be a sequence 
such that iiJn ~^ in id j which then implies that a;„ —> ljq in L^. We again 
start with the functional G^, which is linear on its variable a; G VF^q , therefore it 
is continuous under weak convergence (by definition of weak convergence). So it is 
also IsCtu , and the following equation holds 

Gd,ii^a) = liminf GJiVn). 

n — >oc 

We only have to show that the functional <jJ ^ a][^{uj,uj) given in Eq. (|3.ip is 
also Isc^ . The first term in the bilinear form a^j, defines a norm in W^j^ , since 
the generalized Korn inequality given in Lemma 2] and the fact that the conformal 
Killing operator is bounded in W]^ imply that there exist positive constants fcoi 
Ko such that 

fco MI2 < ll^'^f ^ Ko \Ml2, Va; e W]^^ 

This last inequality means that the map u i— > H^CwH defines a norm in W]i^, and 
so it is IsCiu , a result proven in the Appendix. Therefore, the following inequality 
holds, 

|l£a;of liminf \\CuJnf. 

n — >oo 

The second term in the definition of the bilinear form aj^ contains the two-index 
tensor IK, which is positive definite and symmetric, therefore the function 

is a continuous and convex functional, and so it is IsCm , a result also proven in the 
Appendix. We then conclude that the functional Jjl is lsc„ . Therefore, Theorem[T51 
in the Appendix in the case U = X shows that there exists a minimizer for Jjj, in 

The uniqueness of the minimizer is a consequence of the strict convexity of the 
functional Ji, which is a general result that, once again, is established in the 
Appendix. We have to show that for all non-zero e W]i^ and all t G (0, 1) 

holds 

JiL{t(:j + {1- t)ui) < tJ^{Cj) + t)J^{yj_), 
or equivalently, as it is explained in the Appendix, we only have to show that for 
all non-zero u), cj G W]if holds 

DJil{uj,){u}) < J]l{u} + uy) - Jil{i^). 
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A straightforward calculation shows that for all cD, a; G holds 
Ji(a> + - JjlM = ^ a2L('i,a>) + ai(u;,a>) + G^{GS) 

where the symmetry of the Robin two-tensor field IK is used to establish the first 
line, and the last line is obtained from the generalized Korn's inequality Eq. (|3.10|) . 
Therefore, the functional Jn^ is strictly convex, hence. Theorem 1 141 in the Appendix 
implies that the minimizer u) is unique. This establishes the Theorem. □ 
The next result shows that the minimum uj of the functional on the space 
found in Theorem [T] is a solution of the Eulcr equation DJn^itj) = 0. 

Theorem 2. (Momentum constraint) Assume the hypotheses in Theorem [H 
Then, the functional J^^ is Gateaux dijferentiable on W]^ and the minimizer uj £ 
W JO is solution of the Euler equation 

where the equation above is the momentum constraint Eq. US. 6\} . 

Proof. ( Theorem l^) It is straightforward to verify that the functional Jjl is 
Gateaux differentiable, and its derivative at an arbitrary element u) S W j,^ is 
given by 

with w := W]o +u). Therefore, the Gateaux derivative -DJjl is the left hand side in 
Eq. (|3.6p . Let uj G W],^ be the minimizer of the functional on the space W]i^. 
Then the following inequality holds, 

DJ^{uj){ui) 0, Vwe t^^^ (3.23) 

For the proof, write down the Gateaux derivative of the functional Jjl at the min- 
imizer UJ, 

DJe.{u;){uj_) = ^lirri^ [Jjl('^ + tLo) - Je.{(^)] ■ 

The element w is a minimizer of Jjl, so Je,{(j^ + tuj) ^ Jjj^{u), which establishes 
Eq. (IS231). This Eq. (jX^ holds for -w, so we conclude that DJjl(w) = 0. This 
establishes the Theorem. □ 

3.4. Results using Riesz-Schauder theory. We present here a proof of ex- 
istence and uniqueness of solutions of the weak Dirichlet-Robin boundary value 
problem for the momentum constraint Eq. p.7p . The proof is based on the Riesz- 
Schauder theory for compact operators, see |51| . The proof is more general than 
the one given in i )3. 31 because it includes the case where meas(9A^jD) = 0, that is, 
the pure Robin case. Riesz-Schauder theory was used for the momentum constraint 
in [28j to develop an approximation theory and corresponding error estimates for 
numerical approximations. 
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Theorem 3. (Momentum constraint) Consider the weak formulation for the 
momentum constraint Eq. |j'.7| ). Assume that the Robin tensor field IK satisfies 
Eq. US. 2\) with positive constant K. Then, there exists a unique solution w G A^'^ 
to the momentum constraint Eq. Iff. 7| ), and there exist positive constants Ci and C2 
such that the following estimate holds, 

IIH|l,2 UWL \\K\\-1,2 + ||b*||-l,2 +Ci ||tb:^|Ui,2,w + C2 \\Wn\\l.2. (3.24) 

Proof. (Theorem\^) First translate the problem from the affine space A^'^ into a 
problem on the vector space Wj,^ with the change of variable w = Wjo +u). Then, 
Eq. (|3.7p has the form: Find a; G solution of 

Aj^uj + G0 = (3.25) 

where G^{uj_) := f^{uj) + ai[^{W]D,uj), and we now consider the operator Ai[^ : 
W^i^'"^ . To find a solution of Eq. (|3.25p is equivalent to show that this 
operator is invertible. Lemma [3] says that the conformal Killing operator C 
satisfies Garding's inequality Eq. (|3.8|) . This implies that the bilinear form ajj, also 
satisfies a Garding inequality, which was proven in Corollary [TJ Then, Theorem 1 121 
in the Appendix implies that the operator An^ is Fredholm with index zero. That 
means dim Naj^^ — codim Rajj^ , which can be described saying that the operator 
is bijective iff it is injective. This property is described in the PDF literature 
as "uniqueness implies existence" . So, in order to show that Au^ is invertible we 
only have to show that its null space is trivial. Consider an element u £ 'W^^ such 
that A^u = 0. In particular A^u[u) — 0, which is equivalent to 

= \\Cuf + {KXx^u, trw«)iv > \\Cu\^ + K ||tr^u||^. 

Both terms must vanish, since the tensor is strictly positive definite and so the 
constant K is positive. From the first term one obtains that m is a conformal Killing 
vector, and from the second term together with u £ W],^ one obtains that trit = 
on the whole boundary dAi and so, u £ Wq . Therefore, tt = in the manifold 
A4 since the bilinear form Ojl is strongly elliptic. So, the null space of the operator 
is trivial, and then is invertible. Finally, it is not difficult to check that the 
estimate given in Eq. (|A.6|1 of the Appendix applied to u; = to — Wjo implies the 
estimate on w given in Eq. (|3.24p . This establishes the Theorem. □ 
Remark. In the case that the hypcrsurface dMio ^ the Robin tensor IK need not 
to be strictly positive. There exists a unique solution to the momentum constraint 
in the case that the constant K in Eq. (|3.2[) is slightly negative, that is, K > — Kq 
for small enough Ko > 0. The proof uses the coercivity of the conformal Killing 
operator £, the inequality p.lOp in Lemma IH instead of the Garding inequality. It 
can be shown that the assumption that K is negative but not too negative implies 
that the bilinear form Oj^ itself is strictly positive on W]^. Recalling that the linear 
form : W j,^ — > M is bounded, then the Riesz representation Theorem says that 
there exists a unique u £ solution to the weak problem with Eq. (13. 6p . In 

terms of the operator A-n^, this statement means that A^ ■ W]o W 

ID ' IS 

invertible. 

3.5. Regularity of solutions. In this Section we state without proof regularity 
results, which can be obtained from the literature, and are applied to the weak 
solutions of the momentum constraint. 
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Theorem 4. (Regularity W ^ A ssume the hypotheses in Theorem and in 
addition assume that the boundary set dA4 is C^'^. Assume that the source func- 
tional f^p and the boundary data satisfy the following conditions, 

K, b* G W-'^P, wjo e W^^P, wl, E W-i'PidMjM, 1), P>2, 

then, the solution w to the momentum constraint Eq. Jg. 7| ) satisfies that w E W^'^ 
and there exist positive constants Ci and C2 such that the following estimate holds, 

IIHIi^p ^ ml \\K\\-i,P + \\b*\\-i,p + ci ll^i^wlUi^p.w + C2 Ikclli.p. (3.26) 

Proof. (Theorem^) We only describe a sketch of the proof. See for example pS] , 
Sec also [TT] for interior estimates only, Theorems in §7 and §8. These results can 
be extended up to the boundary for smooth enough boundaries. □ 
We also present here a result from [T3], stating higher regularity of the weak 
solution of the momentum constraint Eq. Ij3.7p in the case that the data and the 
source function also possess additional regularity. 

Theorem 5. (Regularity W^'^) Assume the hypotheses in Theorem\^ and in 
addition assume that the boundary set dAi is C^. Assume that the source func- 
tional have the form = {bT^uj), and = ibj^uj), while the boundary 

data have the form Wj^{^f in^ ~ {wimt^^in^in , for all u G W]^ . If the following 
conditions hold 

o 

then, the solution w to the momentum constraint Eq. Jg. 7\ ) satisfies that w G W^^'^ 
and there exist positive constants c\ and C2 such that the following estimate holds, 

M2.P < U\L \\br\\p + \\bj\\p + Ci ||-U;w|U,p,w + C2 \\Wb\\2,p. (3.27) 

Proof. (Theorem[5i) We only describe a sketch of the proof, which follows [Hj . 
Vol. II, page 296. It is based on the fact that the momentum constraint bilinear 
form OjL is strongly elliptic and satisfies the supplementary and complementing 
conditions given in [2]. □ 

4. The Hamiltonian constraint 

In this section we fix a particular functional in an appropriate space and 
we then look for weak solutions only of the Hamiltonian constraint Eq. (|2.35p . 
We first develop the weak formulation more precisely in §4.1) and as in §2.31 we 
assume the weakest regularity of the equation coefficients such that the equation 
itself is well-defined. As was the case for the momentum constraint, we will be able 
to use variational methods to obtain existence (and when possible, uniqueness) 
results for the Hamiltonian constraint in this weakest setting. First, we establish 
some preliminary results on generalized local and global barriers (constant sub- and 
super-solutions) for weak solutions in HA.2[ The term local means that the barrier 
does not depend on the coefficient a^, while global means the barrier does depend 
on this coefficient. We summarize the generalized local and global barriers in ij4.3l 
In ^AAl we establish some related a priori L°°-bounds on any VF^^^-solution to the 
Hamiltonian constraint. In §4.5) we then use the barriers from §4.2) together with 
a variational argument, to establish existence, and when possible uniqueness, of so- 
lutions to the Hamiltonian constraint in the weakest possible setting of L°° DW^''^. 
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Due to the lack of Gateaux-differentiability of the nonlinearity in W^'^, the con- 
nection between the energy used for the variational argument and the Hamiltonian 
constraint as its Euler condition is non-trivial, and is established through several 
Lemmas. In i )4.6l we give a second (non-variational) argument for existence, using a 
barriers approach as in most of the earlier work [29[ I30j , which requires additional 
regularity on the equation coefficients. Regularity of solutions is discussed briefly 
in 

The results obtained using variational methods in t J4.5l can be viewed as lowering 
the regularity of the recent result of Maxwell on "rough" CMC solutions in W'^''^ 
for k > 3/2 down to L°° n W^'^. We note that the barrier-based existence results 
for the Hamiltonian constraint equation in §4.61 and the compactness argument in 
§5. II giving existence for the coupled non-CMC system, require higher regularity on 
the equation coefficients. However, we still end up with some non-CMC results for 
the coupled system in weaker settings and in more general physical situations than 
have been previously obtained. These additional assumptions are clearly stated in 
those sections. 

4.1. Weak formulation. Let {A4,h) be a 3-dimensional Riemannian manifold, 
where is a smooth, compact manifold with Lipschitz boundary dA4, and h S 
C^(A^,2) is a positive definite metric. Introduce the bilinear form 

ai^-.W^-^ xW^-^ ^R, ai((/),0) (V0,V0) + (iftr^,/),trjv0)^, (4.1) 

where the Robin function K G L°°(97Wjv,0) satisfies the bound 

k lltr^^ll^ < (A'tr^0, tr^0)jv, V0 G W^'^, (4.2) 

with k being a non- negative constant. Fix the functionals 

a*eW^^^\ a;eW-l^\ aleW^^^\ (4.3) 

The assumption on the background metric implies that the function an is continuous 
on the manifold A4, so the functional e W^^'^ given by 

a*n{^):^{an,^), V^eW'^'^ (4.4) 

is well-defined. Given any two functions (pi, 02 G L°° with Q < (pi ^ (p2, define the 
interval 

[4>i,4>2] := {</i G : (pi ^ (p ^ (p2}, 

which is a closed, bounded set in L°°, and also in L^. Introduce the nonlinear 
operator 

f^p{(P) := [a^cf")* + [aRCpT - (flp^')* " M^'T- (4-5) 

We used the subscript w in /^^ to emphasize that w is not a variable for the 
analysis of the Hamiltonian constraint in isolation from the momentum constraint. 
The functional fwp is the generalization of the functional F defined in Eq. (|2.I2p . 
We remark that the operator defined in Eq. (j4.5p is continuous but not Gateaux 
dijferentiable. It has Gateaux derivatives only along directions in L°°, not in the 
whole space L^. This property of the functional fwF will introduce some technical 
complexity in the use of variational methods for the Hamiltonian constraint (see 
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The weak Dirichlet-Robin boundary value formulation for the Hamilton- 
ian constraint is the following: Fix Dirichlet and Robin boundary data 

< ess inf $d $d e L°°{dMD,0) nW^^-^{dMD,0), $*n £ W~^'^{dMN,0); 

Mo 

(4.6) 

Introduce an extension (pjj of the Dirichlet boundary data as explained in i i2.31 in 
particular, given ipo > on dAioi we can use the Laplace-Beltrami operator to 
harmonically extend 4)o to 4)o such that ^j, > a.e. in M. Given such extension 
function (pu: fix any two functions 4>i, (j)2 G L°° n W'^''^, with the property that 
(/)2 and such that G [^i,?!'2] H W^'^; Introduce the non-principal part 

operator including the Robin boundary conditions, 

U ■■ C L2 ^ W-''^, Um^) Urmf) ~ 0^(triv^), (4.7) 

where the functional fwp is given by Eq. ((43| : Let A^'^ be the affine space defined in 
Eq. (|2.29p . which includes the Dirichlet boundary condition; Then, find an element 
6 [(f>i, 412] n A^''^ solution of the equation 

ai(</),^) + /^(0)M = yipeW^'\ (4.8) 

As was the case earlier for analysis of the momentum constraint, it is convenient 
to express Eq. (|4.8p in terms of operators instead of bilinear forms. Introduce the 
operator 

Also recall that, if given any ip £ [(/)i,02], then fwiip) G W^^'^. Hence, Eq. (|4.8[) 
written in terms of operators is the following: find an element (j> G [0i,'?!'2] H A^-^ 
solution of 

Ar^cl) + U{<j>) = 0. (4.9) 

Lemma 5. Given a smooth vector field w, every smooth function (f> solution of the 
classical Dirichlet- Rohin boundary value formulation for the Hamiltonian constraint 
Eq. h2. 1 6\) is also a solution of the weak formulation with Eq. |^.ig[ ) corresponding to 
the equation coefficients and Rohin data function given by the following expressions, 
which hold for all ip e W^'^ , 

a*^{(£) -.^ {ar,'£), a*{i£) -.^ {ap,i£), 

o-loi^) iaw,v), 'Plfi^^N^p) (0jv, trjv(^)Ar. 

Proof. (Lemma\Bi) The proof is similar to the proof of Lemma [T] and it is not 
reproduced here. □ 
Given any function u G W^'"^, recall the notation 

:~ ess max{u, 0}, := —ess min{it, 0}. 

An element </>_ G W^'^ is called a sub-solution of Eq. (|4.9p iff the function 
satisfies the inequalities 

(0B - 0_)- G W^'^ and - + fM-)] G W^l^\ (4.10) 

An element 0+ G W^'^ is called a super-solution of Eq. (14. 9p iff the scalar function 
(/)-)- satisfies the inequalities 

{(t>n - 0+)+ G W^'^ and [A^</>+ + /»(</>+)] G W^l^^ (4.11) 
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The sub and super-solutions of Eq. (|4.9|) may depend on the choice of . A sub- 
solution is called global iff Eq. (|4.10p holds for every functional G ' , and 
it is called local iff it is not global. 



4.2. Global and local barriers. In this section we show that there exist sub- and 
super-solutions to the Hamiltonian constraint equation (j4.9|) for different assump- 
tions on the equation coefficients. The results in Lemmas [Bl IIOl are generalizations 
to the weak problem of the barriers found in 30J in the case of closed manifolds, 
scalar curvature R — —I, and equation coefficients with higher regularity. The main 
idea of this generalization is to look at candidates for sub- and super-solutions only 
among the constant functions, and not among all functions in [(f>i, (f)2] C where 
< (pi ^ (f>2- This type of approach is reasonable, since in the smooth coefficient 
case there exist sub- and super-solutions which are indeed constants. 

Assume that the background metric h belongs to C^{A4, 2), then the Ricci scalar 
of curvature i? is a continuous function on the manifold A4. Introduce the constants 

:= sup II - , := sup - , (4.12) 

A "p(^) A «^('^) -,o^ 

a := sup II ~ , sup - . (4.13) 

In order that the Lemmas below also hold for the particular case when the Ricci 
scalar R vanishes identically, we introduce the constant :— l + a'j^. Given a two- 
index tensor a € LP{A4, 2) and a vector field w £ W with p = 12/5, introduce 
the functionals a* and a£^„ given by a* ((^) — (cr^, </j)/8 and a'}^^{^p) — ({jCw)'^,(f) /8, 
where (j) € M^d'^- Now introduce the further constants 

a<T — sup II - , K sup II I ~ , 

7a rA^rM^) 

0JV := sup — — , 0o:=sup0o, 



+ 



ll¥'l|l,2 M 



where we recall that the function c/jjj is the harmonic extension of the Dirichlet 
boundary data discussed in tJ2.3l In an analogous way, switching sup to inf, 
introduce the quantities a^, a^, a^, a^, a^, X^, (f)'^ and (f)"^. 

Lemma 6. (Local super-solution R bounded) Consider the weak formulation 
for the Hamiltonian constraint given in ^4■^ Assume that the constants a^, K'^ 
are positive, and denote by 4>w+ the constant 



max I 1 , 



K}- (4.14) 



Then, is a local super- solution of Eq. j^.gp . 



Proof. (Lemma[B[) We look for a super-solution among the constant functions. 
Therefore, let be any constant in [^i,(/)2], with > 0, then the following 
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inequalities hold for every element ip £ W}'"^ 



fwFiMif) = {ar4>l)*{'£) + (aii</'o)*(^) - (ap0o^)*(^) " (a^'^o^)*!^) 

Introduce the polynomial on (jjQ given by 

qiM < 00 - 00 - < 0o ' - 0o"'. (4.15) 

The assumptions that the constants and are strictly positive, while the 
constants Up and are non-negative imply that there exists a unique positive 
root of this polynomial. The proof consists of three steps. First, there exists at 
least one positive root of the polynomial q, because for 0o large enough g(0o) is 
positive, and for 0o close to zero from positive values q{(j)o) is negative, as it can be 
seen from the following expression. 



where the term between brackets becomes negative for small enough 0o. Second, 
this positive root is unique, since the function q is increasing for all 0o > Q^o '■— 
[a^/(5a^)]^/^ (the proof is to verify that > for 0o > Q^o); and the function q 
satisfies the inequality q{(j)Q) ^ f(0o) '■— a^ 0o ~ '^r4'o for all positive numbers 0o. 
Since r{ai) = for ai := [a'^/a'^]^^'^, and ai > ao (so the root of the polynomial 
q must belong to the interval where q is increasing), we then conclude that the 
root of the polynomial q is unique. Denote by 0o the unique positive root of the 
polynomial q. Since 0o > ctQ, then q{4>o) > q{(t>Q) = for any 0o > 0o. The idea 
now is to find an upper bound for the root 0o. The result is going to be the first 
two expressions on the right hand side in Eq. (|4.14p : the remaining two terms on 
the right hand side of Eq. (|4.14p will account for the boundary contributions. 

In the case that 0o ^ 1 (which could be verified, for example, by explicit evalu- 
ation), then choose a candidate for super-solution to be 0u,-|- = 1. In the case that 
00 > then there exists an upper bound for this root, as can be seen from the 
following argument. Given any 0o ^ 1, then the following inequalities hold, 

(0o)"+i^i ^ -(0o)-("+i) ^ -1 ^ -(0o)-" -00, 

therefore, this inequality for n = 3 and n — 7 implies that for all 0o 1 holds 



9(0o) > s(0o) := (0o) 



5 



a 



p 



This new polynomial s vanishes at 



and the inequality 0i > ao says that 0i belongs to the interval where the polynomial 
q is increasing. Therefore, we have that 

g(0i) ^ s(0i) = and q(0o) = 0i > 0o. 

So, in this case 0o > 1, choose the candidate for super-solution to be 0u,+ = 0i. 
Then, introducing the constant 

1-1 " ■ ... / ^ 

:= max-^ 1, 
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we have estabhshed that the following inequality holds 

/»f(0o)M^O yipeW^i, V0o^0»+. (4.16) 

We now account for the boundary contributions. The definitions of cj)'^ and K'^ 
imply that for any constant 0o holds 

{[Kcf^o^M^^N^)^ > (i^>o-<^^i)IMIi,2 V^e w^^'t 

In particular, defining the constant 

1/4 



I ' V at \ ' /f^J' 



follows that 5^ (^w+ and the following inequality holds 

([i^</)o-0iv],trjv<^)^ V</)o^^^+. (4.17) 

Adding Eqs. (|4.16p and (|4.17|) we conclude that 

(if</.o,trjv^)jv + /^(</'o)(^) >0 V(^eT4^^-^, V0o^0.,+ . 
Recalling now that any constant 00 satisfies (V0Oj Vt^) — 0, we conclude that 

^L0o(^) + /^(0o)(^)^O V^gM^^^_^, V0o>0.,+ . 
Finally, introduce the constant as given by Eq. (|4.14p . In particular, this 
constant satisfies (00 — = and (t>io^ ^ 4'w+^ the following inequality 

holds, 

[4'D ~ 4>w+ )+ G W^c^ and [Ai0„+ + e 
which establishes that is a super-solution of Eq. (|4.9p . □ 
We now find a global super-solution for the Hamiltonian and momentum con- 
straint Eq. (j2.35|) - (j2.36|l . where global means that the super-solution is independent 
of the vector field w solution of the momentum constraint Eq. ()2.36|) . This global 
super-solution is a generalization suitable to our weak setting of the super-solution 
derived in ^30j. We use the same idea as in i j4.21 that is, we look for super-solutions 
only among the constant functions. The "near-CMC" assumption on the trace of 
the extrinsic curvature made in [30j to construct the super-solution is still present 
here, although in weaker norms. 

Lemma 7. (Global super-solution R bounded) Consider the weak formula- 
tion for the Hamiltonian and momentum constraints given in '>\2.S\ . Assume that 
the numbers (o^ — Ki) and K'^ are positive, where the constant Ki is defined in 
Eq. \4.2i\l . Denote by 4>+ the constant 



:— max < 1, 



^n}, (4.18) 



< -Ki 

with the constant K2 defined in Eq. \4.22^ . Then, the constant 0+ is a global super- 
solution of Eqs. (M^-^MB- 

Proof. (Lemma^) Consider the weak formulation in H2.3I Let (pQ be any constant 
in [01,02] with 01 > 0, then the following inequalities hold, 

/^(0o, w){^) = (a;05)((^) + (a* 0o)(<£) - (a;0o =^)(^) - «0o ')(<£) 

= a*{^P) <t>l + aUv) 00 - a*{ip) 0o ^ - al{(p) (jy^'' 

^ (a^ 0g - 4 00 - a; 00 3 - 0o ||^||i,2. (4.19) 
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The number in the last term is bounded when w is sohition of the momentum 
constraint Eq. (|2.32p with any source (p G [(/)i,(/<2]- For the proof, start with the 
definition of in §4.21 where a'^{(p) = (a^,(^) with G i^/^ and ip G W^'^. 
Then, the foUowing inequahties hold 

{a^,(p) s$ ||au,||| ||<^||6 sS Cs ||a^||| ||^||i,2, 

where we used the imbedding W^''^ C L^, and where Cg is the positive imbedding 
constant that relates the norm of these spaces. The inequality above holds for all 
if G W^'^, and in particular holds for the supremum in that space, hence 

^ Cs lla^ll 6. 

5 

Now, the definition of the function a^,, standard inequalities and the notation p = 
12/6 show that, 

||a„|| e = i Ik + Cw\\l ^ i (||a||^ + \\£w\\l) ^ i [\\a\\l + cc • 
In ^'SA\ and H'S.5\ it is shown that there exist positive constants ci and C2 such that 

^ ii0iiLiib:ii-i,p + iife;ii-i,p+ciiiii^wiUi,p,w+c2ikBiii,p. 

Then, the bound for the number can be written as 

^Ki 11011^ +K2, (4.20) 
where the constants Ki and K2 are given by 

-.^ AcsCcWKlWp, (4.21) 
K2:=| |k||^ + c,c^(||b;il^,p + c?||&^f_i_^^^ + c^||«,^||?_^). (4.22) 

Introducing this expression in Eq. (|4.19p . one finds that for all 4>o constant and all 
0, both in [01, (f>2], it holds that 

Now, evaluate this expression at 0o = = 02 • The result is 

/p(02, w) ^ [{a^, - Ki) - a^i 02 - < 02^' - K2 02-'] Mi^2. (4.23) 
The assumption in Lemma [7] implies that — Ki > 0, so the polynomial 
g(02) (a^ - Ki) 0^ - 02 - 0^^ _ ^-7 

has the same form as the polynomial q introduced in Eq. (|4.15p . Therefore, the 
analysis done on the polynomial q in the proof of LemmalHlholds for the polynomial 
q, in particular, q has a unique positive root 0o. The remainder of the proof involves 
finding an upper bound for 0o, and the argument is almost identical to the one given 
in the proof of Lemma 

In the case that 0o ^ 1 (which could be verified, for example, by explicit evalu- 
ation), then choose a candidate for super-solution to be 0-|_ = 1. In the case that 
00 > 1, then there exists an upper bound for this root, as can be seen from the 
following argument. Given any 0o ^ 1, then the following inequalities hold, 

(0o)"+i>i ^ -(0o)-("+i) ^ -1 ^ -(0o)-" -00, 

therefore, this inequality for n = 3 and n = 7 implies that for all 0o ^ 1 holds 
g(0o) ^ S(0o) := (a;^ - Ki) (0o)-^ - (a^ + + K2)0o. 
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This new polynomial s vanishes at the point 



1/4 



which belongs to the interval where the polynomial q is increasing. Therefore, we 
have that 

q{4>i) ^ s(0i) = 

So, in this case (ftg > I, choose the candidate for super-solution to be (j)^ — (pi. 
Then, introducing the constant 

7 r, r 5^ + < + K2 ll/4- 
(p+ := max-^ 1, 

we have established that the following inequality holds 

fA<p2,w){^)^Q y^eWl^l, ycj)2^h. (4.24) 

and for all vector field we W^'^ solution of the momentum constraint Eq. (|2.36|) 
with source function G [0, (/)2]. 

As in the proof of Lemma [SI what remains is to account for the boundary con- 
tributions. The definitions of and K"^ imply that for any constant 02 holds 

-0iv],tr^^)^ ^ (i^>2 -<^^^) |M|i,2 y^^Wl^l- 
In particular, defining the constant 



1/4 0^ 



follows that (f>^ ^ 0+ and the following inequality holds 

([X02-0iv],tr^^)^ ^0 V(/)2^0+. (4.25) 
Adding Eqs. (|4.24p and (|4.25p we conclude that 

and for all to £ W"^'^ solution of the momentum constraint Eq. (|2.36p with source 
function (j> g [0,02]- Recalling now that any constant 02 satisfies (V02, V(^) = 0, 
we conclude that 

Ai^cj)2[<£) + f{<i>2,w){^)^Q y^eWl^l, V02^0+. 

Finally, employ now the constant 0+ as given by Eq. (|4.18p . This constant satisfies 
the conditions {4>d — 4'+)^ — and 0+ ^ 0^, so the following inequality holds, 

{<l>o - 0+)+ e Wl^l and + /(0+, w)] e W-'^^ 

for all 10 S PF^'^ solution of the momentum constraint Eq. (|2.36p with source 
<j) S [0,0+]. This establishes that 0+ is a global super-solution of Eqs. (|2.35[) - 
(1236)) . _ □ 

Consider now the particular case of a background metric h G C^{M,2) having 
a strictly negative Ricci scalar of curvature, that is, 

^-<<0 V(p e W'^i. (4.26) 

ll^l|l,2 

In this case it is possible to obtain a global sub-solution of Eq. (|4.9p . 
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Lemma 8. (Global sub-solution for R < 0) Consider the weak formulation for 
the Hamiltonian constraint given in j j^. J[ Assume that the constants , K^, 
and (f)"^ are positive, and the Ricci scalar R satisfies Eq. {4^.26\j . Denote by (f>- the 
constant 

Then, 0_ is a global sub-solution of Eq. 

Proof. (LemmalS^) We look for the sub-solution among the constant functions. 
Therefore, let 00 be any constant in [01, (/)2], with 01 > 0, then the following 
inequalities hold for every element (p G W^'^, 

(«T0o)*M + {aR(f>QT{'£) 
^ (ar,^)0o + {aR,'£)<Po 

where we used that both functionals a* and belong to the space W^^'^ , and the 
number 0o > 0. Introduce the polynomial 

9(0o) Or (t^l - flfl'/'O- 

There exists a unique positive root for q given by the number 0_ :— (a^/a^)^^^, 
and for all < 0o < 0- the corresponding values q{(t>o) a-re negative. Therefore, 
the following inequality holds 

fMMi^) ^0 e W^^i, V0O e (O,0_], Va; e W-]:\ (4.28) 

We now account for the boundary contributions. The definitions of the numbers 
0^ and imply that for any constant 0o, it holds that 

([i^0o-0^],tr^(^)^ (i^>o-<^]C)IMIi,2 y^ewl'l- 

In particular, defining the constant 

it follows that < 0_ ^ 0- and the following inequality holds 

([if0o-0iv],tr^^)^ ^0, V0oe(O,0_]. (4.29) 
Adding Eqs. (|4?28l) and (|4?29)) we conclude that 

{K^o,trN^)N + fw{4>o){^) ^0 e Wai V0O e (O,0_]. 

Recalling now that any constant 0o satisfies (V0o, V(/3) — 0, we conclude that 

AlM^) + UiMi^) ^0 e W^i V0O e (0, 0_]. 

Finally, introduce the constant 0_ as given by Eq. (I4.27[) . In particular, this con- 
stant satisfies (^u — 0_)^ — and < 0- ^ 0_, so the following inequality 
holds, 

(0Z, - 0-)- e W^r,'^ and - [Ai0_ + /^(0_)] e T4^4'^ Va^ G 
which establishes that 0_ is a global sub-solution of Eq. (|4.9p . □ 



A ; t) V ; 
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In the case that the Ricci scalar of curvature is non-negative, then it is not 
clear whether a constant and positive sub-solution to Eq. (|4.9|) exists. The latter 
exists when the conformally rescaled matter energy density p satisfies the condition 

> 0. This result is summarized in the following two Lemmas below. 

Lemma 9. (Global sub-solution for R ^ and > 0) Consider the weak 
formulation for the Hamiltonian constraint given in W4-l\ Assume that the constants 
a^, a^, (f)"^ K^, and (j)'^ are positive. Let cj)^ he the constant 



Then, cf)^ is a global sub-solution of Eq. EPI). 



-,1/4 



(4.30) 



Proof. (Lemma\^) We look for the sub-solution among the constant functions. 
Therefore, let 4>q be any constant in [01,02] with (pi > 0, then the following in- 
equalities hold for every element </? g 

fwFiMi^) = {ar(l)l)*((£) + {aR(j)o)*(ip) - (ap0o^)*((^) - (a^tf)^'^)* (ip) 

^ (ar0i5)*M + MoTi^) - («p0o')*M 



05 + 00 - a; 0o 



P 



\Ml.2, 



where the used the assumption that the functional is non-negative. Introduce 
the polynomial 

g(0o) 0g-f a^0o-a; 0o^ 

which is a non-decreasing function, because its derivative 

g'(0o) = 5a^ 04 + a^ + 3a; 0o^ 

is strictly positive for non-zero 0o. Rewrite the polynomial q as follows, 

qiM ^ (l>o^ QpiM with gp(0o) 00 + a^0o - ftp- 

There exists a unique positive root 0^ of the polynomial qp given by 

1/4 



(4.31) 



4a>^ 



Then, the inequality in Eq. (j4.3ip and the non-decreasing property of q imply that 
the polynomial q satisfies 

g(0o)s^O VO<0oSC0p. 
We then summarize the discussion above saying that the following inequality holds 



V0oe(O,0p], 



(4.32) 



and for all a* G W , 



From this point forward, the proof is identical to the 



proof of Lemma [8l What remains is to account for the boundary contributions. 
The definitions of 0)(r and imply that for any constant 0o, it holds that 

([i^0o-0^],tr^(£)^ (i^>o-0^)lMli,2 y^ewl'l. 
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In particular, defining the constant := min {0p, ^^^j? it follows that < ^ 
(pp, and the following inequality holds 

([if0o-0iv],tr^^)^ ^0, V0oe(O,0_]. (4.33) 
Adding Eqs. and we conclude that 

Recalling now that any constant (po satisfies (Vt^o, V(/?) = 0, we conclude that 

A^M^p) + fUMi^) ^0 e W^i V0O e (0, 

Finally, introduce the constant (/)_ as given by Eq. (|4.30p . In particular, this con- 
stant satisfies — 0-)^ = and < ^ so the following inequality 
holds, 

i<i>n~^-r & T^A'' and - [Ar.<f>_ + fU^-)] G W'l^^ Va^ G Ty^^'', 

which establishes that (p- is a global sub-solution of Eq. (|4.9p . □ 
Again in the case that the Ricci scalar of curvature R is non-negative there 
exists a sub-solution to the Hamiltonian and momentum constraint Eqs. (|2.35|) - 
(|2.36p . when the trace- free divergence-free two- index tensor a is big enough. By 
big we mean that > ctq, with the positive constant ctq given in Eq. (|4.35p . This 
result requires the near-CMC hypotheses present in Lemma [3 and is summarized 
below. 

Lemma 10. (Global sub-solution for i? ^ and > cto) Consider the 
weak formulation for the Hamiltonian and momentum constraints given in 
and assume that the hypotheses in Lemma^hold. Assume that the constants a^, 
0^ K'^ , and (f)'^ are positive, while the constant > a^, with the positive constant 
(To given in Eq. {4^.35^ . Denote by (f>cr the only positive root of the polynomial 
q(j{x) :— a'^x^ + a^x^ — a^/A, where a; G M. Let (j)- he the constant 

0- :-min{0., ^, 0^}. (4.34) 
Then, (/)_ is a global sub-solution of Eq. \2.35]) . 

Proof. (Lemma \1(A) Consider the weak formulation for the Hamiltonian and 
momentum constraints given in i J2.3l The defined in Eq. (|2.24p belongs to the 
space W^j^ , since the function a^, = (tr + Cwf' /% belongs to the space LP/^ with 
p — 12/5. Given any positive number e, the inequality 2|(Tafc(£w)°^| < eCT^ + (£to)^/e 
implies that function satisfies the following inequality, 

8a„ = (T^ + {Cwf + 2aab{CwY^ ^ [i ^ e) - {^^ - l) {Cwf , 

hence, for any number e G (0, 1) the functional a*^ must fulfill the inequality 

a;(^) = (a^, (^) ^ (1 - e) |M|i,2 - ^ - l) II^HIp IMIe V G W'/^. 

For every vector field to G W^'^ solution of the momentum constraint Eq. (|2.36|) 
with source function holds the inequality in Eq. p.26p . therefore there exist 
positive constants C£, ci and C2 such that 



\Cw\\l ^ 44 
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Hence, for all source functions e [0,0+], where the constant 0+ is the positive 
super-solution found in Lemma [3 holds the inequality 



\Cw\\l < ^cl 



b':\\K\\\,- 



+ cl\\wl,\t 



,p,lN 



4\\wid\\1,p 



Introducing the positive constant Ct defined by the inequality \\(p\\e ^ ct ||<p||i.2, and 
the constant ctq given by 



we obtain that 



cl WW 



Wli?-i,p,_sv + '^2 II^Jd||i,p 



(4.35) 



Choose the number e = 1/2, then we get a^ip) ^ (l/2)(a^ - (Jo/2) ||^||i,2- By 
assumption, we know that > ctq, therefore, we conclude that 



(4.36) 



Having established the inequality above we now start looking for the sub-solution 
among the constant functions. Therefore, let 0o be any constant in [01,0+] with 
< 01 ^ 0+ , then the following inequalities hold for every element (p £ Wjj_f_ , 

/u,f(0o)(^) = (ar0o)*(^) + (aii0o)*(^) - (ap0o^)*(^) ~ (a«"0o^)*(^) 

^ 0g + 00 - a; 0o ^ - J 0o ||^||i,2, 

where the used Eq. (|4.36p to obtain the last line above. Introduce the polynomial 

g(0o) := 0o + 0o - 0^ ^ - - i 

which is a non-decreasing function, because its derivative 

7 



/(0o) = 5a^ 0^ + a^^ + 3 a; 00"^ + - 0„-«, 



is strictly positive for non-zero 0o. Rewrite the polynomial q as follows 

1 w 



g(0o) = 00 g<T(0o) with go-(0o) '/'o + 01 



A 





(4.37) 



This polynomial has a unique positive root 0o-, which exists because Qa-iO) < 
and lim^g^oo '?(t('/'o) = while the root is unique because the polynomial is 
an increasing function for positive 0o. So, Eq. (j4.37p implies that q{(f>a) ~ 0, and 
together with the property of the polynomial q being non-decreasing, we conclude 
that 

(7(0o)<O VO<0o«;0<x. 
We then summarize the discussion above saying that the following inequality holds 



/»f(0o)(<^) e W^i, V0O G (0,0, 



(4.38) 



and for all vector field w £ W^'^ solution of the momentum constraint Eq. (|2.36p 
with source function e [0, 0+], where the constant 0+ is the super-solution found 
in Lemma[7l From this point forward, the proof is identical to the proof of LemmaEl 
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What remains is to account for the boundary contributions. The definitions of 0](r 
and K'^ imply that for any constant </)o, it holds that 



In particular, defining the constant (f)_ := min|0(j, 'J^}^ then follows that 
< 4>- ^ (f>a, and the following inequality holds 

([if0o-0iv],tr^^)^ ^0, V0oe(O,0_]. (4.39) 
Adding Eqs. and we conclude that 

Recalling now that any constant (t>o satisfies (V0Oj Vt/s) — 0, we conclude that 

AlM^) + UiMi^) ^0 e W^i V0O e (0, 

Finally, introduce the constant as given by Eq. (I4.34p . In particular, this con- 
stant satisfies {(po — 0-)^ = and < ^ so the following inequality 
holds, 

G ly^-' and - [Ai^cj)^ + e W^l^, 

and for all vector field w G W^''' solution of the momentum constraint Eq. (|2.36p 
with source function e [0, where the constant (f)^ is the super-solution found 
in Lemma U\ This establishes that 0_ is a global sub-solution of Eq. p.35p in the 
interval [0, i^-i-]. □ 

4.3. Summary on barriers. We present in this short Section a summary of the 
various results we have obtained in ^A.2\ for weak sub- and super-solutions. We 
state the constant sub- and super-solutions for each value of the Ricci scalar R. 
We do not state again the definition of the various constants that define the sub- 
and super-solutions, which can be found at the beginning of §4.21 Similarly, we 
do not state again the assumptions on the coefficients and data in the Hamilton- 
ian and momentum constraint equations needed to construct the barriers. These 
assumptions can be found in H4.21 or they can be read out directly from the barri- 
ers expressions, because they are the sufficient conditions that guarantee that the 
barriers are positive and finite numbers. 



:~ max < 1 



R bounded 



i? < 



1/4 (j)' 



?_ := mm 



max 1 , 



■K2 



1/4 (f)^ 



R^O 



{(I 

0^}ifa;^>O, 



} if < > 0, 



>0, 



Regarding the sub-solution for Ricci scalar i? ^ 0, only the constant (f)p has been 



given explicitly by the expression 



[-a- + ^(a-y + Aa-a-) /(2<) 



1/4 
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while the constant (pa has not been given exphcitly, but it has been proven that (pa 
is a finite, positive number. 

4.4. A priori L°°-bounds on Ty^'^-solutions. We now estabhsh some related 
a priori L°°-bounds on any VF^'^-solution to the Hamiltonian constraint equation. 
Although such results are standard for semi-linear scalar problems with monotone 
nonlinearities (for example, see |32j). the nonlinearity appearing in the Hamiltonian 
constraint becomes non-monotone when R becomes negative. Nonetheless, we are 
able to obtain a priori L°°-bounds on solutions to the Hamiltonian constraint in 
all cases including the non-monotone case. The results are based on a new abstract 
result (Lemma [24] in the Appendix) which holds for general semi-linear problems 
in ordered Banach spaces, under very weak assumptions on the nonlinearity (see 
the second assumption in parts [i] and [u] in Lemma [24]) . Monotone nonlinearities 
have the required property, but the property is much weaker than monotonicity 
and is satisfied for more general nonlinearities such as the one appearing in the 
Hamiltonian constraint. The results here generalize the a priori L°°-bounds on 
weak solutions appearing previously in an earlier set of unpublished notesQ and in 
a thesitH to the weakest possible assumptions on the coefficients appearing in the 
Hamiltonian nonlinearity. 

Theorem 6. {A priori bounds on W'^'^-solutions) Consider the weak formu- 
lation for the Hamiltonian constraint given in j j^. j[ and assume the hypotheses in 
Lemma in Lemma [3 and either in Lemma or in Lemma \1(A Assume that 
given a functional G Wj^^'^ , there exists a solution (j> G W^'^ of Hamiltonian 
constraint Eq. {4^.9^ . Then, there exists positive numbers (f)\/ , 0a with 0v ^ 0a ; 
such that 

< 0v ^ 0A, a.e. in M. (4.40) 

Proof. (Theorem\^) Let (p-^ and (pf^ be the (constant) sub- and super-solutions, 
respectively, given in Lemmas [MTOl In the proofs of these Lemmas it is established 
that the function fyj given in Eq. (j4.7|) is monotone increasing for its argument 4> 
satisfying ^ 0^ and for < ^ 0v Therefore, Lemma [24l implies that Eq. (|4.40|) 
holds. This establishes the Theorem. □ 

4.5. Results using variational methods. The Hamiltonian constraint Eq. ()4.9p 
can be written as the Euler condition for stationarity of a real-valued functional in 
a Banach space. Direct methods in the calculus of variations can be used to find 
the points that minimize this functional in particular types of closed sets in Ba- 
nach spaces: closed sets under weak convergence, which we denote here as closed^, . 
When the closed^; set is chosen appropriately, and the barriers found in §4.21 are 
incorporated into the argument, we can show that the minimum is actually a solu- 
tion of the Euler condition for the functional, that is, of the Hamiltonian constraint 
equation. The solution found with this approach requires fewer regularity assump- 
tions on the data, and the resulting solution is weaker (has less regularity) than 
the solution found using barrier methods. The variational structure was exploited 

^ M. Hoist, Weak solutions to the Einstein constraint equations on manifolds with boundary. 
Notes from the 2002-2003 Caltech Visitors Program in the Numerical Simulation of Gravitational 
Wave Sources. 

J. Kommemi, Variational methods for weak solutions to the Einstein Hamiltonian constraint 
on finite domains with boundary. Honors Thesis, Department of Mathematics, UC San Diego, 
2007. 



ROUGH NON-CMC SOLUTIONS TO THE EINSTEIN CONSTRAINTS 



43 



in [28j to develop an approximation theory and corresponding error estimates for 
numerical approximations to the Hamiltonian constraint. 

Let : W^'^ x W^;^ — > R be a bilinear form with action defined in Eq. (|4.ip . Let 
a*, a*n, a*, and a^, the functional defined in Eq. Let ^i, £ n VF^'^ 

be the positive, harmonic extension of the Dirichlet data function (jj^i discussed in 
and (f>lf be the Robin data functional, both data defined in Eq. (14.61) . Finally, 
let (f)i, (f)2 G L°° be functions satisfying < (jji ^ ^2, but otherwise arbitrary, and 
denote by U ■={[(f>i — 0d, 02 — 0o] H W^'^) C W^''^. Then, introduce the functional 

Jl-.U dW^'^ Mip) ■.^^aL{ip,ip)+GU^), (4.41) 

where the functional is given by 

Gw{ip) ■■= 9w{(t>D + - (i^lri^fNV) + aLi<pD,ip), 
with the functional gw{4') having the form 

gM) I K'/'' )*(</>) + 1 + 1 {apr''r{4>) + \ {a^r'r{4>); (4.42) 

while the definition M := M U {+00} U {—00} is explained in the Appendix. 

Theorem 7. (Existence of a minimizer) Consider the functional : U C 
Wl)^ — > M defined in Eq. \J^.4-l\j . Fix a positive extension (fyn G L°° n W^'"^ of the 
Dirichlet boundary data, and fix the Robin boundary data (j)*^ £ 2 •^(^A^jv, 0). 
Fix the functionals a*, a*p, a*^, and the function K satisfying the inequalities > 0, 
^ 0, ^ 0, K'^ > 0. Then, there exists an element G U minimizer of the 
functional Jl in U , that is, 

Jl{^) ^ inf Jl{^). 

Furthermore, if the Ricci scalar R is non-negative in M, then the minimizer ip is 
unique. 

Proof. (Theorem^) The set U is closed in W^-"^, as the following argument 
shows: By contradiction, assume that there exists a sequence {(/"n} G U such that 
(j)n — > 00 in W^''^ but 00 ^ U. In particular 0„ — > 0o in and the inequality 
|||0n|| — ||0o||| ^ ||0ri — 0o|| implies that ||0„|| — > ||0o||- The assumption that the 
limiting element (j)o ^ U implies that there exists a set Aio C Ai with meas(Alo) 7^ 
such that 0o > 02 or 0o < 0i. We consider here only the first case, the second 
one is proven in an analogous way. Then, we have the following inequalities, 

UuWm = ||0n|l5^o + ll0«ll(A1\7Vlo) 
WMmo + ll0n||(A<\A<o)' 

where we used the notation ||0n||A^o •= ll0n||L2(Aio)- Then taking the limit on both 
sides 

hm ||0„||^ sc: ||02|lX,„ + ||0o|| 

< UoWmo + ll0o|l(M\>lo) = WMm, 

from which we conclude that lim„^oo ||0n|| < Il0o||i contradicting the assumption 
above that lim„^oo Il0n|| = ||0o||- Therefore, the set U is closed. In addition, the 
set U is convex, therefore it is shown in the Appendix that U is closed^ . 
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We now show that the funetional Jl is coercive. First notice that the condition 
01 > imphes that the functional with values gw{<P) is continuous and bounded in 
[01 J 02] C , therefore there exists a positive constant Cg such that 

\9w{(t>)\ ^ Cg, V0 e [01,02]- 

Second, the Robin term is bounded as the following calculation shows 



where e is an arbitrary positive constant, and ct is a positive constant that bounds 
the trace operator. Third, the Dirichlet term is also bounded because the bilinear 
form flj^ is bounded, then the following inequalities hold. 



Fourth, the following argument, similar to the one used in the proof of Corollary[Tl 
shows that aj^ satisfies a Garding inequality. Start with a bound on the Robin term 
in fli. 




Vo9 e W 



1,2 



aL{(l)D,'p) = (V0D, V<^) + {KtrN(f)D,trNip)N 

^ (l + Q||i^||oc)||0D||l.2||^||l,2 



These calculations show that 



CgW ^^9 + Ye Il'^^ll-i2,iv + (1 + l|/^l|oo)' 1101,1 



(^i:tr^^(^,tr^^^p)^^ s$ ||X||oo \\tr(p\\l 



^Ct\\K\\^M\\\\7^\\ 



then it is not difficult to derive the following inequality on a^. 



^ii^ii?,2-ii'/'ir-|iiifiiLii^ii'-^iiv^ii?,2 
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where Ca(e) := [l + ||iir||^/(2e)] Cs ||</>2||^ ^'i^d Cs is the positive constant in the 
imbedding L°° C L'^ . Finally, the fifth step is to put all these inequahties together, 

>{\-\^)\Ml2-Cj{el V^£C/CW^1•^ 

where Cj(e) = Ca{^) + Cc{e) . By choosing e positive and small enough the inequality 
above establishes that Jj^ is coercive. Therefore, Jl is proper, and is also trivially 
bounded below by — Cj(e). 

We now show that the functional Jl is Isc^j , and we do it term by term on 
Jl- We start with the term proportional to a^i^, f), which is Isc^, by the following 
three facts: First, the norm in a Banach space is a Isc^, functional (statement proved 
in the Appendix); Second, the compactness of the imbedding W^''^ C L^. These 
two facts together imply the following: Given a sequence {ipn} C W^'^ such that 
(fn ipo in Wl)^ , and so — > (/3o in i we have that 

(V^, V^) = ||^||?,2 - \Wf ^ liminf(||^„||2,2 - ||^„f ) = liminf ||V^„||2, 

n — >oo n — ^oo 

which establishes that this term is Isc^, . The third fact is that the remaining term 
in the bilinear form has the form 

which is continuous under strong convergence and convex (the latter because the 
function > 0), therefore Lemma 1121 in the Appendix implies that this term is 
IsCu, . These three facts then establish that the functional ^ ^^{(p, (p)/2 is Isc^ . 
We now consider the remaining terms in the functional Jj^ which are present in the 
functional Gw The terms in which are linear in the function ip are Isc^, because 
they are continuous under weak convergence, by definition of weak convergence. 
The nonlinear terms are gathered together in the functional g^, and all of them 
except the term (aii4>)*{4>) continuous and strictly convex. Therefore, Lemma fT2l 
implies they are Isc^, . The only remaining term, (ancf))* {(f>), where is not positive 
definite, is also Isc^, , since the imbedding W^''^ C is compact. The proof of this 
statement is the following calculation: Let {(/>„} C W^'^ such that 0„ ^ 0o in W^'"^ 
which then implies that 0„ — > 0o in L^- Then we have 

\iaHM*iM~MnncPn) \ = |(a«,(0g-0^))| 

= I (flfl, (00 + 4>n)i4>0 ~ 4>n)) I 

= (l|afl||oo \\4>o + 4>n\{j 1100 - 0„|| -> as n ^ cx), 

which then establishes that the functional (f> i-^ [aa(j))* {(j)) is continuous under weak 
convergence, and so Isc^, . We conclude that the functional Jl is Isc^ . Hence, 
all hypotheses in Theorem [13] in the Appendix are satisfied by the functional Jl , 
therefore there exists a function ip E U minimizer of Jl- Furthermore, the minimizer 
if is unique in the case that the Ricci scalar R is non- negative in M. The reason 
is that in this case the functional Gw is strictly convex, since all terms in G^ are 
convex and at least one of them is strictly convex. Therefore, Theorem [T4l implies 
that the minimizer ip is unique. This establishes the Theorem. □ 
We now show that the functional Jl defined in Eq. (|4.4ip is Gateaux differen- 
tiable on U along directions in U. We also show that this derivative can be extended 
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along all directions on W]^^ , and that it coincides with the Hamiltonian constraint 
operator defined by the left hand side of Eq. ()4.9p . This technical Lemma is critical 
to connecting the minimizer of the functional Jl found in Theorem [7] to solutions 
of the Hamiltonian constraint, which we do below in Theorem [S] 

Lemma 11. (Jj, Gateaux difFerentiable) The functional Ji^ : U C W^''^ M 
defined in Eq. ^-41^ has Gateaux derivative DJi:^{(p)((p) for all (p ^ U along any 
direction Lp £ U . Furthermore, the map DJl{^p) : U ^ M. can he continuously 
extended for every ip G U into a map DJi:^{ip) : VF^'^ — > K, and this operator is 
precisely the left hand side in the Hamiltonian constraint Eq. j^-^P - 

Proof. (Lemma The Gateaux derivative of the functional defined in 

Eq. (|4.4ip can be computed term by term. By definition of the Gateaux derivative 
it is clear that 

DJl{(p){'£) = aL{^,'£) + Dg^ {(pD + <f) M ~ (j)*N M + , ^) 

= aL{4>,^) + Dg^{(f)){<£) - (i)*N{'£), 

where we introduced the notation (p :— (l)r> + if- By the definition of the functional 
5u, given in Eq. (|4.42p , and by the definition of the Gelfand triple structure described 
in the Appendix, it is possible to compute Dg^ term by term. For example, this 
calculation on the first term is the following: Denote gT{4>) •= {^/^){o-t4'^)*{4'), 
then 

i?5r(0)M = i lim I[{ar{4> + t^fn^ + t^)~{ar4>^n4>)] 
— b t^o+ r — — 

= ^ lim lim ^ [( am, {(t> + tip>f) - (arn,(f>'^)] 

= i lim lim \ [{am, {4> + t'f)'^) ~ (a-Tn, 4''^)] 
= lim (am4>^,V') 

n— i-oo — 

where 4> = (j^o+'f, and this calculation holds for all ip,(peU C Wd' . The limits can 
be interchanged to obtain the third line in the equations above because the sequence 
(aTn,0^) is uniformly bounded in the index n for every element (j> S 02], that 
is, 

\{am,4>')\ < mil \{arn,4>)\ < WUl IKII-O ||0||l,2. 

We need the assumption < 0i ^ (/)2 in order we can do the same calculation above 
for the terms in the functional g^ that contain negative powers of the function (j). 
Then, the principle of uniform boundness can be extended from sequences of linear 
functionals to the sequence that approximates the nonlinear functional g^ above. 
Therefore, the limits in the second line of the expression for Dgr can be interchanged 
to obtain the third line in that expression. See ^20j, pages 52-53, and also see |46) . 
pages 80-81 for a proof of the principle of uniform boundness. References about 
generalizations of this principle can be found in 20], page 82. Let us return to the 
proof of Lemma 1111 The expression on the last line in the inequalities above can 
be continuously extended for all ip E W^''^, as the following calculation shows 

{ar4>'r{ip) < ml KMi ml iKii-1.2 iMii,2. 
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Analogous calculations on the remaining terms in the functional then show that 
Therefore, we have established that 

for all G [/ and all (p e W^j"^ . This last equation establishes the Lemma. □ 
So far the functions (pi and 4>2 that define the subset U in Theorem [7] and 
Lemma fTD can be any elements in with only the condition that Q < (pi ^ <p2- 
Theorem [7] above says that there always exists a minimizer if of the functional Jl 
in the set U. The following result says that if the functions (f>i and (j)2 are sub- 
and super-solutions of the Gateaux derivative DJj^ : U C ^ W^^'"^ , then the 
minimizer is actually the solution of the Euler equation DJ^i^Lp) — 0, and is thus a 
weak solution to the Hamiltonian constraint Eq. (|4.9p . 

Theorem 8. (Hamiltonian constraint) Assume the hypotheses given in Theo- 
rem^ and also assume that either the constant > or the constant > ctq, 
where the positive constant do is defined in Eq. {^.35^ . Furthermore, assume that 
the subset U C W^'^ is defined by (j>i = (p- and cp2 — 4>+) where (p^ and (p-\- are 
any of the sub- and super-solutions of the Hamiltonian constraint Eq. \4-S^ found 
in Then, the minimizer if E U found in Theorem is a solution of the Euler 

equation 

where DJ^ is the Gateaux derivative of Jl, and the equation above is the Hamil- 
tonian constraint Eq. 

Proof. ( Theorem O ) Let </>_ , </)+ be sub- and super-solutions of the Hamiltonian 
constraint Eq. (|4.9p . respectively, and define U := — (po, (p+ — (po] H W^^' ^. Let 
(/? G C/ be a minimizer of Jl on U, whose existence was established in Theorem [71 
and denote <p := (pn-'r'^- We first establish the following result involving a minimizer 
ip of the functional Jl- Given any ^ £ Wl)'^ such that ip + tijj G U for small enough, 
positive number t, the following inequality holds, 

DJL{ip)W^O. (4.43) 

For the proof, compute the Gateaux derivative of Jl at if along ^p, 

DJLiipm = lim i [jLiifi + #) - Jl{^)] , (4.44) 

which is well-defined because we assume that for < t small enough the element 
Lp + t^j G U . The function ip is the minimum of the functional Jl in the set U , 
so Jl(<<2 + tip) ^ JLi'p), which establishes Eq. (|4.43p when the limit t ^ 0+ is 
computed in Eq. (|4.44p . 

We now use the inequality (|4.43p to show that the function p is solution of 
the Hamiltonian constraint Eq. (|4.9p . Let C be any scalar function in the space 
{Ai , 0) , and then introduce the mono-parametric family of functions S(pe as a 
perturbation of the function ip inside the set U, 

Sipe := min{(0+ - (/)o), max{((?!)_ - (/)o), ^ + <} }, 
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where e is a positive, otherwise arbitrary real number. The construction above 
impHes that 5ip^ £ U, and also that the functional Jl is Gateaux differentiable at 
the function ip along the function (5(^e — ip, where the latter statement follows from 

= lim i [Ji(t5(^, + ~ Ji((^)], 

which is well-defined because the set U is convex. Therefore, we can choose the 
particular direction 7p = {Sip^ — ip) in Eq. (I4.43p . which implies 

Oi^DJr.i^)iS^e-v). (4.45) 

It is now be convenient to use the equivalent expression Sip^ — {^p + eC) + C ~ C*^ ; 
where we have introduced the cut-off functions 

C :=max{0,(^ + eC)-('/'+-'/'D)}, C - min{0, ((^ + eC) - - (^i,)}. 

Notice that C^, are non-negative continuous functions belonging to W^'^. The 
inequality in ()4.45ll implies 

D,hi^)iO ^ - [DJUm^) - DM^)iC)]. (4.46) 
e 

We now show that each term on the right hand side in the inequality above ap- 
proaches zero as e approaches zero. The first term on the right hand side in 
Eq. (|4.46p satisfies the following inequalities, 

DMip){C) = DM^)iC) - DM^+ - <f>o)iC) + DM^+ - <f>D)iC) 

> [DMip)^DM<j>+^cf>o)]{C) 

= - 0+)(r) + [UW - fM+)] iC), (4.47) 



where the property that 0+ is a super-solution of Eq. (|4.9[) was used to obtain the 
second line in the inequality above. We analyze the last inequality, term by term. 
We will need the subset 7W C defined as follows 

■.= {xeM:C > 0}. 

This definition implies that meas(A^') ^ as e ^ 0. Then, the term involving the 
operator satisfies the following inequalities, 

A^icP - (f>+){C) = (V[(^ - VC) + (i^tr^[0 - tr^C)jv 
= (V[0-0+],V[0-0++eC])^. 
+ {KtrN[<f) ~ (/)+], tr^lct) -<j)++ 

^ e (V[(^ VC)_^, + e (/^tr^.[0 - (^+],trjvC)9_^^, 

where the subscripts tVC^ and dA4% on the inner products mean the inner 
product on these domains, and where we used that = {(f> — (j>+) + eC on M'^ . In 
order to analyze the term with the functional f^, it is convenient to introduce a 
representation based on Gelfand triple structure, W^'^ C = [L'^]* C W^^'"^, 
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with the functions fwn £ and (j) ^ U. Using this representation it is not difficult 
to establish the following inequalities 

[fM) - ^(0+)] (O = Jim ([/^„(0) - /.„(0+)] , C)m' 
> - lim (|/»„(0)-/.„(<^+)|,r)^. 
^ - lim (|/t„«(0) - /«,n('/'+)|, (0 - 0+)) jw. 
-e lim (|/^„((/i) - /t«r,(0+)|,C) jw. 

^ -e lim (|/u,n(0) - /t««(0+)|,C) vf.!- 

Combining the inequalities obtained for the operator Ai^ and the functional /u,, 
and using them with Eq. (|4.47p . we obtain 

- lim (!/«,„((?!)) - /i„„(0+)|,C) v(« 

An analogous calculation can be performed on the second term on the right hand 
side in Eq. (|4.46p . and the result is 

^J^M(C.) ^ -e[(V[0-0_],VC)_^^ + (iftr^[0 -(/._], tr^C) 

+ lim - /^.n(0-)|,C) VI 

where now we have used the property that the function 0_ is a sub-solution of 
Eq. (j4.9|) . and have introduced the analogous subset M.^ C M. defined as 

{a; e : C > 0}. 

These last two inequalities used in Eq. (|4.46|) imply that 

DJ^{^){C)^o{t), ase^O, 

since both meas(A^'^) and meas(A4c) approach zero as e approaches zero. The same 
calculation must hold for — C, and therefore we conclude that DJ/,(</))(C) = for all 
C G C^{M., 0). This space is dense in the space W^,' , so we conclude that 

which establishes the Theorem. □ 

4.6. Results using barrier methods. We now use the barrier method to show 
that there exist weak solutions to the Dirichlet-Robin boundary value formulation 
for the Hamiltonian constraint equation. The barriers found in ti4.2l are used to 
modify the original Eq. (|4.9p into an equation with a monotone decreasing source. 
This allows us to construct an iteration which converges to a fixed point of a 
particular mapping, which is constructed so that the fixed point also solves the 
Hamiltonian constraint equation. The modification of the original Hamiltonian 
constraint equation is called here a shift of the equation, and imposes a restriction on 
the regularity of the equation coefficients. This shift is not needed in the variational 
method, which is a reason that variational methods are able to produce results with 
weaker regularity, requiring fewer assumptions on the data. 

The barrier method and the variational method are both constructive, in the 
sense that they provide an algorithm to construct the solution, which could be of 
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interest in numerical relativity. In the latter method, one can build algorithms based 
on gradient descent that can guarantee progress (descent) at each iteration. In fact, 
the most effective numerical algorithms for the constraints tend to be a combination 
of these two ideas: global inexact-Newton methods are basically highly-tuned fixed- 
point iterations that maximize their contraction rate, and which are "globalized" 
by enforcing descent in an associated energy functional [28 . 

We therefore include the barrier technique here to give the most complete picture 
of what can be shown using both techniques on compact manifolds with boundary. 

Theorem 9. (Hamiltonian constraint) Consider the weak formulation for the 
Hamiltonian constraint given in Assume that the following conditions hold: 

(i) The coefficients functionals a*, a*p and given in Eq. 114. 3\ l have the form 

OrM = (flr,^), al{(£) = (a^,,^), a*p{l£) ^ {ap,l£), M <£_ G Wo"^, 

where the functions a-r, ap and belong to with p — 3. Fix a positive 

extension of the Dirichlet boundary data, as discussed in and the Robin 

boundary data as follows 

(a) In the case that the Ricci scalar R ^ 0, then assume that either Op > or that 
> (To, where the positive constant ao is defined in Eq. ^.35\ ); In the case 
that the Ricci scalar R < 0, then assume that > 0, ^ and ^ 0; 
(Hi) Assume that the constants a^, K'^ , (f)'^, and 0^ defined in ^JT^are all positive. 
Then, there exists a function 4> G {4>-,4>w+\ H A^''^ C W^'"^ which is a solution of 
Eq. \4-^ j where (j)yj+ is the super-solution found in Lemma and is the sub- 
solution given in Lemma\^for the case R < 0, and given in Lemma\^for the case 
R^O. 

Remark. The proof of Theorem [9] begins by shifting the equation (14. 9p in an 
appropriate way without changing its solutions. Then it is shown that a solution of 
the shifted equation exists iff there exists a fixed point of a certain map. It is then 
established that this map is compact, and thanks to the shifting it is also monotone 
increasing. These properties establish the existence of a fixed point. The last step 
in the proof is to show that this fixed point is a solution of the original Eq. (|4.9p . 
Remark. We note that the weak boundary value problem for the Hamiltonian 
and momentum constraint equations introduced in §2.31 is well-defined for source 
functions satisfying Eq. (I2.22p . that is, Ot, ap, and a^j belong to L^^^. However, 
Theorem [3] above requires that these coefficients belong to L^^^. This extra regu- 
larity is needed to shift the Hamiltonian constraint equation. It is not clear if there 
exists a different shifting procedure that also works for for coefficients in with 
6/5 ^ p < 3/2. This issue is also present in [38], where rough solutions are found in 
\Yk,2 fpj. ^ ^ 3^2; here we are basically asking for higher Lebesgue index p instead 
the higher Sobolev index k in [35] to make the shift possible. 

Proof. (Theorem\^) The assumption dm)) is required in Lemma[5]for the existence 
of the (local) super-solution (1)^+ given in Eq. (|4.14p . The assumptions iflz))- ifm)) are 
required to find a sub-solution: In the case of Ricci scalar R < the global sub- 
solution is given in Lemma [51 In the case that the Ricci scalar i? ^ the global 
sub-solution is given in Lemma [5| in the case > 0, and is given in Lemma [TUl in 
the case > ao. Summarizing, in these cases there exists at least local sub- and 
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super-solutions for the Hamiltonian constraint Eq. ()4.9|) . which is sufficient for our 
needs here. 

The condition > in assumption ifm)) imphes that the constant k in Eq. (|4.2p 
is positive in both cases where dMo = and dMo 7^ 0; so, the operator is 
invertible. 

We now use the sub- and super-solutions , (j)w+ to restrict the domain of the 
functionals f^p and fw to the set C L? . The interval itself and the 

functionals fwPi fw are well-defined, due to the property < 0_ ^ 4>w+- Let 
a e L^/^ be the function given by 

a:=b(j)l^^ar + \aii\+5-^ap + 7-^ayj, (4.48) 

and introduce a function s G i^/^ g^gh that (s — a) G L^^^. Then, define the shifted 
operators 

Al : W''^ ^ W-'^\ ^i0M :=A^0M + (s0,^), (4.49) 

/::[0_,0^+]ci2^M^-^'^ /:(0)(^):=/„(0)(^)-(s0,^). (4.50) 

First, note that the operator is well-defined, since for all s G L^/^^ a^d 0, 
G the generalized Holder inequality (p — 1, pi — 3/2, P2 = Pa = 6, in the 

notation given at the end of §l.f p implies 

(S0,^) < ||s||3/2 II0II6 IMIe ||s||3/2 II0II1.2 |M|l,2, 

where we used the fact that the imbedding W^'^ — > is continuous with imbedding 
constant c > 0. The functional is also well-defined because the shift term can 
be bounded as follows 

{S(j),ip) < ||s||6/5 II0IIOO ll^lle c||s||3/2 II0IIOO ||^|ll,2. 

In fact, the first inequality shows that the shift on the functional is well-defined 
for the shift function s G L^^^. However, we have just seen that the shift on the 
operator A]^ is well-defined only for the shift function s G L^^"^ C L^/^. 

This operator A^ is invertible since A^ is invertible (due to the hypothesis K'^ > 
0) and since the function s is non- negative (see for example |25j for a proof). This 
shifted operator A^ satisfies the maximum principle, a result shown in Lemma [20l 
in the Appendix. Therefore, Lemma [H] in that Appendix shows that {AD^^ is a 
monotone increasing operator. 

Second, note that the function s satisfies (s — a) G L^^, which implies that the 
operator is monotone decreasing. The latter means that given functions 02, 
01 G [0_, 0^-1-] with 02-01 G i+ , the functional satisfies -[/^(02) - /^(0i)] S 
W^^^ . The proof of this property is the following: given such functions 02 and 0i, 
compute 

{m2) - = (A«(02) - /.(0i))(^) - (5[02 - 0i],^) 

= ([/«)f(02) - /«,f(0i)],^) - (s[02 - (/-l],^) 
= (a4(02)' - (01 )'],^) + (ai5[02 - 0l],^) 

- (s[02 - 0i],^) - (ap[(02)-' - (0l)"'],^) 
-(a^[(02)-'-(0i)-T,^). (4.51) 
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Now, the conditions < ^ 02 and 0i, 02 G [0-i0iu+] imply the fohowing 
inequaUties, 

4 

ihf - ('^l)' = (E('^2F(0l)'-^)(02 - 01) 

^5(0^+)^(02-0i), (4.52) 




(4.53) 



(4.54) 



These inequahties and Eq. (|4.5ip imply 

imh) - s] (02 - 01 ),^), 

where a is given in Eq. (I4.48p . The choice s £ L^^^ and s ^ a implies that 
(/:('^2)-/:(0i))(^) s=:0 V^Giy^t 

which establishes that is monotone decreasing. 

Having introduced the shifted operators and we now remark that a 
function £ [0_, 0^,+] fl W^''^ is solution of Ai^cf) + fw{4>) = iff is solution of 
A^cj) + f.^{<j>) = 0. So far we have the following structure: 

: [0_,0„,+] CL^^M/^^'^ 
I : W^^^ L^, 

where / is the identity imbedding, which is a compact map. Therefore, the operator 
: [0-,0„,+] CL'^ T^(0) -/ (A^)- 1/^(0), (4.55) 

is well-defined. Both the operator (A|^)~^ and the functional — are monotone 
increasing, therefore the operator T*, is also monotone increasing, a result that is 
proven in Lemma 1221 Furthermore, this operator is compact, because it is a 
composition of continuous maps and the compact imbedding / : W^ '^ — > (for 
example see [20j page 486, Theorem 4, and also see the imbedding Theorems in [1] 
chapter VI). We established that the functions 0_ and 0,^,+ are sub- and super- 
solutions of Eq. (14. 9p . respectively. Therefore, Lemma in the Appendix shows 
that these functions 0_ and satisfy the inequalities in the order given by L^, 

0- «^ Twi(f>-), (t)w+ ^ T^(0u,+). 

Since the order cone in is normal, all the hypotheses in Theorem [15] in the 
Appendix are satisfied. Thus, there exists G [0_, 0^,+] C a fixed point of T^. 
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We now show that the fixed point <j) satisfies that (p G [<j>-,(j)w+] H W^'^. This 
resuh is a consequence of being bounded in W^'^ . Indeed, given any function 

ip G [4>- , 4>w+\ we have that 

||t:M||i,2 - II - (AD-V:MIIi,2 ^ II/:MII-i,2. (4.56) 
RecaUing the definition of the functional that is, 

fw{f)i'£) = + (Oflt^,^) - (ap^"^ (^) - (ay,Lp-\i£) - 0^(triv^) - {sLp,i£), 

we have the foUowing inequahties, 

+ |0^(trAr^)| + (s,(^) 0^+ 

+ Ct||(^;r||_l_2^Ar + S^0u;+ ||^||l,2, 

where is defined in an analogous way as in Eq. (j4.12p . and Cj is a positive 

f^l|i,2,]V 



constant such that ||trjv</5||i 2 at ^ Ct||((5||i_2 for all ip S W^'^. Therefore, introducing 



the constant 

:= [a^ 4>l+ + (flfl + s^) + < + < + ct||(^;||_ i,2.iv] , 
we have the inequality 

sup — II II ~ ^ c™, 

which yields to the desired inequality 

||T^MIIi,2 ^A:^ V(^e 

with hjji] — Ci^Cyj. Therefore, the fixed point point cf) G (/lu,-!-] satisfies 

Mia = WT^mia ^ 

which establishes the property (f> £ <j>w+] H W^'^. Therefore, we can apply the 
operator on both sides of the equation <j) = T^{(j>), and then the fixed point 
function (j) satisfies both the shifted and the non-shifted Hamiltonian constraint 
equations. The latter establishes that the function is a solution of the Hamiltonian 
constraint Eq. ((ij)) . □ 

4.7. Regularity of solutions. The following result states that when the regularity 
of the boundary data agrees with the equation coefficients regularity, the solution 
obtained by barrier methods is actually more regular than is stated in Theorem 1^1 
Note that the Proposition [1] below does not apply in the case of the solutions found 
by variational methods. In this latter case the coefficients in the functional 
belong to W^^'^, so the bootstrap argument mentioned below does not apply. 

Proposition 1. (Regularity) Assume the hypotheses in Theorem\^ and in ad- 
dition assume that the boundary dAi is C^. Assume that the boundary data satisfy 
4>*i^{trMip) = ((/)]v, tr]v(^)jv for all ip € Wjj and the following condition holds 

Then, the function (j) G H W^'^ solution of the Hamiltonian constraint 

Eq. {l§ belongs to the space ^^2,(^/2)^ 



54 



M. HOLST, J. KOMMEMI, AND G. NAGY 



Outline of the Proof. ( Proposition [][) A proof can be based on linear elliptic 
estimates (see for example Theorem 9.11 in [25]) and a standard bootstrap argu- 
ment. □ 

5. Coupled system 

Here we combine the results for the individual constraints derived earlier to es- 
tablish a new non-CMC result for the coupled system. In §5.11 we use the global 
barriers found in M.'2\ to establish existence of non-CMC solutions to the cou- 
pled constraints through fixed-point iteration and compactness arguments directly, 
rather than by using the Contraction Mapping Theorem as was done in the original 
work of Isenberg and Moncrief in [30] . 

It is interesting to note that for the main result on the non-CMC coupled system 
in [30], the near-CMC condition on the trace of the extrinsic curvature is actually 
used twice: once to obtain the global super-solution, and a second distinct time to 
construct a contraction for using the Contraction Mapping Theorem to get exis- 
tence and uniqueness. Here, a weak version of the near-CMC condition must also 
be employed in §4.2l to drive a global super-solution for the Hamiltonian constraint 
in our weaker setting. However, by using a compactness argument for the coupled 
system in il5.ll rather than the Contraction Mapping Theorem, we avoid the sec- 
ond use of the near-CMC condition. If a global super-solution can be constructed 
without the near-CMC assumption, then our compactness argument would give 
existence of solutions to the coupled system in the fully general "far-from-CMC" 
case. What our proof technique gives up is uniqueness of solutions to the coupled 
system, which comes for free with existence when the contraction argument is used 
as in [50] . 

5.1. Existence of weak solutions. This section is dedicated to establishing ex- 
istence of solutions to the weak Dirichlet-Robin boundary value problem (|2.35|) - 
()2.36p . In the case that R < there is no condition on the matter fields other than 
p ^ 0, but in the case that the Ricci scalar i? ^ it is required that either > 
or > (To, with the positive constant ctq defined in Eq. (|4.35p . The equation co- 
efficients are required to have stronger regularity than those previously required in 
Secs. l4.1 land l4.6l Although our problem formulation is different (compact domains 
with boundary), our results can be viewed as extending the result in [30] to weaker 
solution spaces, and extending their result for R — —1 to scalar curvature having 
any sign. The work presented can similarly be viewed as extending the CMC results 
on rough solutions in [38] to the non-CMC case, for compact domains with bound- 
ary. As remarked earlier, the "near-CMC" assumption required for the Contraction 
Mapping Argument in [30] is not required for the compactness argument below. 

Theorem 10. (Non-CMC) Consider the weak formulation for the Hamiltonian 
and momentum constraints defined in Assume the background metric h G 

C^(A^,2) and that the following conditions hold: 
(i) Fix a number p > 3 and denote by q := 6p/(3 +p). Fix source and boundary 
functions 

(l>D&W^^^, 0; eIy-^•P(aA^^^,0), w^^W^^\ w*^eW--^'\dM,M,l). 
with {p^ — j ■ j) G int{L'^^^) in the case j ^ 0; 
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(ii) In the case that the Ricci scalar R of the background metric is non-negative, 
then assume that either the constant a^p > Q or the constant > ctq, with the 
positive constant erg defined in Eq. {4^.35^ ; In the case that the Ricci scalar R 
is negative, then assume that a)^ > 0, Cp ^ and 0; 

(Hi) Assume that the function Ot and the constant Ki defined in Eq. satisfy 
that a!)^ — Ki > 0; also assume that the constants , (j)"^ , and (j)'^ are all 
positive. 

Then, there exists a solution 

6+p 

of the weak Dirichlet- Robin boundary value problem for the Hamiltonian and mo- 
mentum constraint Eqs. \2.35\) - ^2.8b]) . where is the super- solution found in 
Lemma^ and (j)^ is the sub-solution given in Lemma\^for the case R < 0, and 
given in Lemma\^for the case R^ 0. 

Proof. (Theorem\l(A) Notice that the definition of the numbers p and q satisfies 
that 3 < g < p. The assumption t £ L'^ indicates that Or — S L''/^, which 

imphes that the hnear functional a* given by as a* (if) — {or , tp) for aU ip e W^'^ 

is a well-defined element a* € W^^'^. The proof of the latter statement is based in 
the Holder inequality, which implies 

Wriv)] =^ l|ar||(g/2) ll^ll(g/2)'; 

since {q/2) = 3p/(3 + p), the relations + j^j^^y = 1 ^^'^ p ~^ ¥ ^ ^ imply 
that {q/2)' — 3p'/[3 — p']. Now the coefficient p > i implies that p' < 3/2, so we 
conclude that the imbedding W^'P C L^'/^) is continuous (see Corollary 7.11 
in §7.7), and then there exists a positive constant Cs such that 

\al{tp)\ ^ c, ||a^|l(,/2) ll^lli.p', 

which establishes that a* G W^^'P. The coefficient functions Op — np/A and = 
belong to L'^^, and the same holds for the coefficient acw — {Cw)^/8 whenever 
the vector w 6 W^''^. A similar argument as above shows that the functionals 

:= (op,^), al{ip) (a^,^), al^{^) := {acw^^p), G W^p 

are well-defined elements in W^^'P. The choice of the Dirichlet and Robin boundary 
data and the Gelfand triple structure reviewed in the Appendix imply that the 
functional / defined in Eq. (|2.27p is a well-defined map 

The assumption that the function j G L'^^^ implies that j G L^^'^ with (r/2) — 
3(7/(3 + q), since for g > 3 holds that 3 < r < g. Hence, the functional b*(u;) := 

(Kj,a;) for all a; G W^"^ is a well-defined element b* G W^^''. The proof is again 
based in the Holder inequality 

\b*{i^l)\ ^ 1^^11(^/2) l!^±^ll(r/2)', 

The condition g > 3 implies the inequality q' < 3/2 and the relations ^ + ^ = 1 ^i"i 
7r72T TrJW ~ ^ imply that (r/2)' = 3q'/[3 — g']. From the latter relation and the 
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inequality q' < 3/2 we conclude that the imbedding W^ *^ C L^''/^) is continuous, 
so there exists a positive constant Cs such that 



|6*(W)| SC Cs \\j\\{r/2) 



which establishes that b* G Wj^^''' . The assumption (Vt)* e W n^'"^ implies that 

br{ii£) (|Vr,a;) for aU a; G W]^^' is a well-defined element b*^ e The 
choice of the Dirichlet and Robin boundary data and the Gelfand triple structure 
reviewed in the Appendix imply that the functional / defined in Eq. (|2.28|) is a 
well-defined map 

/:[</)!, 02] CL°°-> W-^''\ 

The regularity assumptions on the equation coefficients and the assumptions ([np- 
Imj) are sufficient conditions to establish the existence of a global sub-solution 0_ 
for the Hamiltonian constraint Eq. (|2.35p . In the case that the R < Q this result 
is proved in Lemma [H and in the case that i? ^ this result is proved either 
in Lemma [9] in the case > 0, or in Lemma [10] in the case > ctq- The 
regularity assumptions on the equation coefficients and the assumptions ifm)) are the 
sufficient conditions to establish the existence of a global super-solution (f)+ for the 
Hamiltonian and momentum constraint Eqs. (|2.35p - (|2.36p . which was established 
in Lemma [71 

We now use these global sub- and super-solutions of Eqs. (|2.35p - (|2.36p to define 
the domain [(^-,0+] C L°° of the operators / and / given in Eqs. ([2?27)) - (|2?28)) . 
The inequality K'^ > in assumption fmj) implies that the operator Ajj, : VF^'^ — > 
Wo defined in Eq. (lOSl) is invertible, which was established in Theorem[3l The 
regularity result in Proposition [T] implies that the operator : W^'^ — > 
with the same action as defined in Eq. (|2.33p is also invertible. Regarding the 
Hamiltonian constraint equation, we introduce the same shifting done in the proof 
of Theorem [21 that is, fix a function s € L"^/^, given by 



s 



and introduce the shifting operators 
Al : W' P ^ W^''^, Alcj^i^) := A^cj,{ip) + {s^, ^), 

r ■■ [^Z-^-^^+l C X t^i'" ^ W-''P, r w){^) -.^ /(<^, w){ip) - is(b,^). 
Then, we have the following structure, 

(Al)-^ : W^^'P ^ A^ P C W^'P, 
I : W^'P L° 



- OO 



/: [0-,0+] Ci°°^ W^'^'^, 



where the map / : W^'^ L°° is the identity imbedding, which is compact for 
p > 3. Therefore, the following operators are well-defined, 

S : 0+] C ^ 5(0) -(A^)-i/(0), 

: [0_,0+] C ^ L^, - 1 {Al)-^ f , w) . 

Since we can choose in Theorem [21 to be the constant 0+ found in Lemma [3 
then this Theorem [51 and the regularity results in Proposition [H imply that exists a 
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fixed point G [4>-, (j>+] W^'^ of the iteration 

Therefore, the sequence {(/)", w"} given by (fp — (f)^, uP ~ Wj^, and 

is well-defined, where 0" S [(/)_, (/)-!_] fl W^^'P is a fixed point of the operator T^™, 
with n G N. By definition and by TheoremlHl each element in this sequence satisfies 
the equations 

+ /((/.",«;") =0, +/(</)"-!) =0. 

We will show that the sequence {(/>",«;"} C n VF^-p x W^'"^ is bounded. 

The proof is as follows. First, the elliptic estimates for the momentum constraint 
given in H'6.11 which imply, that there exists positive constants Ki, K2 such that 

Second, a calculation similar to the one performed after Eq. (|4.56p changing norms 
in W~^'^ with norms in W~^'P for p > 3 implies that 

||T:„(0")||oo^co||T:„(0")||i,p<fc„., 
k^^ := CqCl (t)\ + + S!;,) (f)+ + a'^ (pZ.^ 

where we have introduced the constants a^, a^, a^, and a^„, which are defined in a 

similar way as in Eqs. (|4.12p - (|4.13p changing the norms in W^'^ by norms in W^'^ . 
We have also introduced the number given by 

The dependence on n in kw^ is due to the term a^„. However, Eq. (|4.20|) implies 
that can be bounded for all n G N, and we obtain the inequality 

||T^„(^")||oo<fc0, 

fco cqCl [K <j)\ + (5^ + 0+ + h'; (j>Z^ 
+ (Ki^V' + K2) + coU*n\\~i,p,n] > 

SS.ax := 5 <f>lK + + 3 ^ a; + 7 || (Ki^V' + K2). 

This establishes that the sequence {</)", lo"} is bounded in W^-p x The latter 

space is a reflexive Banach space, so the sequence {</>", w"} has a weakly convergent 
subsequence, that is, there exist elements G W^'P and w G PF^'* such that 

(/)"^- ^ G M^^'P, k;"^ ^ w g W^^'*. 

The imbeddings W^'^ — > L°° and PF^'^ — > are compact since 3 < q < p, which 
implies that 

^"^ ^ (/) in vP^ ~* w in , 



58 



M. HOLST, J. KOMMEMI, AND G. NAGY 



that is, the convergence is strong in these spaces. We now first note that (p G 
[(/)-,(/)+], since the sequence (j/^^ ^ in the supremum norm and the interval 
is a closed set in Second, we now show that indeed 

lo"^ ^ w in W^^\ 

that is, the sequence {w^^} converges strongly to w in W^'"^ . The proof is to 
show that } is a Cauchy sequence in W^'"^ . This is shown by the following 
calculation, where we rename {0"^ , w"^} simply as {</>", lo"}. Then, we obtain. 



11^" -^™lli,,^ci||6;ii_i,, I -(rr 

5 

-1,9 



ciii&:ii_i,,||[^(0")^(0") 
<ciii6:ii_i,, [^ii^-iiLii-^' 



5-j 



m ||5-j 



in im 



3=0 



which leads to 



Since 0" is Cauchy in L°°, we have established that is a Cauchy sequence in 

The final step in the proof is to verify that cj) and w satisfy the constraint equa- 
tions l|235 | -(|236 | . Since G [0-,0+], the function T^(0) G W^^'^^ C is well- 
defined. What we have to show is that 



The first equation in (j5.1|) can be written conveniently as follows: 

r:(0) - = (r:(0) - Tixr)) + {Kir) - K.{4>n) + (0" - 

therefore, 

The first term on the right hand side above satisfies the inequalities 

= coii(Ai)-i(r('/',«')-r('^",^"))iii,p 

^coc^||r(0,t«;)-r(0",t.)||_i,p, 
and the last line above can be bounded as follows, 

iir (0, w) ~ nr, ^ iku,^ - {rnoo 

+ \\auup\\r'-ir)-'\u 



(5.1) 



(5.2) 
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Notice that 



and the same holds for 



relatfons together with Eqs. (|4.52p - (|4.54p imply that 

\\n<i>,w)-nr,w)\u..p 



a* and a^. These 



(5.3) 



Theorem ini requires s ^ a, with a given in Eq. (|4.48p . Choosing s = a we obtain, 

11/(0, w) - fir, s; U - 

with 

7- 



1 14 "ill- 



The number can be bounded independently of w since we have the inequality 
^ Ki (j)^ + K2, which is obtained in a similar way as the inequality in Eq. (|4.20p . 
just changing the norms used in that result to the appropriate norms needed here. 
Therefore, the following bound holds. 



with the constant kx given by 



5 ^Xa^ + 5^ + 3 ■ 



+ 7S^(Ki0V'+K2) 



(5.4) 



(Although we do not need to exploit this fact here, note that from the definition 
of kx it can be seen that there always exists source functions and boundary data 
small enough such that ^ fcr < 1, a condition that implies that this map is a 
fcy-contraction, as is it defined in [57 page 17.) 

The second term on the right hand side in Eq. (|5.2p satisfies the following bounds 



<co||r:(</)")-T:„(<^")||i,, 



^coCi(/)_ ||a^ 



Recall that the functional can be expressed as a^if) — (a^, ip) for all (p G Wjj'^ , 
with flu, — (a + CwY /%. We then conclude that e L"?/^, with q — 6p/(3 + p), 
which implies that for p > 3 we have the inequality 3 < q < p. Hence, there exists 
a positive constant Cs such that ||a^||_i^p ^ Cs||a. 

Now, the function be written as 

1 



mil (9/2) • So we have the inequality 

- aw'-\\(q/2)- 



- a^r. = - [(cr + Cwf - (o- + CvPf] 



- [2a + C{w + w'')] C{w - vP), 
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with each factor in L'^, so a simple case of the generaUzed Holder inequality (see 
the last part of §l.ip implies that 

||a^. -a^'. 11(9/2) ^ \\2(j + C{w + 11,'')% \\C{w-w'')\\q. 

Then, we have the further inequalities 

\\a^ - a^A\(,/2) ^ \ (2|k||, + \\C{w+w^)\\,) \\C{w - w'^)\\, 

<f [lkll, + C£ (K14+K2)] \\w-w''\\,,,. 
Denote k2 = coClCsC£0I [||(T||g + cc (Ki 0+ + K2)]/4, then 

-7^:(0")iioo < 11^" - him- 

Therefore, the inequality in Eq. (j5.2p implies 

\\n[<i>) - 0IIOO ^ fcT 110 - 0"||oo + k2 \\W-Uf'\\i.q + W - 0IIOO, 

and all the terms in the right hand side approaches zero when n approaches infinity, 
so we conclude that T^{(f)) = (f). 

We now show that the second equation in (|5.ip also holds using the following 
argument. We begin with the convenient representation 

5(0) -w = [S{^) - Sir)] + {w"+^ - w). 

This gives 

\m) - Hii,, ^ \m) - sir)\k, + 11^"+^ - hu.,- (5.5) 

The first term on the right hand side can be bounded as follows, 

11^(0) - s{r)\\i.g = \\-{A^)-\f{ci>)-f{r))l^^ 
11/(0) -/(</>") ILi,, 
<c^l|6;ii_i,,ii06-(0")6|u, 

where to get the last line we used the product property elements in L°° and elements 
in W^ '^ ^ which is discussed in the Gelfand triple part of the Appendix. Recalling 
the identity 

5 

and that 0, and 0" G [0_,0+], one finds 

11^(0) -5(0")||l,,^fc3 110- 0"|loo, 

with fcs — 6c2L II ''r 11-1,9 0+- Finally, inequality (|5.5p and the inequalities above 
imply 

11^(0) - Hli.g ^ fcs 110 - 0"||oo + ||t«"+' - Hli.g- 

The right hand side in equation above approaches zero as n approaches infinity. 
Therefore we conclude that 5'(0) = w. This result establishes the Theorem. □ 



ROUGH NON-CMC SOLUTIONS TO THE EINSTEIN CONSTRAINTS 



61 



5.2. Regularity of solutions. A bootstrap type argument shows that the regular- 
ity of weak solutions is actually related to the minimum regularity of the equation 
coefficients and of the boundary data. 

Proposition 2. (Non-CMC Regularity) Assume the hypotheses in Theorem \l(A 
assume that the boundary set dA4 is C'^ , and recall the parameters q = 6p/(3 + p) 
and p > 5. If the extension (po of the Dirichlet boundary data and the Robin data 
4'N for the Hamiltonian constraint equation 12.31]) satisfy 

then the solution (j) e n W^'^ and w G W^'' of the weak Dirichlet- 

Robin boundary value formulation for the Hamiltonian and momentum constraint 
Eqs. ^MS^-^MW found in Theorem [M also satisfies (j) £ [(t)-,cj)+] n V|/2,(9/2) 

Outline of the Proof. (Proposition^^) The proof can again be based on linear 
elliptic estimates (see Theorem 9.11 in [55^, and [T^, Vol. II, page 296) and a 
standard bootstrap argument. □ 

6. Summary 

In this article, we considered the conformal decomposition of Einstein's con- 
straint equations introduced by Lichnerowicz and York, on a compact manifold with 
boundary. We began by developing some basic technical results for the momen- 
tum constraint operator, and then established existence and uniqueness of W^'^- 
solutions to the momentum constraint (with conformal factor as fixed data) using 
variational methods. Among the technical results we established were generalized 
Korn inequalities for the conformal Killing operator on a compact manifold with 
boundary. Lemma 21 which does not appear to be in the literature. An alternative 
invertibility argument for the divergence of the conformal Killing operator is given 
using Riesz-Schauder theory, which yielded similar results in the case where the 
Dirichlet part of the boundary is non-empty. In both cases, the assumptions on 
the data were quite weak so that standard techniques cannot be used to establish 
additional regularity. 

We then considered the Hamiltonian constraint (with momentum vector as fixed 
data); using order cones in Banach spaces, we derived weak sub- and super-solutions 
to the Hamiltonian constraint. These can be viewed as non-trivial generaliza- 
tions of the barriers constructed previously in the literature to a setting with 
much weaker assumptions on the data. We also establish some related a priori 
L°°-bounds on any ly^'^-solution to the Hamiltonian constraint (TheoremlH]). Al- 
though such results are standard for semi-linear scalar problems with monotone 
nonlinearities (for example, see [32]), our results hold for a class of non-monotone 
nonlinearities that includes the Hamiltonian constraint nonlinearity and appear 
to be new. The generalized sub- and super-solutions are subsequently used to- 
gether with variational methods to establish existence (and uniqueness when scalar 
curvature i? ^ 0) of solutions to the Hamiltonian constraint in L°° n W^'^. Our 
arguments allowed the scalar curvature R to have any sign; the case of non-negative 
R required the additional assumption that either the matter energy density or the 
trace-free, divergence-free part of the extrinsic curvature be positive. Again, we 
made very weak assumptions on the data so that standard techniques cannot be 
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used to establish additional regularity of the solutions. Due to the lack of Gateaux- 
differentiability of the nonlinearity in the space W^''^, the connection between the 
energy used in the variational argument and the Hamiltonian constraint as its Eu- 
ler condition for stationarity was non-trivial, and was established through several 
Lemmas. Although our problem formulation is slightly different (bounded domains 
with matter) , the final result for weak solutions of the Hamiltonian constraint could 
be viewed as lowering the regularity of the recent result of Maxwell [38] on "rough" 
CMC solutions in W^^'^ for k > 3/2 down to L°°nW^^'^. We also gave an alternative 
non- variational argument using the more standard barrier methods, which requires 
more assumptions on the data, and yields essentially the Maxwell result for our 
problem formulation. 

We then combined the weak solution results for the individual Hamiltonian and 
momentum constraints to establish an existence result for the coupled system in the 
case of nonconstant mean curvature, through fixed-point iteration and compactness 
arguments rather than through the Contraction Mapping Theorem as used in the 
original 1996 work of Isenberg and Moncrief. This result requires more regularity 
than that needed for the results established for the individual constraints, with 
solutions for the conformal factor in W^'^ for p > 3, and momentum vector in 
W^''^ for q = 6p/(3 + p), but still extends the existing theory for the system in 
two ways. First, although our problem formulation is somewhat different (bounded 
domains with matter), the results could be viewed as extending the 1996 result of 
Isenberg and Moncrief on nonconstant mean curvature with Ricci scalar R ~ —1, 
to weaker solution spaces, and to cases where the Ricci scalar R can have sign. 
Second, again although the problem formulation is different, the result could be 
viewed as extending the recent rough solution work of Maxwell from the CMC case 
to the non-CMC case. 

It is interesting to note that for the main result on the non-CMC coupled system 
in [30], the near-CMC condition on the trace of the extrinsic curvature is actually 
used twice: once to obtain the global super-solution, and a second distinct time to 
construct a contraction for using the Contraction Mapping Theorem to get existence 
and uniqueness. By using a compactness argument directly rather than the Con- 
traction Mapping Theorem, we avoid the second use of the near-CMC condition. If 
a global super-solution can be constructed without the near-CMC assumption, then 
our compactness argument would give existence of solutions to the coupled system 
in the fully general "far-from-CMC" case. What our proof technique gives up is 
uniqueness of solutions to the coupled system, which comes for free with existence 
when the contraction argument used as in |30| . 

The variational approach used for the Hamiltonian constraint in the article 
should allow for the treatment of the case where the coefficient of the leading 
nonlinear term (p^ becomes slightly negative, by using the Mountain Pass approach 
as in the recent work of Hebey, Pacard, and PollackH The variational approach 
presented here also for the momentum constraint might make possible the combined 
variational treatment of systems involving the Hamiltonian and/or momentum con- 
straints as part of a large variational system. Finally, if the existing non-constant 
sub- and super-solutions in the literature for the Hamiltonian constraint can be 
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extended to our less regular setting, it would allow for weakening some of the as- 
sumptions on the signs of the coefRcients used for our results aboveH 

Although our presentation was for 3-manifolds, the results in the paper remain 
valid for higher spatial dimensions with minor adjustments, and the techniques we 
employed should extend to other cases such as closed and (fully or partially) open 
manifolds through use of techniques such as weighted Sobolev spaces. 
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Appendix A. Some tools from nonlinear functional analysis 

A.l. Gelfand triples. A readable reference for Gelfand triples is §17.1 in [51]. 
See also §23.4 in [SS], where they are called evolution triples. The vector spaces 
{X, H, X*) form a Gelfand triple iff the space H is a Hilbert space, X is a 
reflexive Banach space with dual space X*, and there exists a continuous imbedding 
I : X ^ H such that I{X) is dense in H . It can be shown that: li I : X ^ H \s 
continuous and I{X) dense in H, then the dual map /* : H* — > X* is continuous; 
in addition, since X is reflexive, I*{H*) is dense in X*. For the proof see [58] . 
page 417. Denote by R : H* H the Riesz map defined as follows: given an 
element h* S H* the element Rh* e H is given by h*{h) = [Rh* ,h)H for all heH. 
It can be shown that this map is a bijection, therefore it is invertible, and its inverse 
satisfies that for all ft, e the element R~^h G H* and the following equation holds, 
R^^hiK) = {h,h)H for all h€H.A Gelfand triple is usually denoted as 

X H = H* ^ X*. 

Gelfand triples are useful to study weak formulations of elliptic PDE. An example 
of a Gelfand triple is given by the Sobolev spaces {W^'^ , L'^ ,W~^-p ), with p' = 
p/(p-l). 

The property of a Gelfand triple we are most interested in is that I*{H* ) is dense 
in X*. This property together with the existence of the Riesz map imply that for 
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all X* G X* there exists a sequence C H such that the elements x* G X*, 

defined as x*^{x) := {hn,Ix)H for all x G X*, satisfy 

X* = lim <, in X*. (A.l) 

n — ^oo 

The elements x* can be written in terms of the Riesz map R and the imbedding / 
as follows, x* = I*R^^hn. The definition of the map /* implies that for all h* G H* 
holds 

rh* G X\ I*h*{x) = h*{Ix) Vx G X. 

Therefore, Eq. (|A.1[) can be expressed as follows: for all x* G X* there exists a 
sequence C H such that 

a;*(a;) = lim (hn,Ix)H Vx G X. 

n — >oo 

Definition 1. (Product) Let X — ^ H = H* X* be a Gelfand triple, and let 
V be a vector space such that there exists an imbedding X : V ^ H . Furthermore, 
assume that for all v ^ V and h G H there exists a map f , /i i— > {Xv)h G H , where 
the element {Iv)h in H satisfies that for every v G V there exists a positive constant 
c^ such that 

{{lv)h,h)^^c,\ih,h)H\ yh,heH. 

Then, given any x* G X* define the map v,x* (vx)* G X* as follows: 
{vxYix) := lim x*J{Jv){Ix)) VxeX. 

n — >oo 

An example of the situation above is the following: Let the Gelfand triple be 
given by the Sobolev spaces of scalar valued functions {W^'P,L^,W~^'P ), let the 
subspace V — L°^\ and let the map v,h vh be defined as pointwise multiplication 
a.e. in the domains of v and h. Then, all the properties in Def. [T] are satisfied. 
Indeed, the name "product" for the definition above originates in this example. 
This product given in Def. [T] satisfies the following property: 

WivxyWx' ^ Cv \\x*\\x' Vx* G X* and Vw G V. 
The proof is the following calculation: 

\{{Iv)h,,,{Ix))J^ 



{vx)*\\x' = lim \\{vxn)*\\x' = lim sup 



n — >oo 



Oyixex \\x\\x 



noticing that \(Xv)hn, {Ix)) ^\ ^ c„|(/i„, {Ix))h\ holds for all x G AT, then it also 
holds for the suprcmum in x, which leads us to the following inequalities, 

|(/i„,(/x)) I 

{vx)*\\x'^Cy hm sup — ^ Cv lim ||x*||x* 

"^°°o#xex mWx 



We now show that the map v,x* <—> (vx)* is well-defined in the sense that it is 
independent of the sequence x* that approximates x*. The proof is the following: 
Let {x^ } and {x~ } be sequences in AT* such that as rt ^ cx) holds 



Xr X m A 



Introduce the sequences {{vx)*^^} and {(wx^ )*} in A* such that 
(vxhj* ^ (vx)*, («%J* ^ (ra)* in X*. 
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Then, the following argument shows that 



^ e + c.„ Ija;;;;^ -il \\x* + 



e 



2e + Cy[\\xl^ - x*\\x* + \\x* - i^^^i,^.^ 
^2(l + c,)e, (A.2) 
where the second inequality comes from the following one, 

\lf \* ( ~ \*\\ \{v[hn~KiIx)J ^ 

Wv^xhj -y^hj lU' = — s^c^Wxj^ - x~ Wx". 

o=ixex \\x\\x 

Since e ^ as n oo, we then conclude from Eq. (|A.2p that {vx)* = (vx)* . 

A.2. Garding inequality and Riesz-Schauder theory. We recall now (without 
giving a proof) the well-known result of Riesz and Schauder. Standard references 
for this result are in [5S] page 76, in [FT] page 166, and in page 372. 

Theorem 11. (Riesz-Schauder) Let X be a Banach space, K : X ^ X be a 
linear and compact map, andix : X ^ X be the identity map. Then, the following 
statements hold: 

(i) d\TaN(Xx-K) = dim7V(IJ^_K.),• 
('^^j Given f G X the equation 

[Ix -K)x^f (A.3) 



has a solution iff x*^{f) ~ for all a;* G N(^x* - 



K' 



(Hi) //dimiV(-x^_x*) = 0, then the condition a;*(/) = is trivially satisfied for all 
elements f € X , hence for every f G X there exist a unique element x G X 
solution of Eq. IIA.3\) . Furthermore, the operator {Ix ^ K)^^ , whose existence 
is asserted here, is linear and bounded. 

(iv) If d\mN(x* -K') > and the condition x'^{f) — is satisfied, then the solu- 
tions X of Eq. liA.3\) are not unique, and given any solution x then 
is also a solution, with a;„ G -^(Ij^-k)/ 

We now use the Riesz-Schauder Theorem above to show whether a linear equa- 
tion involving a bounded bilinear form satisfying Garding's inequality has solutions. 
Let {X, H, X*) be a Gelfand triple, as it is defined in the previous subsection of 
this Appendix. Introduce a bilinear form a : X x X -^M. and consider the following 
problem: Given an element /* G X* find an element x £ X solution of the equation 

aix,x)^f*{x) \^xeX. (A.4) 

It is convenient to reformulate this problem in terms of operators instead of bilinear 
forms. Introduce the operator A : X ^ X*, with action Ax{x) := a{x,x) for all 
x,x G X. Then, the problem above has the following form: Given an element 
/* G X* find an element x € X solution of the equation 

Ax = /*. (A.5) 

It is also convenient to introduce the Banach adjoint operator A* : X ^ X* 
defined as A*x{x) :— Ax{x) for all x,x£X. We are identifying X with its double 
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dual space X** . Let Na, and Na* be the null spaces of the operators A and A* , 
respectively. 

Theorem 12. Let {X, H, X*) be a Gelfand triple, and in addition assume that the 
imbedding I : X —> H is compact. Let A : X —>■ X* be a linear, bounded operator 
satisfying Girding 's inequality, that is, there exist positive constants ko and Kq 
such that 

\\Ax\\x' ^ Ko\\x\\x, ka\\x\\l^\\Ix\\l+Ax{x), WxeX. 

Then, dimA^^ — dimA'^* and Eq. \A . 5\) has a solution x G X iff f*{xn) = for 
all in G Na* ■ Furthermore, the following statements hold: 
(i) If diva N A* = 0, then there exists a unique x € X solution of Eq. IIA.5\) for 
all f* e X* ; Furthermore, there exists a positive constant cq such that the 
following estimate holds, 

\\x\\x co\\Ax\\x' yxeX; (A.6) 

(a) If dim Na-' > and the condition f*{xn) — for all Xn G Na-' holds, then the 
solution X is not unique, since x' :— x + a;„ is also a solution, with a;„ £ Na. 

Proof. (TheoremMSX ) Given the operator A, introduce the operator Ax ■ X — ^ X* 
with action 

Axx{x) :— Ax{x) + {Ix,Ix)h- 
The assumptions that the operator A is bounded and satisfies Garding's inequality 
imply that the operator Ax is bounded and coercive, respectively, hence, invertible. 
Notice that Ax can be written in terms of operators as follows: Ax = A + J , where 
J : X ^ X* is given by J :— I*R^^I, since 

Jx{x) = I*R-^Ix{x) ^ R^^Ix{Ix) = {Ix,Ix)h. 

The Eq. (|A.5|) can be re-expressed as follows: 

Ax^f* ^ {Ax-J)x = f* ^ {Ix - A^\j)x = A-'f*, 

where Xx : ^ ^ ^ is the identity map. Introduce the notation fx A^^f* G X 
and the operator K : X ^ X given by K :— A^^J. So, x is solution of Eq. (IA.5|) 
iff it solves the equation 

iIx-K)x = fx. (A.7) 
Since the imbedding I : X ^ H is compact, and the remaining maps that define 
K are continuous, we conclude that K is compact (for example see [20] page 486, 
Theorem 4). Then, the operator Xx — K is a Fredholm operator of index zero, and 
Theoreminiimplies that dim Nt^^-k = dimiVx^-K* • By construction we have that 

Nxx-K = Na- One can also show that x* e Nj'^^k* iff in ■— (^x^)*a;^ G Na'- 
Due to Ax is a bijection, this shows that 

dim Na = dim Na* . 

Theorem [TT] implies that Eq. (jA.7|) has solution iff x*^{fx) — for aU x* € Nx*^-k* ■ 
This condition can be rewritten as follows: 

= xiifx) = xi{A~^f*) = r = r(5;„) Vi;„ e tv^., 

which is the condition appearing in Theorem [T^ In the case that dimTVA = 0, 
then dimA'^. = 0, and so the condition x'^{fx) = is trivially satisfied for all 
fx £ X. Therefore, Theorem [11] implies that for every element /* £ X* there 
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always exists a unique solution x E X oi Eq. (|A.5|) . This statement defines the 
operator A^^ : X* X, and Theorem [Tl] asserts that this operator is linear and 
bounded, the latter property implies that there exists a positive constant cq such 
that 

\\x\\x ^co\\Ax\\x' VxeX. 
This establishes part ^ in Theorem [T^l In the case that dim A''^ > and the 
condition x*(/x) = is satisfied, then Theorem [TT] says that a solution x G X 
exists, and also a solution, where x„ G Na- This estabhsh part ^ 

in Theorem [H □ 

A. 3. Variational methods. These notes follow the main ideas in Chapter 4 of 
Part Two in [33j, and §1, §2 in Chapter 1 in |3H]. An introduction into this subject 
is §7.1 in [39]. The main result of this Section is Theorem [T3l We could not 
find in the literature this result precisely in this form, needed for the Hamiltonian 
constraint problem, so for completeness we included the proof of the Theorem. 

Given a Banach space A", a subset ?7 C A" is called closed under weak conver- 
gence (closed^) iff for all sequence {x„} C U such that x„ ^ a^o in X holds that 
xo e U. Every closed^, set in a Banach space is closed, but the converse statement 
is not true. A particular class of closed sets that are also closed^ are closed and 
convex sets. Given a vector space V , a subset U C V \s called convex iff for all 
X, x G U the elements [tx + (1 — t)x] € J7 for i G [0, 1]. The proof of the above 
statement is based in a result by Mazur (see Theorem 2.2.4 on page 142 in [35]) 
that says: In a Banach space, for every sequence {xn} such that Xn xq there 
exists a sequence {y-m} such that ym xq, where the element y^ are constructed 
as a convex combinations of the a;„, that is, 

m m 

ym ■= ^ Kxn, (a„ ^0, ^ A„ = 1^ . (A.8) 

n— 1 n— 1 

Using this result is not difficult to show that every closed and convex set ?7 in a 
Banach space X is also closed^ , as the following argument shows: given {xn} C U 
such that Xn xq in X, use Mazur's idea to construct the sequence {ym} as a 
convex combination of the Xn such that ym — > xq in X. However, U is convex, so 
{Vm} C U, and it is also closed, so xq G U. This establishes that U is also closed^ . 

Let X be a Banach space, and introduce the functional J : X ^ R. The 
symbol M means that there might exist points xq G U such that there exists a 
sequence Xn xq with lim„^oo J{xn) = oo or equal — oo. If such points exist, 
then J is an unbounded operator with domain Dj strictly included in X. This 
idea is summarized with the notation M := R U {+00} U {—00}. The notation 
J : X ^ R is more convenient than the notation J : Dj C X ^ R in cases 
where the continuity or the differentiability of the functional is not important in 
the situation under study. An example is the problem of finding the local or global 
minimum of a functional using direct methods, which do not include computing 
the Euler equations for the functional. Given any subset U C AT of a Banach space 
X, the functional J : [/ C A" ^ R is called proper on U iff for all {a;„} C U 
such that ||cc„||x — * 00 holds J{xn) — > +00. A particular case of proper functionals 
are coercive functionals, where J is called coercive on U iff there exist positive 
constants cq, ci, such that for all x G U holds J{x) ^ cq \\x\\x — ci. The functional 
J is bounded below on U iff there exists ao G Hi such that J{x) > aa for all 
X £ U. All coercive functionals are bounded below by — ci. 
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The functional J : X ^ M is lower semi-continuous (Isc) at the element xq G 
X iff for all sequence {xn} C X with Xn —* xq in X holds J{xq) ^ liminf„^oo J{xn). 
Any continuous functional is Isc. The functional J : X — > R is called lower semi- 
continuous under weak convergence (IsCi^,) at the element xq ^ X iS for all 

sequence C X with Xn xq in X holds 

J(xo) ^ liminf J(a;„). 

n — >oc 

Given any set [/ C X a functional J : U C X ^ R is Isc (respectively Isc^u) 
on U if it is Isc (respectively IsCu,) on all points in U. Not every Isc functional 
is IsCiu . The later property is a stronger condition on the functional than the 
former property, due to the set of all sequences that converge weakly in a Banach 
space is bigger than the set of all sequences that converge strongly. An example 
of a IsCu, functional, mentioned in [48j , is the norm in an arbitrary Banach space, 
as the following argument shows: Let X be a Banach space, {a;„} C X be any 
sequence such that a;„ xq, then there always exists an element x*^ g X* such 
that xI^^{xq) — Wx^Jlx* ||a;o||x; then we obtain 

l!4olU* ll^^olU ^^Igixo) = liminf 4 (a;„) liminf ||a;„||x, 

and this implies that H^oUx ^ liniinf„_^oo ||a:^n||xj establishing our assertion. The 
well-known case of the norm in a Hilbert space H being Isc^, follows from the pre- 
vious argument choosing x*g{x) — {xo,x)h, where x is any element in H. 

Theorem 13. (Existence of a minimizer) Let X be a reflexive Banach space, 
and U C X be a closedw subset. Let J : U C X -^R be a proper, bounded below, and 
IsCw functional. Then, there exists an element xq e U minimizer of the functional 
J in the set U , that is, 

J{xo) = inf J{x). 

Proof. fTheorem\13[) The functional J is bounded below in U, therefore there 
exists a positive constant ao such that for all x G U holds J(x) ^ ap. Then, there 
exists a minimizing sequence, that is, a sequence C U such that J{xn) ao as 
n — > oo. The functional J is proper on U, therefore the sequence {Xn} is bounded. 
The Banach space X is reflexive, which implies that there exists xq G X such that 
Xn ^ Xq as n oo. The set U is closedui , therefore xq G U. Finally, the functional 
J is IsCw and the sequence {a;„} is a minimizing sequence, which imply that 

J(xo) ^ liminf J(xn) = ao ^ inf J(x). 

n — >oo x^U 

The definition of infimum implies J{xo) — inf J{x), which establishes the Theorem. 

x^U 

□ 

Let U d X he a. convex set in a Banach space X . A functional J : [/ C AT ^ R 
is called convex iff for all x, x^ U and t G [0, 1] holds 

J{tx + (1 - t)x) tJ{x) + (1 - t)J{x). 

A convex functional J is called strictly convex iff for all x, x £ U, with x ^ x, 
and t e (0, 1) holds 

J{tx -t- (1 - t)x) < tJ{x) + (1 - t)J{x). 

Besides these main Theorems above, the following Lemma is also needed in the 
proof of Theorem [T] 
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Lemma 12. Let X be a Banach space, U G X be a closed, convex set, and J : 
U <Z X ^ R be a convex and Isc functional. Then, the functional J is IsCw . 

Proof. fLemma llSl ) Let C U he any sequence such that Xn xq in X, 

and denote ao '■= hminf„_>oo J{xn)- Earher in this Section it was shown, using a 
Mazur's sequence, that a closed and convex set in a Banach space is also closedu, , 
therefore xq G U. Using once again Mazur's result, let {ym} C U he a convex 
combination of the elements x„ such that ym — > xq in U. Furthermore, let the 
convex combination of the elements Xn start at n — N for some number iV G N 
instead of n = 1, that is, 

m m 

ym = ^ ] A„a;„, ^ 0, ^ ] A„ = 1 

n=N n=N 

The functional J is convex, therefore, 

m 

J{ym) \nJ{Xn), (A.9) 

n=N 

By definition of the constant ao and after selecting a subsequence if necessary, given 
any positive number e there exists a number N{e) G N such that 

J(x„)<Q;o + e, \Jn^N{e). 

Then, Eq. (fO]) implies 



m 

J{yni)<{^ ^ A„^ (ao + e) = ao + e. 



n=N(e) 

SO by choosing e small enough we conclude that liminfm^oo J{ym) ^ Q^o, or alter- 
natively, 

liminf J(?/„i) ^ liminf J(a;„). 

m — >oo n — *oo 

Finally, recalling that the sequence ym xq and that the functional J is Isc, we 
have that J{xo) < liminfm-,oo J{ym), which together with equation above says, 

J(xq) ^ liminf J{xn)- 

n — >oo 

This equation establishes the Lemma. □ 
For completeness we now state the following result, which establishes the exis- 
tence and uniqueness of minimizers for certain type of convex functionals. 

Theorem 14. (Minimizers of convex functionals) Let X be a Banach space 
and U C X be a closed, convex set. Let J : U cX— >R6eo convex, Isc, 
and coercive functional. Then, there exists an element x d U minimizer of the 
functional J in the set U. Furthermore, if the functional J is strictly convex, then 
the minimizer x is unique. 

Proof. (Theorem\14\) We know that a closed, convex set J7 in a reflexive Banach 
space X is a closed^, set. And a convex and Isc functional J on a closed^ set U is 
IsCu) , result proved in Lemma [T2l Since the functional J is also coercive, then J is 
proper and bounded below. Therefore, Theorem [13] implies that the functional J 
has a minimizer x E U. 



70 



M. HOLST, J. KOMMEMI, AND G. NAGY 



Assume now that the functional J is strictly convex, and assume that there exist 
two minimizers xi, X2 (z U of the functional J, that is, 

Jo J{xi) = J{x2) = inf J{x). 

We will now construct a contradiction. Assume that the minimizers are different, 
^1 ^2, Siiid introduce the elements Xt '■— txi + (1 — t)x2, for t £ (0, 1). The strict 
convexity of the functional J implies 

J{xt) = J{txi + (1 - t)x2) < tJ{xi) + (1 - t)J{x2) = tJo + (1 - t)Jo = Jq. 

We then conclude that J{xt) < J{xi) = J{x2), contradicting the assumption that 
the elements xi, X2 are minimizers of J. Therefore the minimizer must be unique. 
This establishes the Theorem. □ 
We finish this Section with a calculation that is useful to verify whether a 
Gateaux diffcrcntiable functional is convex or strictly convex. Let B^{xo) C X 
be an open ball of radius e centered at the element xq E X . If a convex functional 
J : U C X ^ X is Gateaux differentiable in the convex set U with Gateaux de- 
rivative D J, then, there exists a positive and small enough number e such that the 
following inequality holds 

DJ{x)v s; J{x + u) - J{x) Vx e int([/), \/v G B,(0). (A.IO) 

The proof is the following calculation: Fix a positive number e, then given both 
X G mt{U) and v S Se(0) there exists a small enough number e such that x + v € U. 
The convexity of the set U and of the functional J imply that 

J{t{x + i;) + (1 - t)x) ^ tJ{x + v) + {l- t)J{x), 

which in turn implies 

i [J{x + tv) - J{x)\ s$ J{x + v)- J{x). 

Then, Eq. (jA.lOp follows by taking the limit i ^ 0+ in the inequality above. A 
similar proof establishes the following result: If the functional J is a strictly convex 
and Gateaux differentiable in a convex set U, then there exists a positive and small 
enough number e such that the following inequality holds 

DJ{x)v < J{x + v) - J{x) Vx G mi{U), Vw G B,(0), v^Q. (A.ll) 

A. 4. Ordered Banach spaces. These notes follow the main ideas and definitions 
given Chapter 7.1, page 275, in [57], while some examples were taken from [3] 
and [1^. Let X be a Banach space, be the non-negative real numbers. A 
subset C C AT is a cone iff given any a; G C and a G M+ the element ax Q C. A 
subset X^ C A" is an order cone iff the following properties hold: 

(i) The set X+ is non-empty, closed, and X+ ^ {0}; 
(ii) Given any a, b £ R+ and x, xd X-^- then ax + bxCz Ar_|_; 
(in) If a: G X^ and —x G X+, then x = 0. 

The second property above says that every order cone is in fact a cone, and that 
the set X+ is convex. The space X = is a convenient Banach space to picture 
non-trivial examples of cones and order cones, as can be seen in Fig. [T] A pair X, 
X+ is called an ordered Banach space iff X is a Banach space and X^ C AT is 
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an order cone. The reason for this name is that the order cone defines several 
relations on elements in X, called order relations, as follows: 



finally it is also used the notation u ^ v, u < v, and u <^ v to mean v u, 

V > u, V ^ u, respectively. A simple example of an ordered Banach space is R 
with the usual order. Another example can be constructed when this order on M 
is transported into C°(A^,0), the set of scalar- valued functions on a set Ai C K", 
with 71 ^ 1. An order on C°(A4,0) is the following: the functions u, v € C°{A4,0) 
satisfy u ^ v iS u{x) ^ v{x) for all a; £ A^. The following Lemmas summarize the 
main properties of order relations in Banach spaces. 

Lemma 13. Let X , be an ordered Banach space. Then, for all elements u, v, 
w & X, hold: (i) u ^ u; (ii) If u ^ v and v ^ u, then u = v; (Hi) If u ^ v and 

V ^ w, then u ^ w. 

Proof. (Lemma\l!A) The property that u — m = G X^ implies that u ^ u. 
li u ^ V and v ^ u then u — v £ X^ and —{u — v) G X^, therefore u — v = Q. 
Finally, ii v and v ^ w, then u — t; G Xj^ and v — w G which means that 
u ~ w — {u — v) + {v ~ w) G X^. □ 
Furthermore, the order relation is compatible with the vector space structure 
and with the limits of sequences. 

Lemma 14. Let X, X+ be an ordered Banach space. Then, for all u, u, v, v, 
w G X , and a, 6 G K, hold 
(i) If u ^ V and a ^ b ^ 0, then au ^ bv; 

(ii) If u ^ V and u ^ v, then u + ii ^ v + v; 

(Hi) If Un ^ Vn for all n CzN, then lim„^oo Un ^ lim„^oo ^n- 

Proof. (Lemma \14\ ) The first two properties are straightforward to prove, and we 
do not do it here. The third property holds because the order cone is a closed set. 
Indeed, u„ ^ u„ means that u„— w„ G X^ for all n G N, and then lim„^oo(M„— u„) G 
X+ because is closed, then Property dm)) follows. □ 
The remaining order relations have some other interesting properties. 

Lemma 15. Let X , X^ be an ordered Banach space. Then, for all u, v, w G X , 
and a G K, hold: (i) If u ^ v and v ^ w, then u ^ w; (ii) If u ^ v and v ^ w, 
then u ^ w; (Hi) If u ^ v and v ^ w, then u ^ w; (iv) If u ^ v and a > 0, then 
au ^ av. 

The Proof of LemmafTSlis similar to the previous Lemma, and is not reproduced 
here. Given an ordered Banach space X, X^, and two elements u ^ v, introduce 
the intervals 

[v, u] :— {w £ X : V ^ w ^ u}, (v, u) :— {w G X : u ^ w <C u}. 

Analogously, introduce the intervals [v,u) and {v,u]. See Fig. [T] for an example in 
X = . Useful order cones for solving PDE are those that define an order structure 
in the Banach space which is related with the norm and the notion of boundness. 
These type of order cones are called normal. More precisely, an order cone X+ in a 
Banach space X is called normal order cone iff there exists < a G M such that 
for all u, V £ X with < w < u holds ^ a \\u\\. 



u ^ V iS u — V Cz X^ , 
u^viSu — vCz int(X+), 



u > V iS u ^ V and u ^ 
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Figure 1. The shaded regions in the first picture represents an 
order cone, while the second picture represents a cone that is not 
an order cone. The shaded region between u and v in the third 
picture represents the closed interval constructed with the 

order cone R^, which is also represented by a shaded region. 

Lemma 16. If X, is an ordered Banach space with normal order cone X^, 
then every closed interval in X is bounded. 

Proof. fLemma \lb\ ) Let w £ then v ^ w ^ u, and so ^ w ~ v ^ 

u — V. Since the cone X^ is normal, this implies that there exists a > such that 
\\w — v\\ ^ a \\u — v\\. Then, the inequalities \\w\\ < ||w — ?;|| + ||?;|| ^ a \\u — v\\ + \\v\\, 
which hold for all w E [v,u], establish the Lemma. □ 
Not every order cone is normal. For example, consider the Sobolev spaces W'^'^ of 
scalar- valued functions on an n-dimensional, compact manifold A4, with Lipschitz 
continuous boundary, where fc is a non-negative integer, and p > 1 is a real number. 
An order cone in W'^'P is defined translating the order on the real numbers, almost 
everywhere in Ai, that is, 

:= {u e W'^'P -.u^O a.e. in M}. 

In the case fc = 0, that is, we have W'^'P — L^, the order cone above is a normal cone 
[3l[57]. However, in the case fc ^ 1 the cone above cannot be normal, since on the 
one hand, the cone definition involves information only of the values of u{x) and not 
of its derivatives; on the other hand, the norm in W'^'P contains information of both 
the values of u{x) and its derivatives. Since there is no boundary conditions on dA4 
in the definition of W'^'P, there is no way to relate the values of a function in Ai 
with the values of its derivatives. (In other words, there is no Poincare inequality 
for elements in W'^'P, with fc ^ 1.) 

An order cone X+ C X is generating iff Span(X_|_) = X. An order cone 
Ar+ C X is called total iff Span(X_|_) is dense in X. Total order cones are important 
because the order structure associated with them can be translated from the space 
X into its dual space X*. 

Lemma 17. Let X, Xj^ be an ordered Banach space. If X^ is a total order cone, 
then an order cone in X* is given by the set X^ C X* defined as 

xi:^{u* ex* ■.u*{v)^o VweA:+}. 

Proof. (Lemma\T^) We check the three properties in the definition of the order 
cone. The first property is satisfied because A"+ is an order cone, so there exists 
t; ^ in X+, and then there exists u* ^ in X* such that u*{v) = 1 > 0, so X'^ 
is non-empty. Trivially, G X'^. Finally, A"^ is closed because the order relation 
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^ for real numbers is used in its definition. The second property of an order cone 
is satisfied, because given any u* ^ v* G and any non-negative a, 6 G M, then for 
all u G X+ holds 

{au* + bv*){u) = au*{u) + bv*{u) ^ 
since each term is non- negative. This implies that (au* + bv*) G X^. The third 
property is satisfied because the order cone X^ is total. Suppose that the element 
u* e X^ and -u* G X^, then for aU u G X+ holds that u*{u) ^ and -u*{u) ^ 0, 
which implies that u*{u) = for all u G Therefore, u* G X^ C X* , where 

the super-script Ji in X^ means the Banach annihilator of the set X+, which is a 

subset of the space X*. Therefore, we conclude that u* G [Span(X+)]^. Since the 

order cone is total, Span(X+) = X, that implies [Span(X+)]^ = {0}, so u* — 0. 
This establishes the Lemma. □ 
An order cone X+ in a Banach space X is called a solid cone iff X+ has non- 
empty interior. The following result asserts that solid order are generating. We 
remark that the converse is not true. In the examples below we present function 
spaces frequently used in solving PDE with order cones having empty interior which 
are indeed generating. 

Lemma 18. Let X , he an order Banach space. If X^ is a solid cone, then X-^ 
is generating. 

Proof. (Lemma \1^ ) The cone Xj^ has a non-empty interior, so there exists 
xq G int(X+) and .tq 7^ 0. This means that given any x Q X there exists < a G M 
small enough such that both x+ := xq + ax and x- := xq — ax belong to int(X+). 
But then, x = {x+ —X-)/{2a), so a; G Span(X_|.). This establishes the Lemma. □ 
Here is a list of examples of several order cones used in function spaces. All these 
examples use order cones obtained from the usual order in M. In particular, they 
refer to scalar-valued functions on an n- dimensional, compact manifold M. with 
Lipschitz boundary. 

• Introduce on C'' the cone := {u € C'' : u{x) ^0 Va; G M}. This is 
an order cone for all non- negative integer k. The cone is a normal cone in 
the particular case k = 0. The cone is solid for all fc ^ 0, therefore it is a 
generating cone. 

• Introduce on L°° the cone :— {u G L°° : u ^ a.e. in A4}. This is a 
normal, order cone. It is a solid cone, therefore is generating. 

• Introduce on W''''^ the cone := {u G : u ^ a.e. in M}. 
This is an order cone. It is not normal for /c ^ 1. The cone is solid, therefore 
it is generating. 

• Introduce on L^ the cone := {u G L^ : u ^ a.e. in A4}. This is 
a normal, order cone every real numbers p ^ 1. The cone is not solid, 
however it is a generating cone. 

• Introduce on W'^'P the cone := {u G W'''^ : u ^ a.e. in M}. This 
is an order cone every real numbers p ^ 1. The cone is not normal for 
k ^ 1. The cone is not solid for hp ^ n, and it is solid for kp > n. In both 
cases, the cone is generating. 

A. 5. Maximum principles. We have not seen in the literature an approach to 
maximum principles on ordered Banach spaces in the generality we present it in 
this Section. Let X, X+ and Y, be ordered Banach spaces. An operator 
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A : Da C X ^ Y satisfies the maximum principle iff for every u, v € Da such 
that Au — Av G 1+ holds that u — vE X+. In the particular case that the operator 
A is linear, then it satisfies the maximum principle iff for all u G X such that 
Au G 1+ holds that u G X-^-. The main example is the Laplace operator acting 
on scalar-valued functions defined on different domains. It is shown later on in 
this Appendix that the inverse of an operator that satisfies the maximum principle 
is monotone increasing. The following result gives a simple sufficient condition 
for an operator to satisfy the maximum principle. This result is useful on weak 
formulations of PDE. 

Lemma 19. Let X . X^ be an ordered Banach space, and A : X X* be a linear 
and coercive map. Assume that X^ is a generating order cone, and that for all 
u E X such that Au G X*^_ there exists a decomposition u = u^ — u^ with u^ , 
u^ G Xj^ that also satisfies Au'^{u^) = 0. Then, the operator A satisfies the 
maximum principle. 

Proof. (Lemma\19i) Since the order cone X+ is generating, the space X* is also 
an ordered Banach space. Denote its order cone by X^. The assumption that the 
order cone X^ is generating also implies that for any element u G X there exists a 
decomposition u = u'^ — u~ with u"*", u~ G X^. By hypothesis, there exists at least 
one decomposition with the extra property that Au'^{u~) = 0. Now, by definition 
of the order in the space X* we have that 

Aue X\_ ^ Au{u) ^0 Vw G X+. 

Pick as test function u = u^ . Then, 

^ Au{u^) = - u^){u^) = - Au~{;u^) = -Au~{u~), 

where the last equality comes from the condition Am+(u^) = 0. Therefore, we have 

Au^{u^) ^0 = 0, 

because A is coercive. So we showed that u = u+ G X^. This establish the 
Lemma. □ 

An example is the weak form of the Laplace operator on scalar functions in the 
homogeneous Dirichlet problem on a compact manifold M with Lipschitz boundary. 
Consider the case X = Wo'^ with Y = X* = W'^^^, and X+ = W^:'^ , while 
y"_l_ = W^^'^ . The Laplace operator in this case is given hy A : X ^ X* with 
action Au{v) (Vu, Vw). It is not difficult to check that this operator satisfies the 
hypothesis in Lemma [T9l Therefore, this operator satisfies the maximum principle, 
that is, Au G W^^'^ implies u G W^'^, that is, u ^ a.e. in the manifold Ai. 
This result is in agreement with Theorem 8.1 in [25], where it is stated that: "If 
Au ^ 0, then infy^ u ^ — infg7\/i Here we introduced the cut-off function 

u^ := — min(u, 0) ^ 0. Recalling that in our case the domain of A contains only 
functions that vanish at the boundary, then inf^vi u ^ 0, that is, u in M. 

The following example is again the Laplace operator that appears in equations 
when they are written in weak form, but this time using more complicated operator 
domains due to more complicated boundary conditions in the PDE equation. Let 
{M, h) be a 3-dimensional Riemannian manifold, where is a smooth, compact 
manifold with a Lipschitz boundary dAA, and h G C^{M.,2) is a positive definite 
metric. Assume that the boundary set can be decomposed as follows, dM — 
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OMd U OMn and dMo H OMn — 0- Recall the definition of the Sobolev spaces 

W^o'^ {" e W^^'^{M,M) : tr^u = 0}, W^^^'^ := [W^^T- 

Then, define the operator 

Al : Wl^^ ^ Al<i>(^) := a40,0) + (A.12) 

where is the bilinear form 

: Wl^^ X ^R, a^C^^) := (V0, V^) + (i^ tr^v^, tr^0)^, 

and the Robin coefficient K G L°°{dJ^N, 0) satisfies the bounds 

k lltr^^ll^ < (iftr^0, tr^0)jv, V0 G W''^, (A.13) 

with k a positive constant. Assume that the function s G L^'^ , so the second term 
in the definition of the operator is well defined. 

Lemma 20. The operator A^ defined in Eq. iA.12\) satisfies the maximum prin- 
ciple. 

Proof. (Lemma\2(A) We now verify all the hypothesis in Lemma [T9l The cone 

"1^2 12 12 12 1 2 

Wj^j^ = W_|_' n Wjj is generating in W^' therefore, W^^ ' is also an ordered 
space. The constant k is positive and the function s is non-negative, which implies 
that the operator Al is coercive. Using the usual decomposition of a function u 
into u'^{x) = max7vi(u(x), 0) and u~{x) = — min7vi(u(x), 0), then it is not difficult 
to show that — 0, because the two parts of the decomposition of u 

are defined on non-intersecting parts of M. Therefore, Lemma 1201 follows from 
Lemma [191 D 

A. 6. Monotone operators. Let X, X+ and Y, y+ be two ordered Banach spaces. 
An operator F : X ~> Y is monotone increasing iff for all x, x G X such 
that X ~ X e X+ holds that F{x) - F{x) G 1+. An operator F : X ^ Y is 
monotone decreasing iff for all x, x € X such that x — x Q X+ holds that 
-[F{x) - F{x)] G Y+. The following result is a useful relation between linear, 
invertible operators that satisfy the maximum principle and monotone increasing 
operators. 

Lemma 21. Let X , X+ and Y , Yj^. be two ordered Banach spaces. Let A : X ^ Y 
he a linear, invertible operator satisfying the maximum principle. Then, the inverse 
operator A^^ : Y ^ X is monotone increasing. 

Proof. (Lemma\21\) Let y, y G F be such that y — y G 1+. Then, 

A{A-\y~v))^Y+ ^ A-\y-y)eX+ ^ A'^y - A-'ye X+. 

This establishes that the operator A^^ is monotone increasing. □ 
We are interested in a class of nonlinear problems where the principal part in- 
volves a linear operator A : X Y satisfying the maximum principle, and the 
non-principal part involves a nonlinear operator F : X ^ Y which has monotonic- 
ity properties; problems of this type can be written as follows: Find an element 
X G X solution of the equation 

Ax + F{x)^0. (A.14) 
We now establish some results for this class of problems. 
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Lemma 22. Let X , and Y , 1+ be two ordered Banach spaces. Let A : X ^ Y 
be a linear, invertible operator satisfying the maximum principle. Let F : X ^ Y 
be a monotone decreasing (increasing) operator. Then, the operator T : X ^ X 
given by T := —A~^F is monotone increasing (decreasing). 

Proof. (Lemma\^) Assume first that the operator F is monotone decreasing. 
So, given any x, x e X such that x — x e the following inequalities hold, 

x-xeX+ => -[F{x)-F{x)\eY+, 

^ A{-A-^[F{x)-F{x)\)eY+, 

^ -A-\F{x)-F(x)\^X+, 

^ -[A-^F{x)- A-^F{x)\eX+, 

^ T{x) ~T{x) e X+, 

which establishes that the operator T is monotone increasing. In the case that the 
operator F is monotone increasing, then the first line in the proof above changes 
into X ~ X E X^ implies that F{x) — F{x) £ and then all the remaining 
inequalities in the proof above are reverted. This establishes the Lemma. □ 
The next result translates the inequalities that satisfy sub- and super-solutions 
to the equation Ax + F{x) — 0, into inequalities for the operator T = —A^^F. 

Lemma 23. Assume the hypothesis in Lemma l^M 

If there exists an element x-^- € X such that Ax^+F(x-f-) G 1+, then this element 
satisfies that a;_|_ — r(a;+) G X^. 

If there exists an element x^ £ X such that —[Ax^ + i<"(x_)] G 1+, then this 
element satisfies that — [a;_ — T{x^ j^ G X^. 

Proof. (Lemma\23[) The first statement in the Lemma can be shown as follows, 

Ax+ + F{x+) G r+ ^ A{x+ + A-^F{x+)) e Y+ 

=^ x+ +A^'^F(x+) G X+, 

which then establishes that x+ — r(x+) G X+. In a similar way, the second state- 
ment in the Lemma can be shown as follows, 

-[Ax- + F{x-)] G r+ ^ A{-x- -A-^F{x-)) e Y+ 

-A-^F{x-) G X+, 

which then establishes that — [x- — T(a;_)] G X+. This establishes the Lemma. □ 
The last result can be found as Theorem 7. A in [57j, page 283, and Corollary 7.18 
on page 284. We reproduce it here for completeness, without the proof. 

Theorem 15. (Fixed point for increasing operators) Let X be an ordered 
Banach space, with a normal order cone Xj^. Let T : [x-,x^] d X ^ X be a 
monotone increasing , compact map. //— — T(a;_)] G X+ and x-^- — T{x-^-) £ X^, 
then the iterations 

Xn+i := Tixn), xq = a;_, 

Xn+l ■= T{Xn), Xq = X+, 

converge to x and x G respectively, and the following estimate holds, 

X- ^ Xn ^ X ^ X ^ in ^ x+, V?i = N. (A. 15) 
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For nonlinear problems of the form (jA.14p . one can use Theorem[T5]for monotone 
nonlincarities to conclude the following. 

Corollary 3. (Semi-linear equations with sub-/super-solutions) Let X , X+ 

and Y , be two ordered Banach spaces where is a normal order cone. Let 
A : X ^ Y be a linear, invertible operator satisfying the maximum principle. Let 
x+, E X be elements such that (x+ — a;_) G X^, and then assume that the 
operator F : C X —>■ Y is monotone decreasing and compact. If the 

elements x^ and X-^- satisfy the relations 

- [At_ + F{x-)\ e y+, Ax+ + F{x+) e y+, (A.ie) 

then there exists a solution x G [x-, x+] C X of the equation Ax + F(x) — 0. 

Proof. ( Corollary O ) The operator A is invertible, then rewrite the equation 
Ax + F(x) = as a fixed-point equation, 

X = ^A-^F{x) T{x). (A.17) 

By Lemmaim we know that the map T : X ^ X is monotone increasing. Moreover, 
this operator T it is compact, since is the composition of the continuous mapping 
—A~^ and the compact map F. The elements X- and x+ satisfy Eq. (|A.16|) . 
therefore, by Lemma [231 they are also sub- and super-solutions for the fixed-point 
equation involving the map T. It follows from Theorem [15] that there exists an 
element x € X solution to the fixed-point equation (|A.17p , and this solution satisfies 
the bounds a;_ ^ x x+. □ 

A. 7. A priori estimates in ordered Banach spaces. Many problems of the 
form in Eq. (jA.14[l do not have monotone nonlincarities. However, in the case that 
there exist sub- and super-solutions to Eq. (|A.14p it is possible to introduce a "shift" 
into the equation. This shift transforms a problem that does not have a monotone 
nonlinearity into one that does, without destroying the maximum principle property 
required of the linear part. However, the disadvantage of the shift technique is 
that it requires additional regularity in the equation coefficients than the regularity 
needed for the original equation to be well-defined. On the other hand, it is possible 
to construct arguments leading to a priori order cone estimates on any possible 
solution (whether or not it exists) with very weak assumptions on the nonlinearity. 
Although such results are standard for semi-linear scalar problems with monotone 
nonlincarities (for example, see [32]), our result below holds for a class of semi-linear 
problems with non-monotone nonlinearities and appears to be new. Problems with 
monotone nonlinearities fit into this class, but it also includes a much larger set of 
nonlinearities. (See the second assumption in[i|in Lemma [24] below. ) 

The following result (Lemma[24|below) gives sufficient conditions for establishing 
a priori order cone estimates on solutions to certain PDE-like operator equations 
in ordered Banach spaces. These order estimates can be translated into norm 
estimates in the case that the order cone is normal. (See Corollary [3| following 
Lemma [Ml below.) Note that the bounds established in Lemma [Ml below are not 
necessarily sub- and super-solutions; establishing the bounds by first showing they 
are sub- and super-solutions and then using Corollary [3l would require a mono- 
tone nonlinearity, or use of the shifting technique requiring additional regularity 
assumptions. 
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Lemma 24. {A priori order estimates) Let Y, 1+ be an ordered Banach space 
with a generating order cone Kf. Let F : Y ^ Y* be a continuous map. Let 
A : Y —* Y* be a linear, continuous operator with dimiVyi ^ 1. Assume that there 
exists a subspace X C Y , with an induced order cone — X D 1+, such that 
Ax : X ^ X* , the restriction of the operator A to the space X , is coercive. Let 
u Cz Y be a solution of the equation Au + Fiu) — 0. 

(i) If there exists an element y/y £ Na such that (u — j/a)^ G and for all 
y dz Y such that (y — j/a) G Y-^ holds that F{y)(^{y — J/a)^) ^ 0; Then, the 
solution u dz Y satisfies ^ u £ 1+ . 

(ii) If there exists an element y^ £ Na such that (it — j/v)^ G o,nd for all 
y &Y such that —{y — ?/v) G ^h- holds that F{y)(^{y — 2/v) ) =^ 0. Then, the 
solution u £ Y satisfies u — y^ £ Y-^ . 

Proof. (Lemma\24\) We first show part (Op. Given the solution u £Y , introduce 
and element Uo G 1" be an element such that u — Uo £ X. Second, notice that the 
element {u — y/C) belongs to the space Y , which has a generating order cone F+, so 
we know that there exists a decomposition 

{u - y^) = {u~ y^)+ - {u - yf,Y, 

with both elements (u — y^)^, — y/^)^ £ Y^. The first assumption in ^ says 
that (u — t/a)^ G X^ and so the element (u — t/a)^ is a valid test function for the 
functional 

[Ax{u - ua) + F{u) + Aua] £ X*, 
so we have the following, 

Au{u - yA)^ = - yA){u - yA)^ 
= A(u - yA)^('u - ^a)^ 
= Ax{u - yA)^(w - 

Therefore, we have the following inequalities, 

= Au{u - yA)+ + F{u){{u ~ y^)+) 
= Axiu- y^)+{u - y^)+ + Fiu){iu - y^)+) 
^ ao\\{u-y^)+\\l+F{u){{u-y^)+), oq > 
^ flo - yA)^||x- 

The last inequality implies that {u — yA)^ = 0, which then says that u — y^ ^ 
— {u — y;y)~ , and we then conclude that ^{u ^ y^) £ Y^. This condition can be 
written using inequalities as w ^ ^/a- 

We now prove part The element {u — y^) also belongs to the space Y, which 
has a generating order cone 1+, so we know that there exists a decomposition 

(u - yv) = (w - yv)'^ - {u~ yv)", 

with both elements (m — yv)^, {u ~ y^/)^ £ Y^. The first assumption in part fzj) 
says that {u — y-^)^ G X+ and so the element {u — y\j)^ is a valid test function for 
the functional 

[Ax{u - ud) + F{u) + Aud] £ X*, 
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SO we have the fohowing, 

Au{u - y^y ^ A{u - 2/v)(w - 

= -A{u - yvy{u - ?/v)~ 
= ~Ax{u - y^yiu - 2/v)~- 
Therefore, we have the following inequalities, 

= Au{u - 2/v)" + F{u){{u - yv)") 
= -Ax{u - y^yiu - y^y + F{u){{u - yv)~) 
^ -ao \\{u ~ y-^y\\l + F{u){{u - y-^)-) , ao > 
< -ao - yv)"||x- 

The last inequality implies that {u ~ y^y = 0, which then says that u — y\, = 
{u — y^y , and we then conclude that [u^yy] e y+. This condition can be written 
using inequalities as u ^ y v This inequality establishes the Lemma. □ 

Corollary 4. Let Y , K|_ be an ordered Banach space with a normal order cone 
y+. Let Z be a Banach space, and consider the space W = Y Ci Z , with order cone 
W+ :~ Y-i-ClZ . Lf there exist elements u, y^/, andy^ g W such that ^ yy ^ u ^ y^ 
in the order given by then there exists a positive constant c such that the 

following inequalities hold 

c\\yAzi^\\u\\zi^ \ \\yA\z- 

Proof. (Corollary^) It follows directly from the definition of a normal order 
cone. □ 
As an example. Lemma [24] holds with the spaces taken to be X W^'^, Y ~ 
W^'^, the linear operator taken to be Au(v) — (Vu, Vu), and the nonlinear operator 
taken to be a monotone operator such as F{u) = . Lemma [24] also holds for a 
non-monotone nonlinear operator satisfying the assumptions for the Lemma, such 
as F{u) — ~ 2u^. An example of the space Z where Corollary [4] holds is Z = L°°. 
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